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1. Hypothèse du chaos initial. — La théorie 
moléculaire des phénomènes irréversibles qui se 
manifestent dans les transformations macrosco- 
piques dès qu’elles ne sont pas infiniment lentes est 
actuellement l’objet de nombreuses études. Le 
domaine le mieux exploré de cette théorie est 
dominé par l’équation de Boltzmann. Cette équa- 
tion a donné lieu à d'innombrables applications : 
c’est une matière qui, même aujourd’hui, est loin 
d’être épuisée. 

L’équation de Boltzmann ne concerne que les gaz 
dans un état de densité modérée. Sa déduction sup- 
pose que les forces intermoléculaires sont à courte 
portée. Elle ne tient compte en conséquence que des 
rencontres de molécules prises deux à deux. En 
restant fidèle à l’hypothèse des forces à courte 
portée, il n’a pas été possible d’aller beaucoup plus 
loin que d’y introduire tant bien que mal l’effet des 
rencontres triples. L'extension aux forces à longue 


portée — c’est principalement le problème des 
plasmas très ionisés — n’a pas reçu de solution 
irréprochable. 


Ce succès limité des tentatives de généralisation, 
de l’équation de Boltzmann tient au fait que les 
difficultés de calcul s’accumulent avec les difficultés 
physiques au point de rendre les raisonnements 
assez obscurs. C’est que l’équation de Boltzmann 
elle-même demande des éclaircissements. 

Du point de vue des principes, le problème est de 
concilier l’équation de Boltzmann et les équations 
élémentaires de la mécanique, c’est-à-dire une 
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équation relativement microscopique et entachée 
d’irréversibilité d’une part et des équations fonciè- 
rement microscopiques et tout à fait exemptes 
d’irréversibilité d’autre part. L'utilisation comme 
intermédiaire de l’équation de Liouville donne une 
meilleure délimitation de ce problème. Il n’en est 
pas résolu pour autant. 

L’équation de Boltzmann est une équation de la 
mécanique classique. La mécanique quantique a 
apporté de nouvelles lumières à la mécanique statis- 
tique, ceci essentiellement parce que la mécanique 
quantique, étant probabiliste par essence, s’intègre 
sans hypothèse particulière la mécanique statis- 
tique. Il est donc naturel de rechercher la conei- 
liation de la réversibilité mécanique et de l’irréver- 
sibilité macroscopique dans le cadre de la méca- 
nique quantique. 

A ces raisons de principe il s’ajoute un avantage 
pratique. L’immense effort de développement de la 
mécanique quantique depuis trente ans a doté 
celle-ci de puissants moyens de calculs. Il n’est pas 
douteux, puisque la mécanique classique est une 
approximation de la mécanique quantique, que tout 
moyen qui est bon pour la seconde est aussi bon 
pour la première et que ces progrès peuvent être 
mis à profit dans le domaine de la mécanique sta- 
tistique classique et notamment dans l’étude de 
l’irréversibilité. Mais cet effort de transposition des 
méthodes n’est pas indispensable. Ce qui ne veut 
pas dire que la mécanique quantique seule est sus- 
ceptible d’engendrer l’irréversibilité. 
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Les systèmes mécaniques auxquels nous pensons 
ici ont un grand nombre de degrés de liberté. Leur 
description classique s’effectue à l’aide d’une fonc- 
tion de toutes les coordonnées spatiales, de toutes 
les quantités de mouvement et du temps : c’est la 
densité en phase. La description quantique s’effec- 
tue à l’aide de la matrice densité : celle-ci comporte 
dans la présentation spatiale, le même nombre de 
variables que la densité en phase, mais ce ne sont 
pas les mêmes. La matrice densité dépend en effet 
de deux jeux de coordonnées spatiales et du temps. 
Il est possible, si on le désire, d'adopter les variables 
indépendantes de la théorie classique dans la théo- 
rie quantique. Mais cette modification, qui permet 
de rapprocher les deux points de vue et de procéder 
à l’écriture de l’approximation classique lorsque le 
formalisme quantique est connu, engendre d’inu- 
tiles complications pour le traitement quantique 
lui-même. 

Par des raisonnements, qu’aujourd’hui nous pou- 
vons nous permettre de regarder comme intuitifs, 
Pauli a jeté les fondements d’une théorie quantique 
des phénomènes irréversibles. Prenant comme 1in- 
connue ce qu’on appelle aujourd’hui la population 
des états énergétiques d’un système, il a formulé 
l’équation d’évolution de ces populations. Cette 
équation permet de montrer qu’un système hors 
d’équilibre tend vers l’équilibre thermodynamique. 
Écrite dans un contexte très différent, elle n’est pas 
sans analogie avec l’équation de Boltzmann. 

L’équation de Pauli peut s’énoncer dans le lan- 
gage des matrices densités. La matrice densité doit 
être exprimée non plus à l’aide des coordonnées 
continues, mais à l’aide du système orthonormé des 
états énergétiques du système. Dans cette représen- 
tation, les éléments diagonaux de la matrice densité 
ne sont pas autre chose que les populations dont il 
vient d’être question. Quant aux éléments non 
diagonaux, on s’en débarrasse par des artifices : 
ils ne figurent pas dans l’équation de Pauli. 

Les progrès dans la voie ouverte par Pauli ont 
été tentés dans deux directions différentes : la 
première a consisté, dans des recherches d’un carac- 
tère général, à trouver une justification des arti- 
fices précités. C’est dans cette ligne qu’il faut situer 
les recherches de Vau Hove [1]. 

La seconde a été de moindre ambition, plus utili- 
taire. Les phénomènes de la résonance magnétique, 
soit électronique, soit nucléaire, sont typiquement 
quantiques. Lorsqu'ils sont observés dans un milieu 
macroscopique, un cristal, un liquide, certains 
effets irréversibles, la tendance du système à 
atteindre l’équilibre thermodynamique se mani- 
feste avec éclat. Les spécialistes de ces ques- 
tions [2] ont donc été conduits pour leurs propres 
besoins à imaginer la théorie des phénomènes qu’ils 
étudiaient. Une de leurs conclusions, qui a ici une 
importance majeure, est qu’il n’est pas question 
d'effacer, fût-ce par des raisonnements ingénieux, 
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les éléments non diagonaux :ils sont indispensables 
à la description des phénomènes. Du point de vue 
de la résonance magnétique, le problème n’est pas 
de justifier l'équation de Pauli. Il s’agit de trouver 
une équation plus large. 

Il semble que cette attitude doive gouverner 


aussi les recherches de portée générale. Il ne faut 


pas rechercher une équation d’évolution des popu- 
lations. Il faut rechercher une équation d’évolution 
de la matrice densité. Que dans certains cas les 
éléments non diagonaux s’annulent ou soient sans 
importance, ce point ne saurait faire de doute étant 
donné le succès de l’équation de Pauli en maintes 
applications. Ii n’en reste pas moins qu’il faut voir 
la question sous un angle différent. 

Les recherches pour fonder l’équation de Pauli 
reposent généralement sur l’arbitraire que comporte 
la définition des états propres d’un système : ces 
états propres ont chacun une phase arbitraire. Il 
en résulte que, lorsqu’on prend ces états pour base 
d’une certaine représentation d’une matrice, les 
éléments diagonaux sont bien définis cependant 
que les éléments non diagonaux ne le sont pas. Une 
moyenne prise sur ces phases arbitraires annule les 
éléments non diagonaux. Cette propriété est la 
source des raisonnements qui ont pour objet de 
minimiser le rôle des éléments non diagonaux. 
Notre opinion sur le rôle des éléments non diago- 
naux nous conduit à imaginer tout autre chose. 

Après l’exemple de la résonance magnétique, 1l 
est possible d’en citer un autre : Prigogine et 
Résibois [3] ont examiné comment un gaz quan- 
tique approche de l’équilibre. Ils utilisent à propos 
de chaque degré de liberté les états propres indi- 
viduels de chaque particule tels qu’ils sont définis 
lorsqu'on néglige toute interaction. Ils aboutissent 
à une équation du type de Pauli. Ce qui veut dire 
qu’ils se sont débarrassés des éléments non dia- 
gonaux dont nous avons parlé plus haut : ce résul- 
tat a pu être atteint parce que ces auteurs ont 
traité seulement le cas d’un fluide homogène, ils 
ont traité un problème particulier de l’approche à 
léquilibre que l’on pourrait appeler le problème 
de la « thermalisation » — s’il s’agissait d’un pro- 
blème classique on dirait qu’ils ont montré comment 
un gaz homogène où la répartition des vitesses est 
d’abord quelconque acquiert peu à peu une répar- 
tition maxwellienne. Comme Prigogine et Résibois 
le remarquent eux-mêmes dans une courte note 
(page 814) leur méthode ne couvre pas les phéno- 
mènes de transport — avec leur escorte, viscosité, 
conductibilité thermique, etc. Cette remarque 
renforce l’idée que le problème général n’est pas de 
justifier directement l’équation de Pauli. 

Je me propose d'explorer ici avec quelques détails 
une hypothèse que j'ai déjà utilisée dans plusieurs 
publications [4] et qui a été également préconisée 
par Tchen à propos des plasmas. 

Cette hypothèse comporte deux parties : 


VE a 
7 Fr 
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19 Tout système s’est trouvé constitué à un 


. moment donné par le rapprochement de ses élé- 


ments — éléments qui ont pu être soit micros- 
copiques soit macroscopiques. 

20 Avant le rapprochement, les éléments qui 
étaient séparés n’avaient aucune corrélation entre 
eux. C’est l'hypothèse du chaos initial. 

I faut bien dire que la première hypothèse risque 
de nous obliger à remonter extrêmement loin dans 
le passé du système. La seconde est analogue à 
l'hypothèse du chaos moléculaire posée par Boltz- 
mann dans la déduction de son équation, mais ici 
elle n’est demandée que pour l'instant initial. Notre 
hypothèse du chaos initial ne demande pas d’autre 
justification que le succès de ses conséquences, alors 
que l’hypothèse du chaos moléculaire doit être 
continuellement renouvelée à chaque instant. La 
justification de l’hypothèse du chaos moléculaire, 
dans la ligne de l'hypothèse du chaos initial, réduira 
si On y parvient, cette première hypothèse à n’être 
qu’une approxhnation. 

Ce qu'il faut souligner aussi c’est que nous 
n’essaierons nullement de décrire la formation du 
système, nous contentant seulement de l’idée utili- 
taire qu’à un instant très lointain dans le passé il 
n’y avait pas de corrélation. Pratiquement notre 
hypothèse sera donc moins ambitieuse que les deux 
propositions précédentes, mais aussi moins claire. 

L'hypothèse du chaos initial est une hypothèse 
purement statistique. Van Hove (loc. cit., p. 536) 
a souligné que seule une hypothèse sur la structure 
de l’état initial peut concilier la réversibilité des 
équations dynamiques avec l’irréversibilité des 
phénomènes macroscopiques naturels. Il est clair 
que notre hypothèse satisfait à cette prescription. 

Nous nous proposons d’appliquer ici nos hypo- 
thèses au cas de deux systèmes en faible inter- 
action. Ceci afin de réduire les difficultés mathé- 
matiques. L’inconvénient toutefois est le suivant : 
la division d’un système en deux parties est physi- 
quement assez évidente dans bien des cas, mais s’il 
s’agit de décrire un problème de caractère général, 
elle est à première vue purement formelle. 


2. Les corrélations et l’évolution. — Un système 
isolé, C, est composé de deux parties À et B. 
L’hamiltonien de À est désigné par AT ,, celui de B 
par ÀH,. L’interaction entre À et B est représentée 
par l’hamiltonien #V. Le système global C a donc 
pour hamiltonien : 


4H = K{H1 + Hg + V). 


Pour le moment nous n’imposons pas à ces 
hamiltoniens d’être indépendants du temps. 

La situation du système C à un instant donné est 
représentée par une matrice densité appropriée. Au 
point de vue du langage, le terme de matrice 
évoque un certain type de représentation. Nous pré- 
férons souvent parler d’opérateur densité plutôt 
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que de matrice densité, adoptant ainsi un point de 
vue plus abstrait. Nous désignerons par up cet 
opérateur. Nous posons en général, X et Y étant 
deux opérateurs quelconques : 


ENT AE NX. 


Ceci dit, l’évolution de y est écrite par l’équation 
de Liouville-Von Neumann : 


Du = [A ul (1) 


Si on ne s'intéresse qu’au système À, l’infor- 
mation fournie par y est trop généreuse. Une ré- 
gression permet d'obtenir l’opérateur densité rela- 
tif au système À seul : 


ua = TB. 


L'opérateur Tz consiste à prendre une trace — 
partielle — à travers le système B. De même, le 
système B est entièrement situé par la connais- 
sance de : 

uB = Ta. 


Notons que pour le moment nous ne précisions 
nullement s’il s’agit de systèmes macroscopiques. 

L'évolution de u, et celle de u3 s’obtiennent en 
prenant membre à membre des traces de l’équa- 
tion (1) : 


uA Te ei (2) 


Su = (Mn, us] + TAUV, v] 


Ces équations se sufliraient à elles-mêmes si 
l'interaction V était nulle. 
Introduisons les notations suivantes : 


JE = + H3, 


t 
U exp.— i [ H dt, 


I 


Uo 


I 


D 
EXP. — il H, dt. 


L'équation (1) permet d’exprimer l’opérateur 
densité u à l’instant { en fonction de sa valeur 
initiale : 

u(s) = Uu(0) UT. (3) 

Réciproquement : 

u(0) = UT (0) U. 

Un calcul de traces donne : 

u4(0) = Tr UT ut) U, 
ug(0) = Ta Ut ut) U. 

Écrivons qu’il n’y avait aucune corrélation entre 

les deux systèmes à l'instant imitial : 


u(0) — 14(0) uB(0). (4) 
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On notera qu’une telle relation n’est possible que 
si les deux systèmes ne sont pas susceptibles 
d’échanger des particules de même nature, plus ou 
moins nombreuses, car au cas contraire elle ne 
satisferait pas aux conditions de symétrie qu’im- 
posent bosons ou fermions. 

Maintenant on élimine (0), u4(0), ur(0). 

Il vient : 


u(t) = U(Ts Ut us) U) (Ta Ut po U) UT. (5) 
Il est commode d'introduire l’opérateur unitaire : 
Q(:) = U(r) Vo. (6) 
Son inverse est : 
DORE C)AUIT) 
On élimine alors l’opérateur À de l’expression (5). 
Le jeu des produits ou les possibilités de permu- 


tation dans les calculs de trace fait disparaître éga- 
lement /,. Il reste : 


ue) = Q(+) (Ts Q{) us) Q(1) (T4 Qu) A()) A()T. 
(7) 


Introduisons maintenant l’hypothèse suivant 
laquelle interaction V est faible, 
Q(+) diffère peu de l’unité. On posera : 


Q(t) = 1 — iL{), (8) 
L étant un opérateur hermitique du premier ordre 


dont l’expression est donnée par la théorie des per- 
turbations : 


Lt = JF U,(8) VU, (6): d6. 


Au premier ordre, l’équation (7) se simplifie : 
u = au — iL(t), ua ua] 
+ us TB[L(t), u] + ia TalL{(t), u]. (9) 
Tous les opérateurs densités étant considérés au 
même instant {, on a négligé et on négligera désor- 


mais d’expliciter leur dépendance du temps. La 
différence 


H — HA UB 


est du premier ordre. Il est légitime de remplacer 
dans les deux derniers termes de (9) l'opérateur u 
par le produit 44 Ur, Ce qui donne : 


u = Lau — ilL(), ua us] 
+ 18 TB[L(t), 44 8] + iua TalL(t), ua ul. 
La formule (10) exprime en fonction de l’état de 
chacun des systèmes — au premier ordre — ce que 
sont devenues les corrélations initialement nulles à 


l'instant t. 
Comme 


(10) 


T'ag[L(t), u4 us] — 0 en raison des propriétés des traces 


Due par normalisation 


la trace partielle membre à membre de l’équa- 
tion (10) sur B donne comme 1l convient : 


TAN = 


Détaillons les traces qui figurent dans la for- 
mule (10). Par exemple : 


Tell, ta us] = Tell va uB — ba us L] 
= (Tp Lu) ua — UA(T8 L us). 


On posera : 
HA VB EM 
TB Lu = L: 
Ta Lua = Le 
L'= L— Li — Lz. 
Ilvient donc:  . 


u = v— 1[1}, v]. (12) 

Introduisons l’expression (12) de x dans l’équa- 
tion (2) d'évolution de 4. Nous obtenons, au deu- 
xième ordre, l’équation suivante : 


D ua — [H4 + Va, ua] 
STEP LINE LV NL Co 


On a posé : 
PUR: 


L’équation (13) a pour conséquence, comme il 
convient, la conservation de la trace de pu, : 


£ T4 LARA0 

Une équation analogue à (13) peut s’écrire pour 
le système B. L'ensemble de ces deux équations 
constitue un système autonome et complet pour 
décrire l’évolution du système à l’approximation 
envisagée. 

Étant donné que nous n’avons rien fait d’autre 
que d’introduire une hypothèse sur l’état initial, 
u défini par (12) vérifie de son côté l’équation fon- 
damentale (1), comme il est aisé de s’en assurer par 
simple substitution. Toutefois cette vérification n’a 
lieu qu’au premier ordre. 

Il est possible de présenter de diverses manières 
l'équation d'évolution, ces présentations étant plus 
ou moins commodes suivant les circonstances. 

19 On fait apparaître explicitement L, et La. 
On s’aperçoit alors que le terme en Z; peut s’incor- 
porer dans l’hamiltonien : 


0 : 
ua = [Ma + Va + iT[V, Le] ur, ua] 


— iTz{ VLua 8 + ua uB LV — Vus us L — Lu ur V } 
A if Va La pa + ua La Va — Vaua La — La ua VA }. (14) 


N°6. 
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20 On peut grouper autrement ces deux derniers 
termes de manière à faire apparaître les différences : 


TRI ERL 
ÉENDEEN TE 


Il vient : 


DD ; 
TRS (Ha + Va € iTB[V, La] ur, LA] 
es iT B(V"L'u4 UB :e LA UB L" p” 


er V'u4 UB DZ 25 ÿL2 UA UB 17), (15) 


Cette formule montre qu’une partie importante 
de l’interaction entre les deux systèmes est gou- 
vernée par les «fluctuations » V” et L”. 

30 Les présentations qui précèdent ne diffèrent 
que formellement. Maintenant nous observerons 
qu'un problème de perturbation n’est pas exac- 
tement défini. La décomposition de l’hamiltonien 
du système global en trois éléments 


HAE EC. 


dont le dernier est du premier ordre, peut être rem- 
placée par une autre décomposition où de petites 
quantités sont ajoutées aux hamiltoniens princi- 
paux À, et H3, mais retranchées de la pertur- 
bation, ce qui conduit à une nouvelle décompo- 
sition du type suivant : 


H1 + K1 + 


Hp + KB + (V K 1 K3B). (16) 

3. Corrélations limites. — Nous nous limitons 
maintenant au cas où les hamiltoniens ne dé- 
pendent pas du temps. Nous désignerons les états 
propres du système À par a, b, c..., ceux du sys- 
tème B par «, GB, y, … Cette notation peu symé- 
trique est commode parce que par la suite les sys- 
tèmes À et B jouent des rôles très différents. 
L'élément de matrice de L s’écrit : 


t ; 
Lox,ve(t) = e—i(Hat H4)0 Vaud ei(Hp+Hp)0, do (1) 
0 
Si l’un des systèmes au moins a un spectre con- 


tinu, cet élément de matrice a une limite lorsque t 
devient très grand. Cette limite s’écrit : 


- Vou,ve 
D ie Hat 


L &ÿ 
er He = He 


+ (Ho + He— Ha — Ha)|] Vote: (2) 


Il suffit de reculer suffisamment l’origine des 
temps pour avoir le droit désormais de remplacer 
dans les équations du paragraphe précédent Z(t) 
par L(co) que nous désignerons dans la suite par L 
tout court. En même temps il faut nécessairement 
reculer indéfiniment dans le passé l’application de 
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l'hypothèse de l'absence de corrélations initiales, 
ce que l’on peut faire sans regret. 

Maintenant la formule (10) ou (12) donne une 
expression instantanée des corrélations. L’évolu- 
tion du système est entièrement définie par la 
situation présente de ses deux composantes. 

On notera que le formalisme que nous avons 
obtenu n’est rien d’autre qu’un aspect élaboré de la 
méthode de perturbation appliquée à la méthode 
de variation des constantes. 

Proposons-nous maintenant de calculer les corré- 
lations qui existent entre deux systèmes légè- 
rement couplés qui sont en équilibre thermody- 
namique. Nous ferons les évaluations sur la répar- 
tition canonique à la température 7. Nous pose- 
TONS : 

SRE 


Pour le reste nous emploierons autant que pos- 
sible les mêmes notations. C étant la constante de 
normalisation, l’opérateur densité du système 
s’écrit : 

u = C'exp [— s(H4 + Hp + V)]. 


Les opérateurs u, et u, s’en déduisent par les 
formules régressives habituelles. Nous désignerons 
respectivement leurs éléments de matrice par : 


Ua,b OÙ Hop 


A l’ordre zéro, , u4 et uz sont diagonaux. Des 
calculs élémentaires de perturbation donnent le 
résultat suivant, valable au premier ordre : 


tee Ubb URB — aa Ho Y 
Uau,bB — Ua,b Ua, SE 11 ne Hg sgrfet 1 ax,b8 


Hbb — aa 
af TT Æ 1, à Ho Va. ba 


0 


ns , HBB — Haa 
Du, CVs EE EEE V. 
Haa Oab Hg ie Haa V'aux,añ 


—— SUHaa ab Uouc da » Haa Hu Vax,ax (3) 
a,æ 


Lorsqu'on applique formellement la formule (10) 
au cas de l’équilibre précédent, en tenant compte 
des expressions approchées de pu, et de pu : 

Uab  Jab EXP (— SHa), 
Ua,p dog EXP (— So), 


les termes qui dérivent de 
SA TH) 


disparaissent. Le symbole S devient inutile et les 
deux expressions prennent exactement le même 
aspect, à ceci près que l'expression purement ther- 
modynamique (3)contientleterme supplémentaire: 


— SULaa Ôab Ho da 2 Haa oc Vau,au (4) 
Ne4 


Il existe une autre différence plus profonde que 
dissimule l’uniformité des notations. Dans le cas 
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des régimes variables nous avons supposé que l’un 
des systèmes avait un spectre continu. Dans le cas 
de l’équilibre thermodynamique, le système est 
fermé : il n’a que des états discrets. 

Une première condition est nécessaire pour que 
le schéma de l’équilibre rentre dans le schéma du 


régime variable : c’est que l’une au moins des com- 


posantes du système À + B soit extrêmement 
grande de manière qu’elle ait pratiquement un 
spectre continu. Une seconde condition s’écrit : 

D 


24 Uaa Ua Vaaax = 0. 
a,œ 


Ce qu’on obtiendra par exemple par une décom- 
position appropriée de l’hamiltonien suivant (2, 16). 

Dans la suite nous admettrons que ces deux 
conditions sont remplies, il en résulte qu’alors les 
corrélations obéissent à la même relation à l’équi- 
libre et hors d’équilibre. 


4. Équilibre en présence d’un thermostat. — 
Nous nous proposons d’appliquer les résultats pré- 
cédents au problème de la mise en équilibre ther- 
modynamique de deux systèmes légèrement cou- 
plés. 

Ils ont été mis en contact il y a très longtemps. 
A ce moment il n’existait aucune corrélation entre 
eux. L'opérateur L a évolué assez rapidement et a 
atteint sa valeur limite. Les corrélations ont donc 
atteint l’état d’organisation instantanée prédite 
par la formule (2,12). Mais les opérateurs densi- 
tés 1, et # continuent à évoluer jusqu’à l’équilibre 
total. Du moins l’expérience nous enseigne que 
l’ensemble À + B tend vers l’équilibre si l’un au 
moins des deux systèmes — disons B pour fixer les 
idées — est un système macroscopique. Ce que 
nous supposerons et ce qui satisfait à l’hypothèse 
qu’un des systèmes a un spectre quasi continu. 

Une sérieuse limitation subsiste. Prenons un 
exemple : le système B est une bassine d’eau 
chaude, le système À est une bouteille de lait froid. 
Une ménagère met la bouteille de lait dans la 
bassine, de manière à réchauffer ce liquide. Confor- 
mément à notre description où aucun opérateur 
n'intervient à partir de l’instant initial, elle laisse 
tomber la bouteille de lait dans la bassine. Il est 
clair que cet exemple ne correspond pas très bien à 
l’idée que l’on se fait d’une interaction faible, mais 
ce point est secondaire pour notre argumentation. 
Suivons les événements : on suppose pour sim- 
plifier que la bouteille de lait est tombée d’aplomb 
dans l’eau ! L’eau d’abord mise en mouvement, 
revient assez vite au repos puis le lait s’échauffe 
peu à peu. Les mouvements de l’eau se ralentissent 
non pas parce que l’eau est couplée au système À, 
mais par ce que l’eau est douée de viscosité. Nous 
ne pouvons donc pas prétendre — en général — 
décrire la mise en équilibre de deux systèmes com- 
plexes’par le seul jeu de leurs interactions. Il faut 
attribuer à l’un et à l’autre leur irréversibilité 
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propre. Ceci nous conduit — en admettant bien 
entendu que nous sommes dans la bonne voie — à 
subdiviser À et B en parties de plus en plus petites 
et finalement à descendre jusqu'aux degrés de 
liberté élémentaires. 

On échappe à cette difficulté si aucun des deux 
systèmes, pris isolément, n’a pas le moyen d’être 
irréversible : par exemple le système À est un élec- 
tron non relativiste et le système B est le champ de 
rayonnement. Bien entendu le problème du pas- 
sage à la limite subsisterait et nécessiterait de toute 
manière la considération des degrés de liberté 
élémentaires. 

Nous pouvons heureusement tester en partie nos 
idées sans renoncer à la simplicité en traitant le 
problème suivant : 

Le système À est un système de faible capacité 
énergétique. Il est relativement simple et n’a 
aucune raison de se thermaliser spontanément, Le 
système Best un thermostat. Il est, lui, toujours si 
proche de l’équilibre qu’on ne s’occupe pas de 
savoir s’il s’en écarte ou s’il en revient. Ces pertur- 
bations, dans la mesure où elles existent, se dis- 
persent dans le thermostat et ne reviennent pas 
troubler l’évolution du système A. Ou encore le 
système À est assez délié pour se déplacer dans le 
thermostat plus vite que les perturbations qu’il 
provoque. Le premier cas est celui du moment d’un 
ion paramagnétique placé dans un milieu pourvu de 
spins nucléaires et dont la viscosité ne jouerait pas. 
Le second correspond à la thermalisation d’un neu- 
tron dans un modérateur — encore qu'ici il s'agisse 
d'interaction forte. 

Étant donné que le sort des composantes non 
diagonales de y, nous intéresse, il ne convient pas 
de considérer le second cas, qui correspond à un 
problème classique. Nous supposerons que le sys- 
tème À a un spectre discret. 

Au contraire le thermostat est traité comme un 
système in fini : il a un spectre continu. Nous admet- 
tons que les corrélations, phénomène moléculaire, 
s'organisent suivant le schéma du pas à pas bien 
avant que le système À soit thermalisé : ce phéno- 
mène peut demander plusieurs secondes. C’est un 
phénomène macroscopique. 

Nous utilisons une représentation où 4, et H3 
sont diagonaux. L'opérateur densité ux est dia- 


gonal. On peut poser, T étant la température du 
thermostat : 


Ua — C e—sHa, (1) 


Par contre on ne suppose nullement que u4 soit 
diagonal. Pour le reste nous simplifierons le pro- 
blème au maximum. 

Nous supposerons d’abord que l'interaction a 
une structure simple : elle est « séparable ». V est 
produit de deux facteurs, l’un qui ne concerne 
que À, l’autre qui ne concerne que B : 


Un tE (2) 


Si ere 
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Il en résulte : 
Va = Ua Te Uz us. 


Nous pouvons choisir la décomposition de ma- 
nière que la dernière trace soit nulle : 


DRURUE— 0: 


ce qui a pour conséquence : 


V4 = 0, 

et ce qui fait disparaître le terme gênant (3, 4). En 
même temps les termes en ZL, disparaissent. 

Les jeux permis sur l’ordre des facteurs per- 
mettent ensuite d'écrire les formules suivantes : 

TBlV, Le] us = Ua Tag(Ur L — LU) ua ur, 

U1(Te8 Ug L uB) ua4, 
TB ua uB LV = u4alTs LUB ur) U4, 
U1 AT LU3 uB), 
(TB Uz Luz) ua Uu. 


Il est évidemment commode d’introduire l’opé- 
rateur : 


Haliaur— 


1 


TB Vus uB L 
Ie baUB Ve 


Ja = Tr Up lus, (3) 
qui opère sur À seul. Son adjoint est : 
1: 11} LUIBIUE. (4) 
L’équation d’évolution (2, 20) s’écrit mainte- 
nant : 
À 18 £ 
ua = [Ha —5 iUa Ja + JA Ua) + 
HU TalJa — Ji) ua, val 


1 1 + 
= (Ua Ja + JÂ Us) ta + 43 (Ua Ja + JA Ua) 


2 


— 1 


— Ua 4 JA — Ja pa Ua! (5) 

Pour développer J, il est commode de poser : 

Rob Tr > Ug Upa S(H% H8 - Ha H;) Douce (6) 
LA 


T 


8 Up Hy Je He ee Ha 725 He Hoc 


Cr D, (7) 
œ,B 


Fy et Gy ne dépendent en fait que de la diffé- 
rence 1, — H,. Comme 
Uag Ugo = |Uogl”, 


Fet Gsont réels. Dans la formule (6) permutons a 
et b,œet 6. Il vient : 


Fra = T > Uag Ugo S(Hy + He — Ha — Ho) ueg: 
œ, 


Mais, en vertu des propriétés du facteur de Dirac 
et de la formule (1), on trouve ici : 


URB — est(Æ5—Ha) Ua 
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et par conséquent : 
Fa = eSHÿ—Ha) Fo. (3) 
Nous retiendrons : 
Jab — Uav(Fba + 1Gba), 
Jab = Uab( Far + iGab). 


Nous rendrons le problème tout à fait élémentaire 
en supposant que le système À n’a que deux états, 
a et b. Nous supposerons en outre que U, est 
« antidiagonal » : 


Usa = Ur = 0. 


Nous pouvons nous arranger pour que Us soit 
réel et positif. Enfin nous poserons : 


© — Hp — Ho, 
Fab FT (w), 
Gab Sr G(w), 

I = Usb, 


G — T(J4 — JA)ua. 
Il vient d’abord : 


GT 3) Uab + (Job — Jia) ba 


et : 
d 
baie Uaa —= (J5a | Jia) Uaa 
+ (Job + Jab) Ubb + Gluba — Lab) 
K€ wo (Jap + Jos) Ure 
Le (Jba un Jia) Uaa + G(Uab - Uba)s 
nd ET J ns 
le qi tab — KT © ap — (Ja + Jab) Hab 
+ (Jav + Ja) Ua + Clubs — Ua) 
| | er ne 
Ke Uba = —1iK © ba (Jap + Ja) Hoa 


+ (Jos + Jab) Las + Bla — uv). 


La trace de 4 est une constante, qu'il faut 
d’ailleurs prendre égale à l’unité : 


Uaa + Ubb = Ale 


Nous adopterons de nouvelles variables indé- 
pendantes : 
X = Uab + Ua 
Y = ilblab — Hva), 


I 


ULbb 7 Haa: 
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Les équations d'évolution s’écrivent maintenant : 


KE L — — 9{F(o) + F(—w)] Z—2[F(o) — F(— o)] 
+ 2Y{ (Go) + Gt ©)] X + [F(a) — F(— o)] Y Ii 
(12) 
rene = (13) 
dé 
ET = — |[K—? © + 2G(o) —2G(— ©)] À 


— 9[F(o) + F(— o)] Y 

— 22 {[G(o) + GK o)] X + [F(o) — F(— o)] Y }. 
(14) 

Il est clair que le système est au repos lorsque : 


F(— ©) — F(v) = 10 


L ; ht 
FLO RERO à L 


PENZ) 


Prenons comme variable indépendante Z — Z, 
au lieu de Z : les termes linéaires de l'équation (14) 
s’en trouvent sensiblement modifiés : 


pa dZ = 2) _ 


d Fo) (2 — 20) 


— 2[F(o) 


HE (Go) + Dal ro) Fe or}, 


(15) 
en. TAF 
He NOT, (16) 
Sea OU ee Go) F(—w) — G(—o) Fo) , 
K a LL o + A ae 0) ]x 
F{o) F(— o) 


Flo) + F(-— à) 
— 2(2—Z0) { [G(o) + G(—e)] X + (F0) —F(—w)]} }: 
(17) 


Lorsque le système est proche de la position de 
repos précitée, 1l est permis de négliger les termes 
qui contiennent les produits de Z — Z,, X, Y. Alors 
la composante Z tend sans oscillations vers sa 
valeur de repos avec une constante de temps 


11731 = 2ÆA°[F(o) + F(— w)]. 


Les composantes X et Y ont un mouvement 
oscillant amorti. La fréquence des oscillations 
s’écarte, au second ordre, de sa valeur en l’absence 
d'interaction. La constante de temps qui inter- 
vient pour X et Y est différente de la précédente : 

I Flo) F(-— w) 


To no une 


Ces constantes de temps dépendent de w, évi- 
demment, et aussi de la température. Le rapport 
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entre F(w) et F(— w) est donné par la formule (8). 
Il en résulte la formule suivante : 


T, A1 + es)? 
TNT 0 


Lorsque sw est petit — température élevée, 


_ faible différence d'énergie entre les niveaux a et b 


— les deux périodes sont égales. Au cas contraire 7, 
est plus grand que T,. Il n’y a pas lieu de tenir 
cette conclusion pour générale. L'expression du 
rapport T,/T, subit l’influence déjà notée du terme 
non linéaire de l’équation (2, 14). La valeur limite 
de Z a aussi une expression simple : 


1e 


MP A 


Zo = 


L’état dit «de repos » du système n’est donc pas 
différent de l’état d’équilibre que prévoit la thermo- 
dynamique statistique lorsqu'on suppose l’inter- 
action infinitésimale. 

Ce que nous venons d’examiner, ce sont de petits 
écarts à l’équilibre. Lorsque l’écart est plus consi- 
dérable, il n’est plus permis de négliger les termes 
quadratiques. On ne peut plus prétendre non plus 
intégrer exactement les équations d’évolution. 
Mais il est possible de prouver que le système tend 
vers le même état de repos. Posons : 


_ 4H" Go) Fa conte) 
[0] F(o) + F(— ©) ; 


Cette grandeur est petite devant l’unité, puisque 
nous supposons que la perturbation est faible. On 
obtient aisément la relation suivante : 

1 4 (222723) 


Z 7 \ 2 NN tr 2 2 ; 
5 a LZ Zo) = (4 Ï n) À Dé 1h Te 


Cette relation montre que la quantité : 
(REZ) EUR En) EE 
décroît continuellement jusqu’à ce que l’on ait : 


Z = Z, 
Y = 0) 


Cette remarque ne nous apprend rien sur la 
limite de X. Mais, lorsque Z — Z, et Y sont très 
petits, nous sommes ramenés au problème de 
l'amortissement linéaire qui prouve que X tend 
aussi vers Zéro. 

Finalement l’état de repos a tous les caractères 
que l’on peut exiger de l’état d'équilibre thermo- 
dynamique. 

Les termes quadratiques des équations d’évolu- 
tion peuvent être traités comme des perturbations. 
En première approximation ils feront apparaître 
dans les composantes X et Y des oscillations à la 
fréquence fondamentale, mais plus rapidement 
amortis : 

LIRE 
Ts 1013: 
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cependant que les populations des niveaux, repré- 
sentées par Z ont une composante oscillante du 
double de la fréquence fondamentale, avec une 
constante de temps égale à 2/T.. 


5. Conclusion. — Notre étude a montré com- 
ment l'hypothèse de l’absence de corrélations à un 
instant reculé associée à l’existence d’un spectre 
quasi-continu permet de traiter simplement les 
phénomènes de relaxation qui se produisent dans 
un système légèrement couplé à un thermostat. On 
trouvera peut-être qu'elle apporte, dans les do- 
maines qui ont déjà été explorés et notamment 
dans celui de la résonance magnétique, une certaine 
précision dans l’écriture et dans l’art de formuler 
les moyennes. Elle introduit sans artifice l’état 
d'équilibre du système. Notons encore que notre 
étude fait apparaître un terme quadratique qui est 
indispensable pour des raisons de normalisation et 
dont les effets ne sont pas négligeables, même 
quand les phénomènes sont linéaires. D’autre part, 
il n’est à aucun moment question de diagonaliser 
l’opérateur densité du système étudié. 

L'étude n’a pas fait appel à la théorie tradi- 
tionnelle des transitions provoquées par une per- 
turbation. On sait que cette théorie ne s’applique 
qu’à des intervalles de temps relativement courts : 
assez courts pour que le système ait à peine la 
possibilité de changer d’état. Cette limitation n’est 
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pas due à la nature physique des phénomènes, mais 
au caractère un peu simpliste de cette théorie. 
Notre exposé, qui s'applique même si les popu- 
lations subissent des variations considérables, in- 
corpore somme toute le problème des transitions. 

Si l’on essaie d'imaginer un développement des 
idées présentées ici en laissant de côté les appli- 
cations variées que fait encore prévoir sans autre 
modification le problème du thermostat, il apparaît 
qu'il faut s’orienter vers une description plus fine 
des systèmes étudiés : le thermostat dont nous 
nous sommes servis est par trop impersonnel et 
inerte, il n'est pas visqueux. Il convient soit d’in- 
troduire cette viscosité de manière empirique, soit 
de mettre en jeu tous les degrés de liberté du sys- 
tème complexe étudié ou, aussi bien, toutes les 
particules — ou quasiparticules — dont il est 
constitué. Ce sera alors le moment d’aborder les 
problèmes d’entropie. 

Mais un problème plus simple, déjà signalé au 
début du 42 paragraphe, sera l’interaction de la 
lumière avec des particules chargées. 


Note ajoutée à la correction. — N. N. BocoriuBov, 
Problèmes de dynamique en physique statistique, Moscou, 
1943, traduction du Commissariat à l’Énergie atomique, 
1950, a adopté également l’hypothèse des corrélations 
nulles dans l’état initial. 


Manuscrit reçu le 11 février 1960. 
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LA DÉTERMINATION DE LA DIRECTION D'ARRIVÉE 
ET DE LA POLARISATION DES ATMOSPHÉRIQUES SIFFLEURS 


Première partie : Les dispositifs expérimentaux. 


Par Jean DELLOUE, 


Laboratoire de Physique de l’École Normale Supérieure. 


Résumé. — Les explications avancées pour rendre compte des atmosphériques appelés «siffleurs», 
se basent sur la théorie magnéto-ionique pour prévoir un mode de propagation particulier des 
ondes électromagnétiques de fréquence très basse le long des lignes de force du champ magné- 
tique terrestre et à travers les milieux faiblement ionisés de la très haute atmosphère. La portion 
de l’énergie rayonnée vers le haut par un éclair et survivant aux traversées des couches ionosphé- 
riques est canalisée par les lignes de force et ramenée avec une atténuation assez faible à la surface 
du globe dans l’autre hémisphère. La dispersion importante liée à ce mode de propagation produit 
un effet de sifflement caractéristique (Storey, 1953). 

On a cherché à confirmer cette théorie en déterminant expérimentalement les directions d’arrivée 
des siffleurs qui doivent faire un angle assez faible avec celle du champ magnétique terrestre, et 
leur polarisation qui est liée par une formule magnéto-ionique à cet angle. 

La direction d’arrivée est déterminée par la mesure des retards entre les signaux engendrés par 
un siffleur dans deux paires d’aériens identiques placés aux extrémités de deux bases perpendi- 
culaires. Des relais hertziens sur ondes centimétriques sont utilisés pour retransmettre les signaux 
captés. La polarisation « en projection » est donnée par deux cadres verticaux perpendiculaires 
et la polarisation dans le plan d’onde s’en déduit si l’on connaît la direction d’arrivée. 

Diverses considérations ont conduit à utiliser un appareillage sélectif. Les mesures de retard 
sont alors ramenées à des mesures de déphasages. La fréquence 5 500 c/s a été choisie parce que 
les siffleurs y présentent à peu près toujours une amplitude importante tandis que le fond des 
atmosphériques lointains n’a qu’une amplitude faible. 

Les signaux captés par les aériens sont amplifiés par des amplificateurs basse fréquence de 
gain 80 dB à faible bruit de fond, identiques et stables, puis après retransmission éventuelle par 
relais hertziens, par des amplificateurs sélectifs dont les caractéristiques de phase sont semblables 
et dont la largeur de bande est choisie pour obtenir la meilleure réponse aux variations de fré- 
quence dans le temps le plus fréquemment observées. 

Les tensions amplifiées provenant des deux cadres croisés sont appliquées à un oscillographe 
dit de polarisation. Les tensions engendrées aux deux extrémités de chaque base sont appliquées 
à un circuit qui affiche instantanément le déphasage sur un oscillographe sous la forme de l’incli- 
naison d’une droite ou d’une ellipse très aplatie. Une caméra à déroulement continu lent photo- 
graphie simultanément l’oscillographe de polarisation et les deux oscillographes de déphasages 
auxquels une modulation de brillance particulière est appliquée. 

Les essais sur le terrain ont conduit à modifier à plusieurs reprises l’ensemble de mesures. Des 
circuits correcteurs de phase ont été introduits et on a dû séparer les parties «directions d’arrivée » 
et « polarisation » de l’appareillage. 

Les résultats obtenus et leur interprétation feront l’objet d’un second article. 


Abstract. — The explanations which have been put forward to account for whistling atmo- 
spherics are based upon the magneto-ionic theory and suggest a peculiar mode of propagation of 
very low frequency electromagnetic waves along the lines of force of the earth’s magnetic field 
through the feebly ionized regions of the outer atmosphere. That fraction of the energy radiated 
upwards by a lightning stroke which survives the crossing of the ionospheric layers, travels along 
the lines of force and returns to the earth surface in the opposite hemisphere. The high dis- 
persion attached to that mode of propagation accounts for the whistling effect (Storey, 1953). 

Às a check to the foregoing theory, we tried to determine the directions of arrival of whistlers, 
which should make a rather small angle with the direction of the magnetic field, and their polari- 
zation which is related through a magneto-ionic formula to that angle. 

The directions of arrival are determined by measuring the time lag between the signals induced 
by whistlers in two pairs of identical aerials located at the ends of two perpendicular bases. Micro- 
wave radio links are used to re-transmit the induced signals. The ‘* projected ” polarization is 
obtained from two vertical crossed loops, and the plane wave polarization can be deduced provided 
the direction of arrival is known. 

Various considerations led to the choice of a selective apparatus. Time lag measurements are 
then ordinary phase lag measurements. A frequency of 5 500 cps has been selected because the 
amplitude of whistlers is for most of the time large in that part of the spectrum, while the level of 
long distance atmospherics is low, 
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DIRECTION ET POLARISATION DES ATMOSPHÈRES SIFFLEURS 


The signals induced in the aerials are fed to high gain, wide band, low noise amplifiers with 
stable and identical characteristics, and after being re-transmitted over radio links. to selective 
amplifiers with identical phase characteristics and a convenient bandwidth to obtain the best 
response to a signal varying in tone according to the law most often observed. The voltages 
corresponding to the crossed loops are then applied to a ‘* polarization ? oscilloscope. The signal 
corresponding to the ends of each base are fed to a phase comparator which instantaneously displays 
the phase lag as the slope of a straight or slightly elliptic pattern on an oscilloscope. A conti- 
nuous motion slow speed camera photographs simultaneously the polarization oscilloscope and 
the two phase oscilloscopes to which convenient blanking voltages are applied. 

Fields tests of the equipment led to successive alterations of the set-up. Phase correcting 
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networks were introduced and the ‘‘ polarization *” section of the equipment had to be separated 


from the ‘‘ directions of arrival ? section. 


The results obtained and their interpretation will be discussed in another article. 


Introduction. — Le travail présenté dans le 
pages qui suivent est consacré à l’étude de cer- 
taines propriétés d’un type de perturbation électro- 
magnétique d’origine naturelle du domaine des très 
basses fréquences. 

Afin de mieux situer le problème, il nous a paru 
nécessaire de résumer rapidement dans une pre- 
mière partie les traits caractéristiques des diverses 
perturbations électromagnétiques d’origine natu- 
relle et plus particulièrement de celles affectant la 
bande des très basses fréquences, puis de préciser 
l’état actuel des connaissances sur le type de per- 
turbation qui nous intéresse : « l’atmosphérique 
siffleur ». Nous sommes dès lors en mesure d’ex- 
poser l’objet de ce travail : vérifier si les siffleurs 
se propagent bien le long des lignes de forces du 
champ magnétique terrestre en déterminant leur 
direction d’arrivée et si la polarisation observée 
peut s’expliquer par un tel parcours à travers les 
couches supérieures de l’atmosphère. 

La seconde partie est consacrée à l’exposé des 
méthodes de mesures utilisées et à la discussion 
de certains aspects du problème et de leur influence 
sur la conception de l’appareillage : nécessité de 
favoriser le plus possible les siffleurs vis-à-vis des 
autres parasites, problèmes posés par le caractère 
fugitif des siffleurs et leur occurrence impré- 
visible, etc. 

La troisième et la quatrième partie décrivent les 
différents éléments de l’appareillage et les modifi- 
cations successives de l’ensemble de mesures telles 
qu’elles se sont révélées nécessaires au cours des 
expériences. 

L'interprétation, assez délicate, des diagrammes 
observés sur les oscillographes de mesure, tant pour 
la détermination des directions d’arrivée que pour 
la polarisation, est discutée ensuite. Une méthode 
simple de représentation des caractéristiques d’un 
siffleur est proposée. 

Enfin, une dernière partie résume les résultats 
obtenus qui confirment le mode de propagation 
le long des lignes de force du champ magnétique 
terrestre et une tentative y est faite de remonter 
de la polarisation observée aux paramètres affec- 
tant la propagation. ; | 


1. Les parasites électromagnétiques 
d’origine naturelle. 


1-1. — LES PERTURBATIONS ÉLECTROMAGNÉ- 
TIQUES BASSE FRÉQUENCE. — Les décharges élec- 
triques orageuses sont les sources naturelles les 
plus importantes de perturbations dans tout le 
spectre des ondes hertziennes. Elles se manifestent 
sous la forme de claquements brefs (quelques di- 
zaines de millisecondes) ou de crachements lorsque 
l’éclair origine a été complexe. Le niveau moyen 
de ces parasites, assez faible dans la gamme des 
ondes décimétriques et métriques, s’élève très 
rapidement lorsqu'on atteint les ondes décamé- 
triques et de radiodiffusion, mais c’est vers les fré- 
quences basses et très basses qu’on observe les 
niveaux les plus élevés, en particulier aux environs 
de 10 kc/s. Le phénomène est encore très bref et 
l’impression auditive est toujours celle d’un claque- 
ment mais le spectre présente en général dans cette 
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en fonction de la distance. 


gamme de fréquences un aspect remarquable dé- 
pendant de la distance du récepteur à l’éclair 
source (fig. 1) : pour des distances de 2 000 km ou 
plus on constate une absence de composantes dans 
la bande 2 000 à 6 000 c/s. Nous utiliserons d’ail- 
leurs cette particularité ultérieurement, 


516 


Mais, en plus des atmosphériques classiques ap- 
paraissent dans la gamme des très basses fréquences 
des parasites d’un genre tout à fait différent, de 
plus longue durée et à tonalité musicale. 

Le premier type appelé «tweek » est un atmo- 
sphérique musical relativement court dont la durée 
est de l’ordre de 20 millisecondes. Sa forme d’onde 
est caractérisée par une oscillation amortie qui 
suit une impulsion brève. La fréquence de l’oscilla- 
tion décroit d’abord très rapidement puis reste 
sensiblement constante à la limite, aux environs 
de 4 000 à 2 000 c/s. Les tweeks sont produits par 
les arrivées successives des divers échos d’un atmo- 
sphérique formés par réflexions successives sur la 
terre et la basse ionosphère et s’observent prin- 
cipalement au lever et au coucher du soleil. 

Le second type appelé «siffleur » consiste sous 
sa forme la plus courante en un sifflement dont la 
hauteur décroît depuis des fréquences inaudibles, 
rapidement tout d’abord, puis plus lentement pour 
les fréquences basses jusque vers 1 000 c/s. Le 
parasite balaye ainsi plusieurs octaves en l’espace 
d’une ou deux secondes. C’est à ce genre de parasite 
que se rapporte notre étude, aussi ferons-nous plus 
loin une étude approfondie de celles de ses carac- 
téristiques qui sont actuellement bien connues. 

Les autres types de parasites basse fréquence 
sont d’occurrence plus rare et liés à une forte acti- 
vité magnétique : souffle continu affectant la 
bande 1 000 à 4000 c}s pendant des heures, 
« chœur de l’aube » formé de sifflements de hauteurs 
croissantes sur un fond continu de bruits semblables 
à des gazouillis, ou siffleurs ascendants à évolution 
rapide. 


1.-2. — LES PARASITES SIFFLEURS. — Bien que 
leur existence soit connue depuis longtemps, leur 
étude systématique n’a été entreprise que récem- 
ment. Cependant, dès 1918 Barkhausen [1] en- 
tendit « voler les grenades » dans des appareils de 
télégraphie par le sol et conclut très vite que la 
cause de ces parasites était extérieure à l’appareil- 
lage. Par la suite, les mêmes phénomènes se mani- 
festèrent dans les lignes téléphoniques à longue 
distance et dans tout aérien de dimensions assez 
grandes auquel on connectait un amplificateur 
basse fréquence. Eckersley [2] s’intéressa à cet 
aspect de la question et proposa le premier en 1925 
une explication basée sur la dispersion d’une im- 
pulsion électromagnétique dans une atmosphère 
d'ions. Barkhausen [3]émit en 1930 une théorie qui, 
si elle convenait particulièrement au tweek, ne 
pouvait rendre compte des caractéristiques des 
siffleurs. C’est néanmoins lui qui, le premier, vit 
dans la décharge orageuse une origine possible des 
siffleurs. Vers 1933 Burton et Boardman [4], étu- 
diant des câbles sous-marins, publièrent de nom- 
breux résultats sur les tweeks ainsi que le premier 
oscillogramme et la première analyse spectrale d’un 
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siffleur. En 1935, Eckersley [5], reprenant son an- 
cienne idée relative à la dispersion d’une impulsion 
électromagnétique dans l’ionosphère, et supposant 
que la correction de Lorentz, due à la polarisation 
du milieu, n’est pas à faire, mit en lumière une 
possibilité de propagation dans le mode extraor- 
dinaire des ondes électromagnétiques de très basses 
fréquences à travers l’ionosphère et trouva pour la 
variation de la fréquence d’un siffleur dans le 
temps une loi en 1/#?. Cette loi était confirmée 
par l’analyse spectrale de Burton et Boardman. 

Ce n’est qu’en 1953 que Storey [6] reprend la 
question du point de vue expérimental et théo- 
rique et complète la théorie d’'Eckersley par une 
hypothèse très intéressante. Nous allons examiner 
ses résultats en détail. 


1.-2. 1. — Résultats expérimentaux. — Storey 
montre qu'on peut classer les siffleurs forts en 
deux catégories suivant qu’ils sont ou non précédés 
d’un fort claquement atmosphérique. Au moyen 
d’un réseau de goniomètres cathodiques, il découvre 
que ces derniers sont dus à des éclairs et 
proviennent de zones orageuses distantes de 
moins de 2 000 km du récepteur. 

Analysant au moyen d’un spectrographe basse 
fréquence les dispersions de nombreux siffleurs 
enregistrés sur bande magnétique, il vérifie la rela- 
tion proposée par Eckersley : 


== Din? 


D étant appelée dispersion du siffleur. De plus, sur 
un graphique représentant fl? en fonction de #, 
la droite représentative d’un siffleur coupe toujours 
axe des temps à l’instant où s’est produit le fort 
claquement cause présumée du siffleur. Enfin, la 
mesure des dispersions de nombreux siffleurs 
montre que ceux qui sont précédés d’un fort cla- 
quement ont toujours une dispersion approxima- 
tivement double de celle des siffleurs ne jouissant 
pas de cette propriété. 

Lorsque les siffleurs arrivent en groupe, ce n’est 
pas nécessairement le seul fait du hasard. Storey 
distingue les « trains de siffleurs » formés d’un sif- 
fleur et de plusieurs échos (jusqu’à une dizaine se 
succédant à intervalles réguliers avec des ampli- 
tudes de plus en plus faibles et des dispersions 
croissantes. Sur les diagrammes 7 —12(#) les droites 
représentatives concourent en un même point de 
l’axe des temps et, s’il y a eu un claquement atmo- 
sphérique avant le train, il s’est généralement 
produit à cet instant et les dispersions sont 
comme 1-2-3-4 tandis que dans le cas contraire 
elles sont comme 1-3-5-7, Les «siffleurs d’éclair 
multiple » se suivent irrégulièrement à de courts 
intervalles, sont d’amplitudes comparables et ont 
la même dispersion : leur origine est vraisembla- 
blement dans un éclair multiple. Enfin les « pai- 
res » ressemblent à un train de deux siffleurs mais 


Or: 


leurs amplitudes sont du même ordre et leurs dis- 
persions ne sont pas dans un rapport simple. 

Une étude des variations systématiques des pro- 
priétés des siffleurs par des moyennes portant sur 
des périodes de l’ordre du mois montre, outre 
l'irrégularité de leuroccurrence, qu’ils se produisent 
bien plus fréquemment de nuit que de jour et 
que les siffleurs longs se font entendre surtout l’été 
et les courts, l’hiver. 


1.-2. 2. — La théorie de Storey-Eckersley. — La 
loi t — D.f1? résultant de la théorie de Eckersley 
se trouve vérifiée par l’expérience. Il y a donc tout 
lieu de penser que les siffleurs sont bien produits 
par la dispersion d’atmosphériques impulsifs le 
long d’un trajet ionisé. Mais il reste à déterminer 
ce trajet, et en particulier à expliquer comment 
les siffleurs longs, ayant leur source dans une dé- 
charge orageuse locale, reviennent dans la région 
où ils ont pris naissance avec le caractère que nous 
leur connaissons, et pourquoi ils sont d’autant plus 
forts que l’éclair-source s’est produit plus près. 
D'autre part, l’origine des siffleurs courts est à 
déterminer. 

La théorie magnéto-ionique décrit la propagation 
d’une onde plane monochromatique dans une di- 
rection donnée, à travers un milieu ionisé et en 
présence d’un champ magnétique, au moyen de 
l’indice de réfraction y correspondant à une vitesse 
de propagation uw. Dans le cas des siffleurs, il s’agit 
en fait de trains d’ondes qui, lorsque le milieu est 
dispersif, voyagent avec une vitesse de groupe v 
qui est différente de w. Si le milieu est amisotrope, 
le train d’ondes se déplace dans une direction fai- 
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sant un angle « non nul avec la normale à la surface 
d’onde et tel que (fig. 2): 


Dans cette formule 0 est l’angle entre le champ 
magnétique et la normale et & est positif si la nor- 
male est située entre le champ magnétique et la 
direction de propagation. Dans un milieu à la fois 
anisotrope et dispersif, le train d’ondes voyage 
dans la direction de propagation avec la 
vitesse V = p/cos «. 

La formule donnée pour p? par la théorie de 
Appleton-Hartree peut se simplifier dans le cas où 
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les collisions sont négligeables et où l’on admet 
l’approximation quasi-longitudinale de Booker. 
Cette dernière est valable aux fréquences basses 
pour des directions de propagation faisant des 
angles assez grands avec le champ magnétique et 
pour des densités électroniques très variables. Il 
vient alors, en omettant le terme de polarisation 
de Lorentz et compte tenu des valeurs des para- 
mètres à attendre dans le cas qui nous intéresse : 


u? = æfyr avec x — pè[p? — Ne?]mf? 
et y = (pxfp) cos 0 — (Hefmc) cos 0, 
PQ étant la pulsation critique du milieu et pg la 
pulsation correspondant à la gyrofréquence. 


L’indice correspondant à la vitesse de groupe est, 
dans les mêmes conditions : 


= (1/2) (x{yc)"? = (1/2) w 
et l’angle « est défini par : 
tg &« — — (1/2) tg 0. 


Le signe moins montre que la direction de pro- 
pagation est toujours située entre le champ magné- 
tique et la normale à l’onde plane. 

Si, au moyen de ces formules, on construit la 
surface donnant la vitesse de propagation en fonc- 
tion de la direction de propagation, on voit que, 
même pour des valeurs de x/y petites, Les directions 
de propagation possibles sont limitées à un cône 
étroit axé sur le champ inagnétique. Le train 
d’ondes est donc guidé par le champ magnétique 
dans le milieu ionisé (fig. 3). 


X/Y = 5 
HTC 


X/y = 00 


X/Y = 1 


Les indices de réfraction pour des fréquences 
basses sont très grands dans l’ionosphère (de 10 
à 300 pour le mode extraordinaire et parallèle- 
ment au champ magnétique, plus dès que l’on 
s’en écarte). Les ondes émises par un éclair et 
pénétrant dans l’ionosphère suivant une gamme 
d’angles d'incidence très étendue voient leurs di- 
rections de propagation rapprochées énergiquement 
du champ magnétique par l'effet de réfraction. 
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Puis, dans la haute atmosphère, l'effet magnéto- 
ionique continue à focaliser l’énergie dans la di- 
rection du champ magnétique. On peut par ailleurs 
montrer que le trajet suivi est pratiquement in- 
dépendant de la fréquence. 

Storey se basant sur les remarques précédentes 
suppose que les siffleurs sont des ondes produites 
par des décharges orageuses qui ont traversé 
l’ionosphère dans le mode extraordinaire et suivi 
les lignes de force du champ terrestre. A l’ar- 
rivée dans l’autre hémisphère, il reste encore une 
quantité d’énergie appréciable grâce à la foca- 
hisation. La dispersion produit l’effet de sifflement 
caractéristique. 

Un siffleur court capté dans l'hémisphère Nord 
a pour origine un éclair de l’hémisphère Sud dont 
l’atmosphérique ordinaire, propagé entre la terre 
et l’ionosphère, est trop atténué par la longueur 
du trajet pour être détectable. Un siffleur long 
capté dans l’hémisphère Nord a pour origine un 
éclair proche de l'hémisphère Nord (qui produit 
un atmosphérique audible) les ondes ayant par- 
couru l’aller-retour le long de la ligne de force avec 
réflexion sur le sol dans Jl’hémisphère Sud. La 
focalisation joue pour rendre le siffleur long audible 
dans un rayon relativement faible autour de l'éclair 
source (fig. 4). Les trains de siffleurs s’expliquent 


Siffleur court 
reçu ici Nn 
Eclair 


Siffleur. long 


reçu ici 


Éclair 


F1G.4. — Trajectoires d’unsiffleur long et d’un siffleur court. 


par des aller-retour successifs. La dispersion étant 
proportionnelle au trajet parcouru, les trains consé- 
cutifs à un claquement atmosphérique ont des 
dispersions 2D, 4D, 6D, etc.., ceux non précédés 
par un claquement ont des dispersions D, 3D, 
5D, etc. L’occurrence de siffleurs plus nombreux la 
nuit que le jour est lié à l'absorption dans la basse 
ionosphère et l’asymétrie journalière à l’activité 
orageuse toujours plus forte dans l’après-midi. 
Enfin, l’abondance de siffleurs longs l’été est liée 
à la variation saisonnière des orages et celle des 
siffleurs courts l’hiver à l’alternance des saisons 
entre les deux hémisphères. 


N° 6 


1. 2. 3. — Les conséquences de la théorie de Sto- 
rey et l’objet de ce travail. — Une première consé- 
quence importante de cette théorie est de montrer 
qu'il ne faut pas inclure le terme de Lorentz dans 
la formule de la polarisation. La validité de ce 
terme est contestée depuis longtemps et aucune des 


“expériences imagirées pour trancher la question 


ne s’est révélée concluante. 

Une seconde conséquence est l’existence d’un 
grand nombre d’électrons (400 à 1 000 par cm) 
à des distances de l’ordre du rayon terrestre de la 
surface du globe. De telles densités sont nécessaires 
pour expliquer la propagation des siffleurs et les 
fortes dispersions observées. Ce phénomène semble 
d’ailleurs avoir été confirmé par des mesures sur 
la lumière zodiacale. Cette ionisation serait due 
à un apport de matières ionisées en provenance 
du soleil et expliquerait certaines occurrences si- 
multanées de siffleurs et d’orages magnétiques ou 
d’aurores. 

On imagine aisément l’intérêt présenté par les 
siffleurs pour l’investigation de la très haute atmo- 
sphère jusqu’à plusieurs rayons terrestres. La 
structure du champ magnétique terrestre, la gran- 
deur et la variation de l’ionisation à ces hauteurs 
se reflètent dans la dispersion, la forme d’onde 
et la polarisation des siffleurs. Nous nous sommes 
donc proposés tout d’abord de vérifier la théorie 
de Storey-Eckersley en déterminant l’angle d’ar- 
rivée des siffleurs par rapport au champ magnétique 
terrestre et d’entamer l’étude entrevue de la très 
haute atmosphère terrestre en enregistrant simul- 
tanément la polarisation de ces siffleurs afin de 
tenter de la relier à l’angle d’arrivée, On sait en 
effet que la théorie magnéto-ionique permet d’éva- 
luer en même temps que l’indice de réfraction pu, 
Ja quantité À — — h,[h, caractérisant le champ 
électromagnétique de l’onde plane. 


2. Les dispositifs expérimentaux. 


2.-1. — LES MÉTHODES DE MESURE. — La détermi- 
nation des directions d’arrivée est obtenue par deux 
mesures de déphasages. Les siffleurs induisent des 
signaux dans trois aériens disposés aux trois sommets 
d’un triangle rectangle. La station où s’effectue 
les mesures est située au sommet de ce triangle et 
nous l’appellerons station centrale, les deux autres 
étant les stations auxiliaires. Nous appellerons base 
l’ensemble constitué par la station centrale et une 
station auxiliaire. Nous mesurons les déphasages 
entre Je signal de référence capté à la station 
centrale et les signaux des stations auxiliaires. Des 
formules simples relient ces déphasages aux angles 
(azimuth, élévation) définissant la direction d’ar- 
rivée du siffleur, Les fréquences en jeu étant très 
basses, on doit, pour mesurer un déphasage 
notable, utiliser des bases relativement longues 
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(quelques dizaines de kilomètres). Le problème 
majeur est alors de retransmettre les signaux 
captés par les stations auxiliaires à la station 
centrale sans altérer leurs caractéristiques. Ceci 
est effectué au moyen de deux ensembles de relais 
hertziens fonctionnant sur ondes centimétriques. 

D’autre part, les aériens de mesures doivent être 
placés dans des endroits suffisamment dégagés, et 
éloignés d’accidents du sol trop importants afin 
de capter une onde aussi peu déformée que possible. 
Ces conditions, jointes à celles im posées par l’éta- 
blissement des liaisons hertziennes rendent très 
difficile la découverte d’emplacements convenables 
à des distances relativement grandes. 

La mesure de la polarisation se fait à la station 
centrale au moyen de deux cadres perpendiculaires 
verticaux. On obtient ainsi en fait la polarisation 
«en projection », comme nous l’appellerons, mais 
la connaissance de la direction d’arrivée permet, 
en principe au moins, de passer de la polari- 
sation en projection à la polarisation dans le plan 
d’onde. 


2.-2. — LE CARACTÈRE SÉLECTIF DES MESURES. 
— La bande de fréquences où se manifestent les 
siffleurs est assez grande, de 4 500 c/s à 15 000 c/s 
et même plus. Il ne semble donc pas a priori 
s'imposer d’opérer à fréquence fixe. Mais, deux 
caractères des siffleurs nous ont amenés à les ob- 
server dans une bande relativement étroite de 
fréquences. Tout d’abord la caractéristique sif- 
flante du siffleur indique que des portions différentes 
du spectre basse fréquence sont affectées successi- 
vement par le siffleur, tandis que les atmosphé- 
riques classiques se manifestent simultanément 
dans toute la bande de fréquences. En opérant en 
bande large nous capterions donc à un instant 
déterminé toute l’énergie rayonnée par les atmo- 
sphériques à cet instant, mais une portion seule- 
ment de l’énergie totale du siffleur. Or, ces atmo- 
sphériques perturbent les mesures de phases et de 
polarisation : nous aurions donc un rapport signal/ 
bruit (naturel) peu satisfaisant. Cet effet est d’au- 
tant plus gênant que les siffleurs sont, la plupart 
du temps, beaucoup plus faibles que les atmo- 
sphériques présents (principalement ceux d’origine 
proche) : c’est leur variation de tonalité qui permet 
à l'oreille de les distinguer encore nettement lorsque 
les appareils de mesure ne donnent que des indica- 
tions infimes. D’autre part, le désir que l’on a 
d'éviter toute déformation du front d’onde in- 
terdit d’utiliser le relief avoisinant comme pro- 
tection par effet d’écran. 

C’est donc l’idée de protection à l’égard des 
atmosphériques qui nous a conduits dès le début 
à opérer dans une bande de fréquences relative- 
ment étroite. Or, il est possible de choisir cette 
bande de façon à minimiser encore l’action des 
atmosphériques et à renforcer celle des siffleurs. 


En effet, la propagation dans le guide d’ondes un 
peu particulier constitué par la surface de la terre 
— pratiquement parfaitement réfléchissante — et 
les couches inférieures de la basse ionosphère — où 
l’on à plutôt réfraction progressive que réflexion 
métallique — provoque dans le spectre des atmo- 
sphériques ordinaires un affaiblissement d'autant 
plus marqué des fréquences entre 2 000 et 6 000 c/s 
que la distance de propagation est plus grande. 
D'autre part, observation courante montre que 
la plus grande partie de l’énergie des siffleurs est 
localisée dans la partie basse du spectre basse 
fréquence, de 1 500 à 6 000 c/s. C’est donc cette 
région qui est la meilleure. Nous opérons à 5 500 c/s 
pour des raisons techniques : les relais hertziens 
employés ayant une bande passante utile allant 
de 4 000 à 28 000 c/s, nous avons évité de nous 
placer trop près de la limite inférieure 4 000 c/s 
pour éviter les rotations de phase assez importantes 
dues aux filtres passe-haut utilisés. 

L'exploitation de l’appareillage nous a apporté 
une justification supplémentaire de ce point de 
vue : le spectre des siffleurs n’est pas constitué 
par une fréquence pure décroissant dans le temps 
mais par une bande de fréquences plus ou moins 
complexe dont la fréquence centrale dérive dans 
le temps. En bande large, les motifs observés sur 
les oscillographes de mesure seraient inexploitables 
parce que résultant de la superposition de plusieurs 
motifs élémentaires en général três différents dûs 
aux diverses composantes spectrales présentes à 
un instant donné. 


2.-3. — LE CHOIX DE LA LARGEUR DE BANDE. — 
Vis-à-vis des circuits sélectifs utilisés, le siffleur 
se comporte comme un signal de tonalité variant 
des aiguës aux graves pour lequel on peut définir 
une fréquence instantanée : c’est la fréquence de 
la composante de Fourier dont la phase est sta- 
tionnaire par rapport à la fréquence à l’instant 
considéré. Dans ces conditions, on peut montrer 
que la réponse d’un circuit sélectif à une telle tona- 
lité est la meilleure et le rapport signal/bruit est 
optimum lorsque la largeur de bande B est prise 
sensiblement égale à la racine carrée de la varia- 
tion de la fréquence dans le temps. L’emploi d’une 
largeur de bande plus faible donnerait à la sortie 
du circuit sélectif un tension dont l’amplitude de 
crête serait nettement inférieure à celle qu’on ob- 
tiendrait pour une fréquence fixe égale à la fré- 
quence propre du circuit sélectif. De plus, dans 
ce cas, l'enveloppe de la tension de sortie présente 
des oscillations parasites dont l’amplitude est d’au- 
tant plus forte par rapport à la réponse initiale 
que la largeur de bande est plus faible. Il est pré- 
férable d'éviter ces phénomènes. 

A la latitude à laquelle nous avons fait nos ob- 
servations, la dispersion moyenne des sifileurs 
courts était de l’ordre de 50 s—1”?, celle des siffleurs 
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longs de 100 s—V2, Or, si la dispersion est repré- 
sentée par la loi 


t—\D:f 3 on no — — SP? et par suite 


di 
df court < 
[] — 16 000 c/s/s et B opt. — 125 cycles. 
dt]5.500 c/s 
df long 
Ë] — 8000 c/s/s et Bopt. — 90 cycles. 
dt]5.600 c/s 


En fait, la largeur de bande utilisée est, pour des 
raisons techniques, beaucoup plus grande, de 
l’ordre de 500 cycles. Il est difficile autrement 
d'obtenir des réponses en amplitude et phase 
suffisamment semblables pour tous les amplifi- 
cateurs sélectifs utilisés dans l’appareillage. 


2.-k. — RÉALISATION DE LA SÉLECTION DE FRÉ- 
QUENCE. — Il faut, en particulier dans les mesures 
de déphasages, que les caractéristiques de phase 
et d'amplitude des divers canaux de l’appareillage 
soient identiques. Cette condition est plus délicate 
à remplir dans une station mobile et la division de 
deux des canaux en deux tronçons par les relais 
hertziens rendrait difficile le contrôle global de la 
sélectivité de l’appareillage. Nous avons donc pré- 
féré travailler en bande large à l’entrée de chaque 
canal, les parties sélectives de l'équipement étant 
toutes groupées à la station centrale où 1l est 
facile et rapide de vérifier leur bon fonctionnement 
en cours de mesures. Un avantage de cette méthode 
est de pouvoir éventuellement changer de fréquence 
d’étude sans avoir à modifier l’équipement des 
stations mobiles. 

Enfin, l’emploi d’une bande large permet l’écoute 
et l’identification des siffleurs sans équipement 
annexe au poste central et aux postes mobiles 
ainsi qu'une vérification rapide à la sortie des 
relais hertziens de la bonne marche des canaux 
auxiliaires. 


2.-5. — LE MODE DE PRÉSENTATION DES RÉSUL- 
TATS. — Comme nous l’expliquerons plus bas, 
les informations fournies par l’appareillage sont 
disponibles pendant un temps relativement court, 
de l’ordre au maximum de 1/5 seconde, résultant 
de la largeur de bande de l’appareillage, et de la 
« largeur naturelle » des siffleurs. Nous avons donc 
eu recours à des oscillographes cathodiques pour 
afficher ces informations et à un enregistrement 
photographique pour en permettre l’étude a pos- 
Lertor. 

La présentation de la polarisation en projection 
s’obtient simplement sur un tube cathodique dit 
«de polarisation » en formant une figure de Lissa- 
jous avec les signaux provenant des deux cadres 
croisés de la station centrale. 

Mais, pour la présentation de chacun des dé- 
phasages relatifs aux deux bases, nous avons choisi 
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de traduire directement le déphasage par un angle. 
Un circuit décrit plus loin fait apparaître sur l'écran | 
du cathodique « de déphasage » un segment ou une 
ellipse assez fortement aplatie. L’angle de ce seg- 
ment, ou du grand axe de cette ellipse, avec la 
direction de déflection correspondant à l’une des 
paires de plaques du tube fournit la mesure du demi- | 
déphasage entre les tensions induites aux deux 
extrémités de la base correspondante. Cette mé- 
thode est infiniment plus commode et plus précise : 
que celle qui consisterait à mesurer le déphasage à 
partir d’une ellipse de Lissajous ordinaire. 

2.-6.— LE PROBLÈME DE L'ENREGISTREMENT DES 
CARACTÉRISTIQUES DES SIFFLEURS. — L’occurrence 
des siffleurs est totalement imprévisible. D’autre 
part, la réalisation d’un appareil différenciant les 
siffleurs des autres types d’atmosphériques en un 
temps suffisamment court pour que ce renseigne- 
ment soit utilisable et permette de photographier 
les caractéristiques de ces siffleurs est très délicate 
et le fonctionnement peu sûr. Il est donc nécessaire 
de filmer en permanence les oscillographes : on 
utilise une caméra 16 mm modifiée pour permettre 
un déroulement continu du film à la vitesse de 
JC. 

D'autre part, 1l faut utiliser une modulation dis- 
continue de la brillance des cathodiques afin de 
distinguer les motifs successifs traduisant l’évo- 
lution du siffleur. Les impulsions de déblocage du 
wehnelt doivent cependant être assez longues 
(1/100 sec.) afin que le spot puisse décrire une 
tranche de siffleur d’une cinquantaine de cycles 
sur laquelle des atmosphériques d’amplitude équi- 
valente ne produiront que des perturbations par- 
telles ne détruisant pas la cohérence de l’ensemble 
du motif. La fréquence de récurrence de ces im- 
pulsions à été choisie de l’ordre de 30 par seconde 
pour donner une séparation suffisante entre deux 
tranches consécutives. Enfin on a trouvé commode 
de ne pas éteindre complètement le spot entre deux 
impulsions de brillance : on peut ainsi observer 
l’enveloppe des siffleurs ce qui facilite leur localisa- 
tion sur le film. 


3. Les différents éléments 
de l’ensemble de mesures. 


3.-1. — RELAIS HERTZIENS. — La partie de 
l’appareillage qui nous était imposée était 
les relais hertziens : ce sont des équipements 
M X 641 C P fabriqués par la C. S. F., et conçus 
principalement pour l’établissement de tronçons 
hertziens dans des quartes téléphoniques à courants 
porteurs. Ils fonctionnent sur ondes centimé- 
triques, entre 1 700 et 2 000 mégacyeles, avec une 
puissance à l’émission de 300 milliwatts. La por- 
teuse est modulée par des impulsions « en densité » 
dont la fréquence de récurrence moyenne est 
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115 kc/s. Une voie de service utilise la partie 
200-3 000 c/s du spectre basse fréquence transmis. 
Le reste de la bande, de 4 000 c/s à 28 000 c/s est 
réservé aux « courants porteurs ». Le rapport 
signal/bruit est meilleur que 55 db, et des ampli- 
ficateurs basse fréquence permettent d'utiliser à 
Pémission des niveaux allant de 2,5 à 0,025 volts. 

Ce matériel avait plusieurs caractéristiques in- 
téressantes pournos mesures. [l est peu encombrant, 
facile à déplacer et à installer. C’est une considé- 


ration très importante pour les stations mobiles 


qui n'existent pas en permanence, mais doivent 
être rétablies à chaque séance de mesures. De plus, 
sa consommation faible le rend particulièrement 
bien adapté à une alimentation par accumulateurs 
et commutatrice. La fréquence utilisée et le carac- 
tère directif de l’émission évitent toute induction 
parasite dans les cadres de mesures. Ce phénomène 
qui s’est manifesté, par exemple, avec des équi- 
pements émetteur-récepteurs sur ondes métriques 
conduit à une saturation dès Ja première lampe de 
l’amplificateur basse fréquence avec démodula- 
tion et amplification énergique de cette modulation 
par le reste de l’appareillage. La portée maximum 
de ces émetteurs est d’environ 60 km en vue 
directe, c’est-à-dire de l’ordre de grandeur d’une 
longueur d’onde pour les fréquences étudiées. Ils 
sont donc largement suffisants à ce point de vue. 
Enfin, la voie de service permet la liaison entre les 
stations mobiles et la station centrale, et l’exécu- 
tion des opérations de contrôle de l’appareillage 
en cours de fonctionnement. La bande « courants 
porteurs » permet la retransmission en «bande 
large » au-dessus de 4 000 c/s des signaux basse 
fréquence captés, les parties sélectives de l’appa- 
reillage se trouvant toutes à la station centrale. 
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3.-2. — AÉRIENS. — L’étude de la polarisation 
en projection nécessite l'emploi de deux cadres 
croisés. Dans la réalisation initiale chacun de ces 
aériens fournissait en plus la tension de référence 
pour la mesure du déphasage de la tension induite 
dans l’aérien de la station mobile correspondante. 
Les aériens des stations mobiles étaient donc par 
nécessité des cadres. 

Dans un but de standardisation du matériel, 
on a donné les mêmes dimensions et les mêmes 
caractéristiques à tous ces cadres. Les dimensions 
externes sont faibles afin de permettre le transport 
de tout le matériel d’une station mobile dans une 
voiture de tourisme. Les cadres sont carrés, de 
1 mètre de côté, 10 centimètres de large, et com- 
portent deux couches de 20 spires chacune de fil 
15/10 isolé sous plastique. Les deux couches sont 
connectées en série, le point milieu étant relié à 
la masse du cable coaxial bifilaire assurant la 
liaison avec l’amplificateur : ce montage symé- 
trique permet d'éviter au maximum les effets 
d'antenne, Sur le terrain, les cadres sont main- 
tenus verticaux par des haubans sur des supports 
démontables qui les surélèvent d’environ 1 mètre 
au-dessus du sol. 

La hauteur effective de ces cadres à la fréquence 
5 500 c/s est d’environ 0, 5 cm. 


3.-3. — AMPLIFICATEURS € BANDE LARGE ». — 
Les aériens sont reliés à des amplificateurs basse 
fréquence qui doivent avoir un gain important, 
de l’ordre de 80 db, une réponse uniforme dans la 
bande de fréquences étudiée mais aussi une atté- 
nuation importante aux fréquences où des émet- 
teurs ondes longues sont susceptibles de se mani 
fester, et enfin être stable vis-à-vis des variation 
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Fic. 5. — Amplificateur « bande large ». fe 
: CEA RER ? xl: 71e Dre £ if a céle 
C : cadre ; F : filtre ; a, a’ : sortie vers la partie émission d’un relais hertzien ou vers un amplificateur sélectif, 
39 


922 


éventuelles des tensions d’alimentation (surtout 
pour les stations mobiles). 

L'amplificateur utilisé (fig. 5) comporte quatre 
lampes groupées en deux paires séparées par un 
filtre. Les deux premières lampes sont des triodes 
12 À X 7, afin de minimiser le bruit de fond. La 


grille de la première triode est attaquée par le” 


secondaire d’un transformateur d’excellente qua- 
lité et de grand rapport de transformation (120), 
dont le primaire, à basse impédance, est relié à 
l’aérien correspondant. Une contre-réaction assez 
énergique stabilise l’ensemble des deux triodes. 

Le filtre laisse passer la bande de fréquences 
2 500-11 000 c/s. La coupure supérieure évite les 
interférences des postes ondes longues relative- 
ment proches et de fréquences voisines : G. B. R. 
(Rugby) sur 16 000 c/s et F. U. B. (Houilles) sur 
17 800 c/s. On a préféré supprimer les fréquences 
inférieures à 2 500 c/s en raison des nombreuses 
interférences industrielles rencontrées dans cette 
gamme. 

Les deux dernières lampes sont des pentodes 
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E F 40 montées avec contre-réaction et attaquées 
par un diviseur de tension donnant les atténua- 
tions 1/1 1/2, 1/5, 1/10. La sortie s'effectue à 
travers un transformateur d’adaptation d’im- 
pédances. 

Le gain maximum en tension de l’amplificateur 
est de 120 000. Le bruit de fond ramené à l’entrée 
est inférieur à 0,2 microvolts. La saturation dans 
la position 1/1 est obtenue pour 5 microvolts à 
l’entrée. Tous les amplificateurs de ce type ont 
été réalisés aussi semblables que possible, tant 
vis-à-vis du gain aux diverses fréquences que de la 
courbe de rotation de phase en fonction de la fré- 
quence. Cette identité est préservée même pour les 
variations des tensions d'alimentation à attendre 
dans les stations mobiles. 


3.-4. — AMPLIFICATEURS SÉLECTIFS. — Un dé- 
phaseur résistance-capacité d’impédance moyenne 
est attaqué par l’entrée basse impédance à travers 
un inverseur introduisant à volonté un déphasage 
de 1800 (fig. 6). La sortie du déphaseur se fat à 
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Fic. 6. — Amplificateur sélectif. 


D : déphaseur ; I : 


inverseur ; b, b’: 
d’un amplificateur bande large ; ©, ©’: 


haute impédance, une partie de la charge étant 
constituée par un potentiomètre de réglage de 
gain. Ces conditions d’impédance donnent au dé- 
phaseur une réponse sensiblement constante en 
amplitude. 

L’étage sélectif est constitué par une pentode 
E F 40 chargée dans la plaque par un circuit ac- 
cordé et dont la polarisation est ajustable afin 
de régler au mieux le point de fonctionnement de 
la lampe. Les variations minimes des caractéris- 
tiques d’un tube à l’autre sont en effet suffisantes 
pour introduire des anomalies dans la réponse en 
phase et en amplitude. Il est même nécessaire de 
sélectionner ces tubes. L’étage de sortie est formé 
d’une E F 40 montée en cathode follower à très 
haute impédance d’entrée afin d’éviter de charger 
le circuit sélectif. 

Le gain maximum en tension de l’ensemble est 


entrée reliée à la partie réception d’un relais hertzien ou à la sortie 
sortie vers un calcu- lateur de phase ou vers l’oscillographe de polarisation. 


de 100. La plage de variation de la fréquence 
d'accord est de 1 000 cycles autour de 5 000 cys. 
La rotation de phase maximum utilisable donnée 
par le déphaseur est de 1000. Le réglage de l’en- 
semble des amplificateurs sélectifs pour l’obtention 
de courbes de réponse identiques est de loin le plus 
délicat de l’appareillage. Heureusement, il est rela- 
tivement stable. 


3.-5. — CALCULATEURS DE PHASE. — Étant 
données les deux tensions 


V;, = à, cos ot Ve = a2 cos (ot + o) 


on forme les tensions somme et différence et on 
avance (ou diminue) de +/2 la phase de la tension 
somme par exemple. On effectue ensuite un Lissa- 
Jous entre la tension ainsi obtenue et la tension 
différence, Un calcul simple montre qu’on observe 


. Ne 


alors sur l’oscillographe une ellipse dont le grand 
axe fait un angle 6/2 avec une paire de plaques, 
et dont le rapport des axes est lié homographique- 


ment au rapport des amplitudes des deux tensions 
initiales : 
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S1 les deux tensions initiales ont même ampli- 
tude (ai — 4) on observe une droite faisant avec 
cette paire de plaques un angle w/2. Si les deux 
tensions sont presque de même amplitude, l’ellip- 
ticité observée est un moyen très sensible de véri- 
c We 5 É Nu + " nl . 
fier l’égalité des tensions. Enfin, si l’une des ten- 
sions est nulle, Pellipse observée devient un cercle : 


entrées reliées aux amplificateurs sélectifs ; 


+ . nn . . . 
Eh Eh; Ev 64: Sorties vers les amplificateurs horizontal et vertical de l’oscillographe de phase. 


l’aspect de l’oscillographe fournit donc un moyen 
instantané de contrôle des réglages. 

Les tensions somme et différence sont obtenues 
au moyen d’un montage où les transformateurs 
sur pots Ferroxcube ont été réalisés avec un très 
grand soin pour avoir des rapports de transforma- 
tion égaux (fig. 7). Ces transformateurs sont atta- 
qués par deux cathodes follower à 12 À X 7. Les 
tensions somme et différence sont appliquées, après 
passage à travers deux autres cathodes follower, 
Pune à une paire de plaques de l’oscillographe de 
déphasage, l’autre à une 6 À U 6 montée en in- 
tégrateur Miller. La tension prélevée dans la plaque 
est en quadrature avec celle injectée dans la grille, 
L’intégrateur est réglé à 5 500 c/s mais son fonc- 
tionnement est correct dans une bande de 
500 cycles. Les erreurs introduites par ce calcula- 
teur dans les mesures de phase sont inférieures 
aux erreurs de lecture sur l’oscillographe. 


3.-6. — LES OSCILLOGRAPHES DE PHASE ET DE 
POLARISATION. — Les tubes cathodiques utilisés 
sont du modèle 5 S P II À à deux canons et à 
fluorescence bleue. L’un des tubes est utilisé pour 
les deux mesures de phase, l’autre, dont un seul 
canon sert, pour la polarisation. Un petit tube 
néon placé entre les deux cathodiques permet de 
localiser les siffleurs sur le film : l’opérateur l’allume 
quelques instants lorsqu'il perçoit l’un de ces pa- 
rasites. 
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F16. 8. — Amplificateur d’oscillographe 

et circuit de modulation de brillance. 
ff” : entrée reliée à un amplificateur sélectif ou à un calcu- 
lateur de phase ; xx’ : sortie vers les plaques de l’oscillo- 
graphe ; C : cadrage ; G : gain ; E : équilibrage ; www” 
vers les wehnelts des oscillographes ; D : durée, N : niveau. 
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Les sensibilités des deux paires de plaques de ces 
tubes cathodiques étant voisines, on utilise Je 
même modèle d’amplificateur sur les deux voies 
(fig. 8) : l'attaque des plaques se fait par deux tubes 
E F 80 montés en amplificateur continu différentiel. 
Un étage préamplificateur continu à E C C 81 


comporte les organes de réglage de gain et d’équi- 


librage. La bande passante s’étend du continu 
à 150 kc/s et la rotation de phase avec la fréquence 
est négligeable jusqu’à 70 kc/s. La sensibilité est 
de l’ordre de 1 volt crête à crête pour la déflection 
totale. 

Enfin un dispositif de modulation de brillance 
permet de moduler les whenelts des trois canons 
avec des ondes carrées de durée et de période 
réglables fournies par un multivibrateur à 6 J 6. 


4, Les versions successives de l’ensemble 
de mesures et leur mise en œuvre. 


4-1. — LA PREMIÈRE VERSION. — Dans le but de 
réduire au maximumlenombre d’élémentsnécessai- 
res, on a d’abord adopté un montage où les aériens 
de polarisation en projection servaient également, 
en association avecles cadres respectifs des stations 
mobiles, aux mesures de déphasages. La disposi- 
tion spatiale des cadres vus à la verticale était donc 
celle représentée ci-contre (fig. 9a). Le diagramme 
de l’appareillage était alors assez simple : quatre 
amplificateurs large bande seulement étaient uti- 
lisés. Les portions des quatre canaux comprises 
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F1G. 9. — Version I. 


LÉGENDE DES DIAGRAMMES DE L’'APPAREILLAGE. 


Fic. 9, 11, 12. 
A : amplificateur bande large. 
RH : relais hertzien. 
(@ : correcteur de phase. 
D : déphaseur. 
S : amplificateur sélectif. 


CP : calculateur de phase. 

OD : oscillographe de déphasage. 

OP : oscillographe de polarisation. 
YYY : contacteurs de calibration. 

m,mMm, : commutateurs de mode. 

81823 : réglages de gain des oscillographes, 


bp To: 


réglages d'amplitude. 


entre les aériens et les sorties des amplificateurs 
sélectifs étaient identiques, à part l’introduction 
dans deux d’entre eux des deux relais hertziens dont 


le gain global en tension était choisi égal à l'unité | 


(fig. 9) 


Un défaut très gênant qui s’est révélé à l’usage 


provenait de la rotation de phase introduite par 
les relais hertziens. La complexité des circuits de 
ces derniers provoque une variation de déphasage 
de 1000 par 1 000 cycles au voisinage de 5 500 c/s. 
Étant donnée la largeur de bande (500 c/s) utilisée, 
on observe une rotation du diagramme de dépha- 
sage très importante lorsqu'on applique un signal 
d'amplitude constante et de fréquence variable à 
l’entrée des deux canaux correspondants. 

D'autre part, le rapport signal sur bruit de l’ap- 
pareil s’est révélé insuffisant, la majorité des sif- 
fleurs observés à l’époque se distinguant à peine 
du bruit de fond.  : 


4-2 — DEUXIÈME VERSION. — La sensibilité 
de l’équipement a été facilement augmentée en 
utilisant des transformateurs de couplage aériens- 
cadres de rapport de transformation beaucoup plus 
élevé. On ne pouvait envisager d'augmenter les 
dimensions des cadres en raison des servitudes 
d'espace aux stations mobiles. 

Pour compenser la rotation de phase introduite 
par les relais hertziens, on a inséré aux points xx 
du précédent diagramme (fig. 9), deux correcteurs 
dont la caractéristique de phase était aussi proche 
que possible de celle des relais hertziens. Les con- 
tingences techniques et le désir d’éviter une trans- 
formation complète de la réalisation matérielle de 
l’équipement ont conduit à utiliser des correcteurs 
purement passifs. Ces circuits, tout en provoquant 
un déphasage qui varie rapidement avec la fré- 
quence, doivent avoir de plus une réponse en am- 
plitude constante dans toute la bande passante 
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F1G. 10. — Circuit correcteur de phase. 
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utilisée : ceci a pu être obtenu en couplant deux 
cellules en T comme indiqué sur la figure 10. Le 
problème est rendu très complexe par les varia- 
tions importantes de l’impédance caractéristique 
dans la bande passante. Le couplage par condensa- 
teur des deux cellules influe aussi beaucoup. 
D’autre part, de légères retouches sur les divers 
éléments de chaque cellule permettent de modifier 
légèrement l’allure de la courbe de phase du cor- 
recteur et de la rendre aussi semblable que pos- 
sible à celle des relais hertziens. 

Le contrôle de l'efficacité de ces correcteurs 
s'effectue en rassemblant tout le matériel en un 
même point, les stations mobiles auprès de la 
station centrale, et les relais hertziens fonctionnant 
sur des distances de quelques mètres. La seule 
différence avec le fonctionnement réel réside alors 
dans le déphasage dû au temps de propagation 
entre l’émetteur et le récepteur de chaque relais 
hertzien. Les quatre aériens étant débranchés, on 
injecte à leur place au moyen d’un générateur basse 
fréquence les quelques microvolts nécessaires et 
on fait varier la fréquence. On doit alors observer 
sur les deux écrans de déphasage, deux segments 
qui varient d'amplitude mais ne tournent pas. Les 
correcteurs de phase permettent d’obtenir ce ré- 
sultat à mieux que 1° près dans la gamme 
5 500 + 200 c/s, et à mieux que 2° près dans la 
gamme 5 500 - 500 c/s. Un petit effet résiduel 
de distorsion du segment représentatif subsiste 
néanmoins du côté des fréquences basses, mais 
pour des amplitudes vraiment très fables. Vu 
l’aspect complexe des oscillogrammes de siffleurs, 
on peut considérer que le défaut est suffisamment 
corrigé. 

Cependant, deux nouveaux défauts sont apparus 
à l’usage. Tout d’abord l’augmentation de sensi- 
bilité a permis de constater que l’emploi de deux 
paires d’aériens de directions perpendiculaires pour 
la mesure des déphasages n’était pas satisfaisant. 
En effet, si les tensions induites dans les deux 
cadres croisés de polarisation ne sont pas du même 
ordre de grandeur, le diagramme de déphasage 
correspondant à la plus faible des tensions est 
difficilement lisible. Autrement dit, l’état de pola- 
risation des siffleurs peut rendre l’une ou l’autre 
des mesures de déphasages impossible. L’observa- 
tion ayant montré à l’époque que les cadres d’orien- 
tation Est-Ouest étaient les plus souvent défa- 
vorisés, nous avons décidé de faire les deux me- 
sures de déphasages à partir de trois aériens 
orientés Nord-Sud. D'autre part, chacun des deux 
correcteurs étant réglé pour compenser au mieux 
les défauts du relais hertzien correspondant, les 
deux canaux de la polarisation en projection 
n'avaient plus des caractéristiques de phase iden- 
tiques. Même en se résignant à adopter pour les 
deux correcteurs une formule moyenne qui 6om- 
pense au mieux les deux canaux hertziens, il n’était 
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jamais possible de les rendre suflisamment iden- 
tiques pour que les mesures de polarisation n’en 
souffrent pas. 


4-8. — TROISIÈME VERSION. — La nouvelle 
disposition spatiale des cadres devint alors celle 
représentée sur la figure 114. Le diagramme d’ap- 
pareillage en était sensiblement modifié et perdait 
son aspect symétrique. Un amplificateur sélectif 
supplémentaire était relié à l’amplificateur large 
bande Nord-Sud de la station centrale par l’inter- 
médiaire d’un seul correcteur de phase et consti- 
tuait ainsi le canal de référence pour les deux 
mesures de déphasage (fig. 11). 
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Mais cette nouvelle version s’est révélée rapide- 
ment peu satisfaisante car les contigences de réa- 
lisation ne permettaient pas de coupler l’ampli- 
ficateur large bande Nord-Sud au correcteur de 
phase par une lampe d'isolement. Ce dernier ayant 
une impédance d’entrée extrêmement variable en 


Frc. 12. — Version IV. 
À gauche partie « déphasages », à droite partie € pola- 
risation », 
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fonction de la fréquence, réagissait sur le canal 
de polarisation sur lequel il était branché et rendait 
les mesures de polarisation inutilisables. 


4-4, — QUATRIÈME VERSION. — Nous avons 
alors résolu de séparer complètement la partie 
«polarisation » de la partie « déphasages » abou- 
tissant ainsi au diagramme de la figure 12. Nous y 
voyons un cadre, deux amplificateurs large bande 
et un amplificateur sélectif de plus que dans la 
première version. Mais cette complication du.ma- 
tériel n’affecte pas les stations mobiles et l’indé: 
pendance des mesures, et donc celle des réglages, 
est assurée, Cet ensemble fonctionne parfaitement. 


4.-5. — LES OPÉRATIONS DE CONTROLE. — 
4-5. 1. — Contrôle partiel (réalisable même en cours 
de mesures). — Il est possible de contrôler rapide- 


ment le fonctionnement de la partie de l’appareil- 
lage située au delà des entrées des amplificateurs 
sélectifs : des contacteurs indépendants, dits de 
calibration, permettent d'isoler complètement 
chaque voie, ou bien de la relier à un générateur 
basse fréquence. D’äutre part, deux commuta- 
teurs dits de mode m, et m,, permettent d'observer 
sur les oscillographes de phase, soit l’une ou l’autre 
des tensions de sortie des amplificateurs sélectifs, 
soit un diagramme de Lissajous entre ces deux 
tensions, soit le fonctionnement normal à travers 
ie calculateur de phase. Enfin les oscillographes 
sont munis de réglages de gain g, £o £a. 

On effectue sur la partie « déphasages » les opé- 
rations suivantes : 

a) On règle successivement les amplificateurs 
sélectifs en mettant les contacteurs de calibration 
sur la position « générateur basse fréquence » et 
en observant les tensions de sortie au moyen des 
contacteurs de mode, puis on les retouche en les 
observant deux par deux en Lissajous pour obtenir 
un segment d'amplitude variable mais de direction 
fixe quand la fréquence varie. 

b) On met le contacteur de calibration de la voie 
commune sur la position « générateur basse fré- 
quence », les deux autres commutateurs étant sur 
«isolé » et les commutateurs de mode sur «calcula- 
teur ». On doit alors observer sur chaque oscille- 
graphe un cercle et, au besoin, on amène l’ellipse 
observée à devenir un cercle en agissant sur les 
réglages de gain de l’oscillographe correspon- 
dant (g, ou g»). 

c) On met les trois contacteurs de calibration 
sur la position « générateur basse fréquence » et 


N°°6 
les deux contacteurs de mode sur « Lissajous » : 
on observe deux ellipses qu’on transforme en 
droites au moyen des déphaseurs à l’entrée des 
amplificateurs sélectifs. On amène ensuite ces 
droites à faire un angle de 45° environ avec les 
directions des paires de plaques des oscillographes 
au moyen des potentiomètres de réglage d’am- 
plitude pp', qq’. 

d) On met enfin les commutateurs de mode sur 
la position « calculateur » : il apparaît une droite 
ou une ellipse très allongée qu’on transforme en 
droite en retouchant légèrement p ou q. Ces droites 
peuvent être horizontales ou verticales suivant la 
position des inverseurs de déphaseurs. On choisit 
l'orientation la meilleure pour la lisibilité des dia- 
grammes de siffleurs comme expliqué plus loin. 

Par construction, si l’on agit sur les déphaseurs, 
toute augmentation du retard de phase d’un canal 
à relais hertzien par rapport au canal de référence 
de la station centrale fait tourner la droite corres- 
pondante dans le sens des aiguilles d’une montre. 

Le réglage de la partie « polarisation » est beau- 
coup plus simple : seule l'opération a) est à effec- 
tuer. On règle ensuite jes gains de l’oscillographe 
de polarisation pour que le segment observé fasse 
un angle de 450 avec les directions des paires de 
plaques de l’oscillographe de polarisation. 


4-5. 2. — Contrôle total. — Cette opération 
nécessite, comme expliqué au paragraphe « deu- 
xième version » de faire fonctionner tout l’appareil- 
lage au même endroit. Elle permet de contrôler 
la marche correcte des tronçons de l’équipement 
inaccessibles en cours de mesures. En fait, le seul 
réglage est celui de la tension de sortie des relais 
hertziens, ceux-ci ayant un gain différent de l’unité. 
On effectue l’opération d) précédente, en transfor- 
mant l’ellipse allongée par retouche du potentio- 
mètre de gain. Si l’on fait varier la fréquence du 
générateur basse fréquence, on observe, en général, 
en plus de la variation d'amplitude, une très légère 
distorsion du segment, tandis que sa direction ne 
change pas. On retouche, si nécessaire, les dépha- 
seurs pour obtenir des segments de déphasage nul 
exactement verticaux ou horizontaux. Ce contrôle 
total est une opération très longue exigeant la mise 
en batteries de tout l’appareillage. Elle n’est pas 
effectuée fréquemment et l’expérience montre 
qu’elle est heureusement rarement utile car la sta- 
bilité des réglages est très grande. 


Manuscrit reçu le 20 novembre 1959, 
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Résumé. — Dans le but de vérifier si la règle des sommes pour les triplets (152p 3P9,,,,-1s4s 3,) 

et (152p 8P9,,,6-154d D le l de l’héli intensités de 
SPP EU )3,2,1) de l'atome de l’hélium est valable, les rapports des intensités des 

composantes de ces triplets ont été mesurés, d’une part en tenant compte de la réabsorption, 
d'autre part, dans des conditions telles que la réabsorption n’intervienne pas. 

Ainsi, le rapport mesuré des intensités de la somme des composantes proches des raies 4 713 À 
et 4 472 À respectivement à celle de la composante extrême est : 

7,5 + 0,1 au lieu de 8, obtenu théoriquement 

9,7 + 0,1 


et 


Abstract. — In order to verify the sum rule for the helium atom triplets : (152p 3P9,,,0-184s 9.) 
and (1s2p %P9,,,9-1s4d °D3,2,1), the ratios of these triplet component intensities have been mea- 
sured, partly by taking into account reabsorption, partly in such conditions that the reabsorp- 
tion may be neglected. Thus the ratio measured of the sum of the intensities of the two near 
components of the lines 4 713 and #4 472 respectively to that of the remote component are : 

je 0,1 


TOME 6, JUIN 1960, PAGE 527, 


ATEN O 


instead of 8, obtained theoretically. 


Introduction. — Il est connu, depuis longtemps, 
que la règle des sommes d’Ornstein, Burger et 
Dorgelo n’est pas valable pour les doublets de la 
série principale des spectres des métaux alcalins 
(S. Sambursky{1]) et qu’elle n’est pas satisfaite non 
plus dans les doublets du spectre de l’aluminium, 
du gallium, de l’indium et du thallium (0. U. Von- 
viller [2]; S. E. Williams et J. Herlichy [3]; 
R. Payne-Scott [4]; D. Kunisz [5]). Les écarts à 
cette règle augmentent avec le numéro du doublet 
dans la série. D'autre part, ils sont plus grands 
pour les doublets correspondants des éléments plus 
lourds. E. Fermi a expliqué ce phénomène par 
l’interaction spin-moment orbital (/-s). 

Il est connu, d’autre part, que la structure très 
serrée des niveaux triplets de l’hélium est anormale 
et qu’elle ne peut pas être interprétée tout sim- 
plement par des modèles de couplage dans 
l’atome L-S ni 7-j. Ainsi, on a observé de grands 
écarts à la règle des intervalles de Landé surtout 
pour les niveaux plus bas. 

Il est intéressant d'examiner si Ja règle des 
sommes est satisfaite pour les triplets de l’hélium 
et surtout pour les triplets provenant de la combi- 
naison des termes inférieurs. 

Les expériences ont été faites pour les triplets 
suivants : 

Série IT secondaire 

(SD DEP ES PSAS EST) AE 718: 

Série [| secondaire 

(1s2p 3P9,,,0-154d 3D;,2,1) À 4 472. 


Les rapports des intensités pour ces triplets ont 
déjà été en partie mesurés par D. Burger [7]; 
W. V. Houston [9]; J. Brochard, R. Chabbal, 
M. Chantrel et P. Jacquinot [10]. Pourtant on n’a 
pas encore discuté l’influence des conditions d’exci- 
tation du spectre et en particulier la réabsorption, 
sur les rapports des intensités des composantes des 
triplets. Aussi il n’était pas possible de dire si, 
d’après les résultats précédents, la règle des sommes 
est valable ou non. 

Dans le présent travail, on à examiné comment 
les rapports des intensités dans les triplets men- 
tionnés dépendent des conditions d’excitation du 
spectre pour trois différentes sources lumineuses. 


I. Dispositif expérimental. — SOURCE LUMI- 
NEUSE. — On a utilisé trois sources lumineuses : 
1. Cathode creuse refroidie par l’azote liquide, ali- 
mentée par un générateur de tension stabilisée 2 kV 
Elle est représentée sur la figure 1. 2. Cathode 
creuse refroidie par eau, dont le schéma est donné 
sur la figure 2 ; elle est alimentée par le même 
générateur. 3. Tube Geissler refroidi par de l’azote 
liquide, alimenté de la même façon. Il est repré- 
senté sur la figure 3. 


SPECTROGRAPHE, INTERFÉROMÈTRE ET DISPO- 
SITIF D'ENREGISTREMENT. — On a utilisé un spec- 
trographe de la S. G. O. avec un interféromètre 
Fabry-Pérot à balayage par variation linéaire de 
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pression. Les lames de l’interféromètre, fournies 
par Hilger, ont un diamètre de 10 em ; elles sont 
recouvertes de 5 couches diélectriques et ont une 
finesse efficace supérieure à 20. 


Fenêtre d'observation 


Pièce en laiton 
reliée à l'anode 


Cylindre 
en allechart 


Anode. 
| en aluminium 


Cathode 
en laiton 


Fenètre d'observation 


Anode. 
en aluminium 


Vers piège 
et robine 


Cathode creuse 
en laiton refroid'e 
par eau 


TG 


Le dispositif expérimental se compose d’une cel- 
lule photoélectrique Lallemand à 20 étages, d’un 
adaptateur d’impédance et d’un millivoltmètre 
enregistreur Philips, 


Cathode G , È 
en aluminium Fenêtre. 
refroidie d'observation 
RETIEeS Anode 


en aluminium 


TUE 


II. Résultats. — Pour les deux triplets, on a 
mesuré le rapport de la somme des intensités des 


FEG 4. 


raies appelées 7, et Z, sur la figure 4, à l'intensité 
de la raie appelée Z,, soit : 


TS AS2 PE PASSE ST) PRE NT 
+ I (152p SP9-As4s 85,) à 4 713,1 
T,(152p *PI-1sks5S;) 104719,4 


Fu désignerons dans le texte ce rapport par 
A 
(==) et le rapport : 
1 4713 


13(152p °PS-AS4d D, À 4474,5 
+ L(152p SPIA1SL4S D) RER 
1,(152p $P$-1s4d SD.) x 4 471,7 


Cake 

li 4472 

Pour avoir des résultats sûrs, il vaut mieux ne 
pas résoudre les deux composantes fortes des tri- 
plets. On à excité le spectre dans les conditions 
suivantes : 

1) Dans la cathode creuse refroidie par azote 


liquide (fig. 1) remplie d’hélium sous la pression 
de 1,1 et 4,0 mm de mercure. 


P AT 
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2) Dans la cathode creuse refroidie par eau 
(fig. 2) remplie d’hélium (pression : { mm Hg) et 
de néon (pression : 3 mm Hg). 

3) Dans le tube Geissler (fig. 3) rempli d’hélium 
(pression : 1,4 mm Hg) refroidi par azote liquide ou 
à température ambiante. 

On a mesuré les rapports d’intensités 


(£ + I a (£ : 1) 
li 4713 li 4472 


en fonction de la pression du gaz dans la source 
lumineuse et de l’intensité du courant électrique 
dans la source. 


Peer 
a) VARIATION DES RAPPORTS | —5—? un 
li 4713 


T3 + 
——— EN FONCTION DE LA PRESSION DU 
L: 4472 
GAZ DANS LA SOURCE. — On a constaté que : 


1. Pour les petites valeurs de l’intensité du cou- 
rant électrique dans la source, les rapports 


(= —) "1 (= —") 
nn 4713 L 4472 


diminuaient quand la pression augmentait. 

2. Pour les valeurs plus grandes de l’intensité du 
courant dans la source, les rapports précédents 
augmentaient avec la pression du gaz. Les courbes 
expérimentales sont données sur la figure 5. 


hi 
1 


4718 


Pression mm Hg 


les 


Les deux triplets ont les mêmes niveaux infé- 
rieurs 3P9 ,,,. Le niveau supérieur du triplet de là 
série II secondaire est simple (singlet) tandis que 
celui de la série I secondaire est triple mais très 
peu résolu. 

Si la réabsorption n'existait pas, le rapport 


FF +) Bafavë + RofoVa 


serait égal à : ë 
l; gli vi 


OÙ £1; 
4713 : Re : 
4 sont les poids statistiques des niveaux 


ne 
152p SPO1s52p °P 152p 3P£; fus fo Ja t Vus Vas Vs 
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les forces d’oscillateurs et les fréquences des raies 
(1s2p SP$-As4s 3S,), (1s2p 3P9-As4s3S.), 
(1s2p 3P2-1s4s 35). 
Ensuite, sous l’effet de la réabsorption, l’inten- 
sité de chaque raie diminue et devient : 
DANSE 
Te Se la 
La Sala 


où l'est l'épaisseur de la couche absorbante et les 8% 
k —1, 2, 3) dépendent uniquement du pro- 


: 2e? Nz fr TU 
duit Al. Le AE SAT 
(S. E. Frisch[11]) donne l’absorption au centre de la 
raie. V, est la population du niveau inférieur de la 
raie. Si À = 0, $3 = 1. Quand le produit X, L aug- 
mente, S diminue et devient nul lorsque le pro- 
duit X, l'augmente indéfiniment. 

Ainsi en tenant compte de la réabsorption, on 
obtient : 


(nn en 
nn 4713 li | li 


facteur …r K 


ÉÊRE RE _ 82/2 V3 92 
Bafa Si gif VW S1 

Le facteur k est proportionnel à V4 fx /vr. 

Le rapport des intensités des raies dépend alors 
de la population de leurs niveaux inférieurs. 

On peut supposer que pour les petits courants 
dans la source, les rapports de la population des 
niveaux 1s2p SPS(N.,), 1s2p SP9(N;), 1s2p SP(N;) 
sont égaux aux rapports des poids statistiques de 
ces NIVEAUX : 


Ni Ni: Nas = gi: 82:88 


D'autre part, pour toutes les composantes d’un 
seul triplet de la série IT secondaire, f; — constante, 
ce qui donne : 


INR 
kiikoiko= ti DR gigi gs —=1:3:5. 
RE PSN 
Comme 
k: < ko <, Rs, S; > S 9 > 15%, 


(= E =) (:) S3 © (3 S3 
Î; 4713 T; 4713 5; ll 4713 91 
et comme 


BE Ho 
Ds RES EN 


GAL AT 


IIS SE +) < (= de #) 
JL 4713 JT, 4713 


On peut voir de la même façon que : 


is) <(5) 
li 4472 li 4472 


c’est-à-dire que la réabsorption diminue les rap- 
ports des intensités des raies mentionnées. Ceci 
n’est vrai naturellement que si les populations des 
niveaux inférieurs sont du même ordre de grandeur, 
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Comme la réabsorption augmente avec la pres- 
sion du gaz dans la source, les rapports : 


( 3e 2 “ É “e 2) 
li 4713 li 4472 


doivent diminuer. 

Ceci est observé pour les faibles courants dans la 
source, Par contre, pour les courants plus impor- 
tants, les rapports 


() ” (= "A 
in 4713 f Î; 4472 


augmentaient avec la pression du gaz dans la 
source. 

Pour expliquer ce phénomène, il faut tenir 
compte du cas où les populations des niveaux 
1s2p 3P$, 1s2p 3P? et 1s2p $P$ diffèrent beaucoup. 

Pour les valeurs élevées du courant électrique 
dans la source, une transition interdite du niveau 
1s2p 3P£ au niveau 1s? 18, a lieu, ce qui diminue 
la population A; 

Les collisions thermiques des atomes du gaz ne 
peuvent pas compenser ces différences des popu- 
lations V,, NV,, N,, car la durée de vie de l’atome 
sur le niveau 1s2p *P? est d’un ordre de grandeur 
plus petit que l’intervalle de temps entre deux 
collisions dans nos conditions expérimentales. 

C’est pourquoi, pour une assez grande valeur du 
courant dans la source, il peut se faire que N, <N,; 
RETENIR ER 


Il peut donc se faire que : 


(2 & 1 + (£ + 1 
Î; 4713 un 4713 


Dm di 
l; 4472 li 4472 


c’est-à-dire que pour les grands courants, l’absorp- 
tion des raies Z, et Z, augmente avec la pression du 
gaz dans la source, tandis que la raie Z, n’est prati- 
quement pas absorbée ; c’est pour cela que les 


rapports 
Het) 
l' ans li ae 


augmentent avec la pression du gaz dans la source. 


et 


Ia +1 
Î; Jan 


| AVEC L’INTENSITÉ DU COURANT 
4472 


b) VARIATION DES RAPPORTS | ET 


(£ + 
li 
DANS LA SOURCE LUMINEUSE. 
La figure 6 représente les courbes expérimentales 


T3 +1 
. à PA . 
des rapports (sr en fonction du courant 
1 
électrique dans la source, 


N°6 


Pour la raie (1s2p 3P0-1s4s %S), cette courbe aug- 
mente avec le courant plus rapidement pour des 
pressions du gaz plus élevées dans la source. Par 
contre, pour la raie (1s2p 3P9-1s4d 3D) le rapport 


9 Z, "t 
( Des diminue quand le courant augmente 
li 4472 
la+la 
pa 2) TG 1,4 mmHe 800°K 4718 (8) 
cc 1,1-1,2 mmHe 4718 (1) 
8 47 
7 F TETE) OR He Ne ie 
H N 2 + e PE 
ce He Net) pe 4472 (2) 


__— 


\ 
NTG 14 mmHe 4472 (3) 


70 m 


F16c. 6. 
C. G. : Cathode creuse. 
T. G. : Tube Geissler (3). re 
(1) : Cathode creuse refroidie par l’azote liquide. 
(2): Cathode creuse refroïdie par l’eau. 


Ceci peut s’expliquer ainsi : quand le courant 
électrique augmente dans la source lumineuse, 

I. Le nombre d’atomes excités croît, ainsi que 
les populations des niveaux 1s2p BP 6 Ve IVe 
D’après la discussion précédente, c’est pour cela 

I3 +1 Re 
que les rapports (= — diminuent. 
1 
IT. L'intensité de la raié interdite 


(SAS 162 np) 


augmente, ce qui cause la diminution de la popu- 
lation NV, du niveau 1s2p %P?. La réabsorption de 
la raie (1s2p P9-1s4s 3S.) est environ 5 fois moins 
probable que celle des raies (1s2p 3P9-1s4d 3D,.,). 
C’est le raison pour laquelle, lorsque le courant aug- 
mente, la raie (1s2p Pî-1s4s 38.) peut être très 


peu absorbée et le rapport (ste) croit avec 
qi 4713 
le courant tandis que les raies (1s2p 3P9-1s4d 3D,.,) 


T3 + IT, 

Des CU 41472 
diminue quand le courant croit dans la source lumi- 
neuse, 

On peut observer que l'influence de la réabsorp- 
tion sur les rapports des intensités des composantes 
des triplets de l’hélium est semblable à celle décrite 
par Z. Les et H. Niewodniczanski [11] pour les 
triplets du cadmium, 


sont encore absorbées et le rapport 


MEN 6 


II. Conclusions. — Dans le but de vérifier si la 
règle des sommes est valable pour les triplets 
(152p ?P0-1s4s 35) et (1s2p 3P0-1s4d D) de l’hélium, 
il faut mesurer les rapports des intensités des com- 
posantes de ces triplets dans des conditions où la 
réabsorption n’intervienne pas. Ceci est possible 
lorsqu'on remplit la source lumineuse avec un 
mélange d’hélium et d’un autre gaz rare (par 
exemple du néon), ce qui diminue énormément la 
probabilité d’absorption. Dans ces conditions on 
obtient les rapports : 


I + +) 
———— 7 01 
V4 4472 ; 
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231 


qui ne dépendent ni de lintensité du courant 
électrique dans la source, ni de la pression du 
gaz. 

Ces valeurs montrent de petites déviations par 
rapport à la règle des sommes, qui donne pour ces 
triplets : 


SÉRIE 
( à IE 


Ces déviations sont probablement causées par 
lPinteraction spin-moment orbital. 


LAERUTS Manuscrit reçu le 12 janvier 1960. 
(=) = DIN) TI 0,1 
l; 4713 
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REVUE DES LIVRES 


Joueur (M.), Traité d’Électricité théorique, tome III. 
Principe de relativité et lois générales de l’Électromagné- 
tisme et de l’Électrodynamique. (Un vol. br. de 179 p., 
16 X 24 em, Gauthiers-Villars, Paris, 1960.) 

Ouvrage d’une très grande clarté et d’un très grand 
intérêt. L’exposé est fait systématiquement du point de 
vue de la Relativité restreinte, ce qui rend simple et facile 
des questions habituellement très ardues, telles que les for- 
mules de Liénard, sur les champs d’une charge ponctuelle 
en mouvement. Les développements mathématiques sont 
concis et parfaitement corrects ; tout l'ouvrage est d’une 
lecture agréable et aisée, et il s'adresse au public scienti- 
fique dans son ensemble. Voici l’ordre des matières : 
Chap. L.: Transformation de Lorentz-Poincaré et espace- 
temps de Minkowski. A la fin, dynamique relativiste et 
applications (Charge ponctuelle dans les champs uniformes) ; 
Chap. II : Le Champ électromagnétique et les équations de 
Maxwell (potentiels retardés, charges ponctuelles en mou- 
vement, leur interaction, etc.) ; Chap. III. : Le Champ 
électromagnétique dans la matière et électrodynamique 
des corps en mouvement ; Chap. IV : États quasi-perma- 
nents. Généralités sur les forces électromotrices d’irduc- 
tion et sur les courants induits ; Chap. V : Ondes électro- 
magnétiques. Énergie électromagnétique et rayonnement. 
L'ouvrage se termine sur les notions de masse électro- 
magnétique et de moment cinétique électromagnétique. 
L'examen de l’électromagnétisme classique, introduction 
nécessaire aux théories quantiques, se trouve donc complet. 


J. WINTER. 


Zaman (J. M.), Electrons and Phonons (théories des phéno- 
mènes de transport dans les solides). (Un vol. relié de 
554 p., 16 X 24 em, Clarendon Press : Oxford Univer- 
sity Press, 1960, prix : 84 shillings.) 

L'objet est l'étude générale des phénomènes électroniques 
et vibratoires dans les solides ; métaux, semi-conducteurs 
et isolants. La division de l'ouvrage est très claire, et cela 
le rend accessible malgré l’immense quantité de phéro- 
mènes couverts. 

Ce livre, malgré son étendue, en raison du grand dévelop- 
pement des théories du solide, résume dans beaucoup 
d’endroits d’autres ouvrages (tels que la « dyramique des 
réseaux cristallins » de Born). Les notions de base sont 
celles d’onde électrorique quantifiée, dars un réseau, et 
celle, beaucoup moins claire, de phoron, Car le phonon 
n’est pas une onde élastique quantifiée, comme le précise 
l’ouvrage actuel. C’est la mécanique quantique ordinaire 
(modèle 1930 et non modèle 1955) qui gouverne l’interac- 
tion et les transferts. L’ouvrage est dore accessible à un 
public, sérieux certes, mais réarmoirs éler du. I. Phorons 
(dynamique des réseaux, chaleurs spécifiques). II. Élec- 
trons (zones de Brillouin, surfaces de Fermi). III. Inter- 
action phonon-phonon (forces anharmoniques, expansion 
thermale). IV. Interaction électron-électron. V. Interaction 
électron-phonon (discussion du principe adiabatique.) 
VI. Diffraction par les imperfectiors du réseau. VII. Théo- 
rie formelle du transport. VIII. Conduction d’un réseau. 
IX. Conduction électronique dans les métaux (conductivité 
électrique et thermique). X. Mobilité dans les semi- 
conducteurs. XI. Effets superficiels et effets de dimension. 
XII. Phénomènes de transport dans les champs magné- 
tiques. 

J. WINTER. 
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ABSORPTION NUCLÉAIRE DES PHOTONS PAR l?C ET *’Al 


Par G. TAMAS, J. MILLER, C. G. SCHUHL et C. TZARA, 
Section de Physique Nucléaire à Moyenne Énergie, C. E. N., Saclay. 


Résumé. — L’absorption nucléaire des photons par 1?C et ?7Al a été mesurée par transmission 
à l’aide d’un spectromètre à effet Compton. Les résultats obtenus sont comparés aux prévisions 
de diverses règles de somme. Certains faits intéressants semblent se dégager de cette comparaison 


et il importe d'accroître la précision des mesures. 


Abstract. — À Compton spectrometer has been used for measuring by transmission the total 


absorption nuclear cross-section of 1?C and ?7AI. Its Lio 
Some interesting facts appear, but it is necessary to increase the precision 


of various sum rules. 
of the measurements. 


La section efficace totale d’absorption nucléaire 
des photons est encore mal connue et les résultats 
actuels, surtout pour les noyaux légers, sont sou- 
vent en désaccord. Cette grandeur est intéressante 
car elle intervient dans les règles de somme où elle 
est reliée à des caractères assez généraux de la 
matière nucléaire, par exemple la proportion des 
forces d’échange entre nucléons. Moyennant cer- 
taines hypothèses (transitions dipolaires seulement, 
forme de l’opérateur d'interaction non relativiste), 
la connaissance de la densité locale des états excités, 
de la forme de leur fonction d’onde et de celle de 
l’état fondamental détermine la probabilité d’ab- 
sorption totale. Donc il est plus facile d'interpréter 
les résultats de l’absorption totale que ceux de 
chacun des processus qui la composent. 

La section efficace d'absorption nucléaire des 
photons peut être mesurée en étudiant tous Les pro- 
cessus possibles (y, y), (v, v’}, (r, p) (, n) … à 
l’énergie considérée ou par des mesures de trans- 
mission [1]. Cette dernière méthode offre un avan- 
tage considérable sur toutes les autres : il suffit de 
connaître la densité superficielle des noyaux de 
l'échantillon pour obtenir une valeur absolue de la 
section efficace d'absorption nucléaire et électro- 
nique. Cette dernière doit être soustraite de la 
valeur mesurée de la section efficace d’absorption. 
La section efficace électronique est en principe cal- 
culable exactement. Mais des difficultés pratiques 
surgissent si on veut tenir compte exactement de 
l'effet d’écran pour les processus de création de 
paires dans le champ du noyau ou des électrons, des 
corrections radiatives aux divers processus d’ab- 
sorption (*). Pour l’absorption électronique nous 


(*) Le D'H. Olsen a obtenu des résultats prometteurs en 
calculant les corrections radiatives au rayonnement de 
freinage. Ces travaux permettent d'envisager que, dans 
un avenir pas trop éloigné, les coefficients d'absorption 
électronique seront calculés avec une précision de l’ordre 
de 10—*% (Communication privée du D' L. G. Maximon), 


The results are compared with the predictions 


“ 


utilisons les résultats numériques de G. White [2]. 
Une vérification expérimentale de l’absorption cal- 
culée est possible dans la région où l’absorption 
nucléaire est nulle (en dessous des seuils d’émission 
de particules). Au-delà, il faut faire confiance au 
calcul. 


Dispositif expérimental. — La source de photons 
que nous utilisons est constituée par le rayonnement 
de freinage des électrons de l’accélérateur linéaire 
de Saclay dans 0,6 mm de platine. Les électrons, 
dont l’énergie moyenne est de 32 MeV, avec une 
largeur à mi-hauteur de 10 %, sont répartis en 
500 impulsions de 2 us chacune. L’intensité 
moyenne du courant d’électrons est 10 mA. Le 
faisceau de photons est collimaté, comme l’indique 
la figure 1, avant l’absorbant, par 5 em de plomb 
et 5 cm de bismuth et, après l’absorbant, par 
80 em de plomb. Un capteur d’électrons en alumi- 
nium de 4 em d’épaisseur, situé avant le premier 
collimateur, absorbe les électrons émergeant du 
platine et mesure l'intensité du courant d’électrons ; 
une chambre d’ionisation à parois épaisses de gra- 
phite [3] contrôle le faisceau de photons à la sortie 
du dispositif expérimental. 

L’échantillon absorbant peut être escamoté par 
un système tournant télécommandé qui permet de 
passer rapidement d’une mesure directe sans absor- 
bant à une mesure avec absorbant sans arrêter le 
faisceau de l’accélérateur. Le carbone a 700 mm de 
long et une densité de 1,68. Quant à l'aluminium, 
sa longueur est 350 mm et sa densité 2,7. 

Le spectromètre Compton, inspiré d’un modèle 
russe de petites dimensions [4], analyse, par un 
champ magnétique, les électrons émis par effet 
Compton vers l’avant et détectés à leur passage à 
travers les détecteurs A et B (voir fig. 1 et 3). I] 
fera l’objet d’une description ultérieure plus 
poussée. 


Rappelons que l’énergie des photons incidents Æ 
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\source 0.6 mm. Pt 


Fic. 1. — Schéma du dispositif expérimental 


est reliée à l’impulsion p de l’électron émis par la 
relation valable à £, > mc? : 
mMoc? 

se 

Pour un photon d'énergie £ tombant sur le ra- 
diateur où se produit l’effet Compton, le nombre 


d'électrons qui passent par les détecteurs A et B 
est donné par la relation 


EF, cp — 


RE) AE We Fit aps ae 
4 8Tpo 2Y +1 2u + (emly) 


où N°: nombre d’Avogadro, 

Z : numéro atomique du constituant du ra- 
diateur, 

À: nombre de masse du constituant du ra- 
diateur, 

ro : rayon classique de l’électron, 

VOTE Tree, 

m  : épaisseur du radiateur en g/cem?, 

z : perte d’énerge de l’électron dans le ra- 
diateur, en m,c?/g 

2e  : hauteur du radiateur, 

21 : largeur du radiateur, 

eo : rayon de la trajectoire moyenne dans le 
spectromètre. ; 


2ue, : largeur des détecteurs A et B. 
La réponse spectrale théorique d’une raie mono- 


Em 


——— 


E 


r(-2-En) KA) se 
So 
Fic. 2. — Réponse spectrale d’une raie monochromatique 


du spectromètre Compton. 


chromatique est représentée par la figure 2 et Ja 
résolution relative en énergie sera 
s — AEIE = p + (emly). 

Ces relations sont valables si (Z2/A)m < 240, 
c’est-à-dire si la distribution angulaire des électrons 
émis est pilotée par l’effet Compton plutôt que par 
la dispersion angulaire due à la diffusion multiple 
dans le radiateur. Les conditions de l’expérience 
sont les suivantes : 


m — 50 mg de polystyrène 
e— 10" mm 
PL 25 mm 


£o — 400 mm 
200; = 4: mm, 


Dans ces conditions, et en tenant compte de la 
diffusion multiple dans le premier détecteur nous 
obtenons : 


R(y) = 6,4 X 103 


î L l 1 
1 + (30/y) (: "4 + 7.) 
et 
s — 5 x 10—35[1 + (60/Y)]. 


Pour F=120,2 MeVansi=t,22 100% 

Il faut en plus tenir compte du rendement de 
détection qui réduit sensiblement A(y). Pour des 
mesures de transmission la connaissance du rende- 
ment n’est pas nécessaire, seule importe la stabilité 
dans le temps. Nous avons tenu à en donner lex- 
pression afin de concrétiser les possibilités instru- 
mentales du spectromètre. 

La détection du passage des électrons en À et B 
(fig. 1 et 3) s’effectue par deux scintillateurs plas- 
tiques de type Ne102, le premier de 0,2 mm d’épais- 
seur, le deuxième de 2 mm. Deux photomultipli- 
cateurs 56 AVP collectent les photons par l’intermé- 
diaire de deux conduits de lumière en plexiglass sur 
lesquels les scintillateurs ont été collés. Les électrons 
qui circulent dans la chambre sont convenablement 
collimatés, en particulier pour éviter l’effet Ceren- 
kov dans le plexiglass. 

Les impulsions des photomultiplicateurs atta- 
quent deux circuits à coïncidence Cet C, placés en 
parallèle (fig. 3). Sur le circuit C, l’une des impul- 


COÏNCIDENCE 2 


ras CHELLE 2 


RETARD 5:10 


DM.A 
ÉOHTÉTATEES è SEE ES — . 
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Fic. 3. —{Schéma de la détection. 


sions est retardée de 5.108 s de sorte que le cir- 
cuit C, donne le nombre des coïncidences fortuites. 
Les coïncidences sont enregistrées dans deux 
échelles de comptage rapide de temps de réso- 
lution 0,1 us étudiées et construites par les soins 
de M. Mey du Département d’Électronique. Les 
circuits à coïncidence sont du type Graham Bell et 
Petsch (*). 

La largeur à mi-hauteur de la courbe de réso- 
lution de ces circuits est 2r — 1078 s. 


(*) Nous tenons à remercier M. Armbruster de l’Institut 
de Recherches Nucléaires de Strasbourg, qui nous a fourni 
un prototype de ce circuit, 
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L’induction magnétique dans le spectromètre est 
mesurée par l’intensité du courant circulant dans 
ses bobines en suivant. toujours le même cycle 
d’hystérésis et a été étalonnée par la résonance 
nucléaire si bien que l’énergie est définie à mieux 


que 1 %. 


Conduite des mesures. — Soit V5 le nombre 
d'électrons comptés dans la voie sans retard, N6 le 
nombre d’électrons comptés dans la voie avec 
retard, NV le nombre de positons, sans retard 
(obtenu en inversant le champ magnétique du spec- 
tromètre), V5 le nombre de positons, avec retard, 
sans absorbant dans le faisceau, et N—, NT, N*, 
N'T les mêmes nombres lorsque lPabsorbant est 
inséré dans le faisceau, le spectre incident restant 
le même. 

La section efficace d'absorption est alors donnée 
par 


À 1 IN 
OGtot — N el Log N 
où 
NN EN) y) 
et 


N = (N——N—) —(N+ — N'+) 


Si l’on exprime ot en millibarn, 
pour le carbone-12 


Otot — (169 + 0,4) Log (N,IN), 


pour l’aluminium-27 


Otot —= (474 + 1) Log (N,IN). 


Chaque point a été repris plusieurs fois. Ainsi, 


pour le carbone-12, nous avons fait 5 mesures en 


moyenne à chaque énergie et pour l’aluminium-27, 
4 mesures. Les différents résultats étaient bien 
reproductibles entre eux, aux incertitudes satis- 


tiques près. L'erreur relative sur ox due à la géo- 


métrie non idéale est inférieure à 2 °/, et, dans la 
zone où passe le faisceau de photons, la variation 
relative de la densité du graphite est inférieure 
à 1 % ce qui n’affecte pratiquement pas le coeffi- 
client d'absorption. Nous avons constaté l’existence 
d’une perte de comptage croissant avec l'intensité 
incidente et en avons tenu compte. 

Les erreurs systématiques proviennent, d’une 
part, des variations de l’intensité et de l'énergie du 
faisceau d'électrons et de l’évolution de l’élec- 
tronique entre la mesure directe et la mesure avec 
l’absorbant interposé, et, d’autre part, de l’incer- 
titude (1 %) sur le champ magnétique due à l’hys- 
térésis entre deux mesures sur un même point. 
L’intensité de la chambre d’ionisation varie de 
moins de + 1,25.10—2 autour d’une valeur moyenne 
et celle du courant d'électrons de + 102, Quant à 
l'électronique, ses variations au cours d’une mesure 


LE 


N° 6 


entrainent des variations de comptage inférieures 
à 1072. Nous représenterons l’ensemble de ces 
mesures par A’(o+). Ces erreurs sont ajoutées qua- 
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dratiquement aux erreurs statistiques. En utilisant 
les notations précédentes, nous avons pris comme 
erreur statistique sur oct, l'expression suivante : 


A 1 
Astatist. (o tot) 


— — + 1+ 
9 (Nabs UE Nabs SE Nabs SE Nas) 


_ Wei NÉ 


V2 CERTES 


1e 2 
N àbs 


En dehors des incertitudes énumérées ci-dessus, 
il subsiste encore une incertitude systématique qui 
provient, d’une part, de la correction de perte de 


comptage, d'autre part, des erreurs de mesure sur 
la densité de l’absorbant. Pour le carbone, cette 


incertitude est de T mb et pour l'aluminium, 


de 10 mb. 


Résultats. — Sur les figures 4 et 5, nous avons 
représenté, pour les deux éléments étudiés, la sec- 
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Fic. 4. — Sections efficaces d'absorption 
de photons pour l?C. 


tion efficace d'absorption totale mesurée et la sec- 
tion efficace d'absorption électronique calculée [2] 
ainsi que la section efficace d'absorption nucléaire 


Onucl — Otot — el. 


Les barres d’erreur représentent A(o+) telle que 


A(Gtot) — V'[A/(otot)]? + [Astat(otot)]?. 


L’incertitude systématique signalée à la fin du 
paragraphe précédent n’est pas figurée. De plus, 
POUF Ga, 11 Y à une incertitude supplémentaire qui 
vient de ce que la section efficace d’absorption élec- 
tronique n’est pas calculée à mieux de 1 à 2% près. 

Dans la zone où il n’y a pas (ou peu) d'absorption 


=. log Nadir. 
mr Jÿt "Sés 
mb 


E MeV 


c c, 


nucl.= rotal.— Sa. 


120 


Fic, 5. — Sections efficaces d’absorplion 


de photons pour ?7AL 


nucléaire, nous constatons que les points expéri- 
mentaux sont plutôt au-dessous de la courbe théo- 
rique. La méthode de moindres carrés nous donne 
une correction : 


: es ; 1e 
Ac = — Au = — 4,3 + 2,6 mb pour (@ 


AGei CEE AGuuel — — 14 + 9 mb 


pour ?7Al. 


936 


Remarquons que Aou fou © 1,5 %. 

Cette valeur trop petite de 6 peut être due à une 
erreur systématique d’origine non décelée, ou bien 
à ce que a est plus petite que la valeur calculée. 
Dans un domaine d’énergie ne s’étendant pas trop 
au-delà du début de la « résonance géante », nous 


COrrigerONnSs Gnna en rajoutant la valeur de (-- Aa) . 


calculée ci-dessus. 

En ce qui concerne la forme de cu, 11 semble que 
l'absorption ne se localise pas dans une région 
d’énergie particulièrement étroite et qu’au-delà de 
la « résonance géante » vue par de nombreux 
auteurs, les photons sont encore absorbés d’une 
manière appréciable. Pour le carbone-12, dans la 
région de la «résonance géante » (20-23 MeV), nous 
trouvons une section efficace d’environ 25 mb, ce 
qui est en accord avec les résultats trouvés jusqu’à 
présent. Quant à l’aluminium-27, il semble que la 
courbe présente plusieurs maximums comme dans 
les résultats de Mihaiïlovie et coll. et ceux de 
Ziegler [1]. Le maximum d’absorption nucléaire, 
de 80 mb, est atteint à une énergie de 22,5 MeV. 


Diseussion. — Vers 20-25 MeV, la résolution 
instrumentale est de 200 à 300 keV, mais les 
mesures sont distañtes d’environ 1 MeV. Nous 
ferons l'hypothèse que omna Varie peu rapidement 
dans cet intervalle ainsi que les réactions (y, n) le 
montrent. 

Dans ces conditions, pour le carbone : 


amie ont iet 
four reilÆ) — (200 + 35) MeV.mb. 


La règle d'absorption dipolaire de Levinger et 


Bethe 
(ee) 27°? e? ñc NZ 
pe î Onna(E) dE = 


et celle de Gell-Mann, Goldberger et Thirring [5] 
jusqu’au seuil de photoproduction de mésons 


: ue 2x? e? Ac 2 
56 = ff cue(E) dE = ER LE (1 + 01 ) 
donnent le même résultat si x — 0,5 

6o = 60 — 250 MeV.mb. 
Pour l’aluminium, nous obtenons 


(4 + 0,8x) 


28MeV < 
jou GauallÆ) dE = (580 + 100) MeV. mb, 


alors que six = 0,5 
6 66 — 570 MeV.mb. 
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Les valeurs expérimentales des sections efficaces 
intégrées sont, ou très voisines ou supérieures aux 
sections efficaces intégrées théoriques. Des mesures 
ultérieures plus précises devront être faites avant 
de conclure sur la validité de la règle de G.G.T. 
(qui devrait être une limite supérieure mais où le 
coefficient 0,1 (A2/NZ) n’est pas très sûr), ou de 
ER se 

La règle de somme 


LR POSER REP EAN 
ia dE 


É* < R3 
5. AA Rÿ > 


est valable, comme l’a remarqué Foldy [6], si la 
fonction d’onde de l’état fondamental du noyau 
(A, Z) est symétrique en les coordonnées des 
neutrons et des protons. Elle permet de trou- 
ver < R2 >, Si < R£ > est très différent des 
valeurs données par les mesures de distributions 
des charges par diffusien élastique d’électrons, nous 
pouvons être assurés que cette condition n’est pas 
remplie. Remarquons que o& test peu dépendant 
de ce qui se passe à haute énergie. 
On trouve 


cm (120) 28 EM SE D 


o—1(27Al) — 31 + 7 mb, 
d’où 
VER? D 9 = (1,7 + 0,15) 10=1#cm, 
VE RE y = (2:52 + 0:25) 40—15 em, 


alors que les résultats de la diffusion élastique des 
électrons donneraient (1,75 -- 0,02) 1071 cm pour 
le carbone et (2,5 0,03) 10% cm pour l’alumi- 
nium en prenant une densité uniforme et un 
rayon r, égal à 1,2.1071% cm, et après soustraction 
de l’effet de la distribution de charge du proton. 

La condition de Foldy, rappelée ci-dessus, semble 
remplie pour ?C et ?7AI à la précision près de nos 
résultats. 

Nous tenons à remercier M. J. Thirion pour l’in- 
térêt qu'il à porté à ee travail, MM. Barreau et 
Malé pour la fabrication et la mise au point des 
divers appareillages utilisés et l’équipe de gestion 
de l’accélérateur linéaire de Saclay. 


Manuscrit reçu le 19 janvier 1960. 
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LES GROUPES DE TRANSFORMATION ET LES ÉTATS PROPRES 
DANS LES THÉORIES DU NEUTRINO 


Par K. H. TZOU, 


Institut Henri-Poincaré, Paris. 


Résumé. — On identifie les états propres d’une particule de masse propre nulle de spin 4/2 
engendrés par le groupe d’inversion & dans la théorie à quatre composantes ainsi que dans celle à 
deux composantes. On tiendra compte de la constante de mouvement y;, la chiralité, qui est 
étroitement liée à l’invariance dekjauge M. D'autre part, du fait d’un arbitraire relatif au choix 
de la définition des solutions dans le cas de masse propre nulle, la conjugaison de charge donne 
lieu à l’invariance de jauge C en cas de particule neutre. On démontre que l’invariance de jauge M 
est une invariance parfaite, mais l’invariance de jauge C ne l’est pas, si l’on précise les nombres 
quantiques du moment cinétique, de l’hélicité et de la chiralité. 


Abstract. — In the case of a particle of spin 1/2 and of zero rest mass, we identify the proper 
states generated by the inversion group $ in the four-component as well as in the two-component 
theory. The constant of motion y;, chirality, whichis intimately connected to gauge invariance M, 
is taken into consideration. On the other hand, owing to an arbitrariness in the choice of the 
definition of solutions in the case of zero rest mass, charge conjugation leads to gauge invariance C, 
when a neutral particle is considered. It is shown that gauge invariance M is a perfect invariance, 
but gauge invariance C is not, if one takes into account angular momentum, helicity and chirality 


quantum numbers. 


1. Introduction. Les invariances de jauge. — 
Dans deux articles précédents [1], [2], nous avons 
étudié en détail des conséquences de l’invariance 
de la théorie de l’électron de Pirac par rapport au 
groupe de Lorentz généralisé comportant les inver- 
sions P, T, la conjugaison de charge Cet l’inversion 
de masse M. Nous y avons établi, entre les états 
propres d’un électron, toutes les relations engen- 
drées par les transformations de ce groupe, en par- 
ticulier par celles du groupe d’inversion. 


G—{1, P,T, M, C, et combinaisons 1. 


La masse propre nulle constitue toujours un cas 
singulier dans les théories des particules et nécessite 
une formulation particulière qui n’est souvent pas 
un simple cas spécial des théories générales pour 
masse propre non nulle. Dans la théorie du spin 1/2, 
une particule de masse propre nulle doit obéir à 
l’équation de Dirac suivante : (avec # — c — 1) 


(de + ieAbr, t)) F{e, r, 1) 
no (De ced(r t) le, rt} 0, (1) 


OÙ à, — d/0t, d, — ddr, etc., et «y, Gr, Ga Sont 
trois matrices 4-4 anticommutatives. Il s’agit ici 
d’une théorie à quatre composantes. Le manque 
d’un terme (absence du terme de masse) dans cette 
équation crée, au point de vue de la représentation 
du groupe de Lorentz, un certain arbitraire qui se 


manifeste dans deux invariances de jauge de la 


théorie. La première, que nous appellerons inva- 
riance de jauge M, est liée directement à la masse 


nulle. En fait, l'équation d’ondes (1) est invariante 
par rapport à la transformation 


F(e,r, t) > eme P(e,», 1), (2) 


5 = Mets 


La deuxième, qu’on appellera invariance de 
de jauge C n’est valable qu’en cas de particule 
neutre, neutrino par exemple. Dans ce cas, l’équa- 
tion d’ondes (1) est invariante par rapport à la 
transformation (*) 


Fr, i) = Qc FX(r, 1), (3) 
Q, étant une matrice 4-4 déterminée par 
Qc a* O0 — ©. 


Du fait de l’absence du terme de masse dans 
l’équation de Dirac, la théorie à quatre compo- 
santes est une représentation réductible du groupe 
de Lorentz. Elle peut se réduire en fait en deux 
équations identiques et irréductibles, ayant cha- 
cune deux composantes indépendantes. La réduc- 


(*) L’invariance par rapport à ce type de transformation 
est évidemment valable aussi dans le cas d’une particule 
de masse propre non nulle. Mais C engendre toujours un 
changement d’état d'énergie (+ <—> —) [1], on ne peut 
donc pas, à ce propos, parler de l’invariance de jauge 
quand la masse propre n’est pas nulle. Lorsque la masse pro- 
pre est nulle, ilexiste un certain arbitraire dans le choix de la 
définition des solutions. Grâce à cela, la transformation (3) 
peut bien être considérée comme une transformation de 
jauge. Pour le détail, voir les paragraphes suivants. 
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tion se réalise justement en raison de l’invariance 
de jauge ÆZ, car on peut décomposer 


1 
LE - = (4 & Y5); 


où n + Ÿ, comme Ÿ elle-même, satisfont aussi à 
l’équation d’ondes (1), mais n’ont chacune que 
deux composantes indépendantes. 

Selon (4), la réduction de la théorie à quatre 
composantes conduit à la théorie à deux compo- 
santes. L’équation d’ondes dans cette dernière 
théorie est 


(de FE ieA(r, t)) ÿe, , t) 
+ s.(d — iæA(r, t)) ble, r,t) = 0, (5) 


OÙ S4, S2, S4 sont trois matrices 2-2 anticommu- 
tatives et 4 est un spineur à deux composantes. On 
sait que cette théorie irréductible à deux compo- 
santes n’est pas une représentation complète du 
groupe de Lorentz généralisé, car elle n’est pas inva- 
riante par rapport à P et C séparément. Évidem- 
ment, cette théorie ne comporte plus d’arbitraire 
lié aux deux invariances de jauge M et C comme 
dans la théorie à quatre composantes. 

& Dans cette note, nous étudierons, entre les états 
propres d’une particule de spin 1/2 de masse propre 
nulle, les relations engendrées par le groupe de 
Lorentz, en particulier par le groupe d’inversion G. 
En raison de la masse nulle, l’opération M perd sa 
signification ; elle engendre l’invariance de jauge M 
dans la théorie à quatre composantes. 

Malgré leur caractère de jauge, il sera intéressant 
d’examiner si les transformations (2) et (3) engen- 
drent, d’un état propre quelconque, encore le même 
état. Nous verrons plus loin que c’est bien le cas 
avec l’invariance de jauge M, mais non avec l’inva- 
riance de jauge C. 


2. Solutions engendrées par le groupe d’inver- 
sion GS. — Lorsque la masse propre est nulle, y, est 
une constante de mouvement. Il sera cependant 
utile de prendre encore en considération l’opéra- 
tion M pour examiner la relation des états propres 
sous l’invariance de jauge M. 

Les propriétés fonctionnelles des Au(x) étant 
prescrites dans [1] et [2], les solutions engendrées 
d’une certaine solution Y(e, r, t) par P, T, Cet M 
sont, dans la théorie à quatre composantes, 


(P) Y{e, r,t) = Qr Y{e, — r,t), (6a) 
(T) Ve, r, t) = Qr YX(e, r, — t), (6b) 
(C) Fe, r, 1) = Qc F*(— 6, r, 1), (6c) 
(M) Fe, r,t) = Qu Y{e, rt). (6d) 


L'opération (W) est définie ici comme une simple 
multiplication par y. Les Q sont des matrices uni- 
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taires et hermitiennes déterminées par les condi- 
tions 


Qp a Op = — a, (7a) 
Qr a* QT = — 9, (76) 
Qc a* Qc = 0, (70) 
Qu a Op = 0. (7d) 


Même à part une phase arbitraire, ces matrices Q ne 
sont déterminées qu’au facteur e"* près, car toutes 
les conditions (7) sont invariantes par rapport à 


A — eivsE Q. 


La théorie à deux composantes est invariante 
sous CP,T et CPT seulement. Les solutions engen- 
drées par ces opérateurs sont 


(1) de, r, t) = or de, r, t), (8a) 

(CP) de, r, t) HO CP d*(— Ce 1 t), (8b) 

(T) de, P,, 1) —IOT g*(e, T, EE )) (8c) 

(CPT) Ÿ{e, r,1) = ooPr dE, tr, — 1)" 1(64) 
Les matrices w satisfont aux conditions 

OL CPT S OLGPT = S, (9a) 

&GP",T S* OCP,T — —S. (9b) 


Pour l'identification des solutions (6) et (8) en- 
gendrées parle groupe d’inversion G,1l faut connaître 
d’abord les constantes de mouvement. Nous étu- 
dierons à ce propos aux paragraphes suivants deux 
problèmes simples : onde plane et onde sphérique 
monochromatiques. 


3. Onde plane monochromatique. — Dans le cas 
d’une onde plane monochromatique, l’héhcité 
69 — 6.p/pest une constante de mouvement. Les 
fonctions propres de l’énergie, de la quantité de 
mouvement et de l’hélicité peuvent s’écrire 


FÆ)(6)(p A u(+)(s){p) ep.rFpt), (10) 


où $ — + 1 sont les valeurs propres de 64. Les 
fonctions w satisfont aux équations 


cp UE (p) = + u ES (p), (11) 

6p UENS(p) = su HO(p), (12) 
Où & — &.p/p. [l est évident que 

uH) 9) (p) & F9 (— p), (13) 


% désignant l’égalité à la constante de normalisa- 
tion près. 

Lorsque la masse est nulle, y, est toujours une 
constante de mouvement (chiralité). Dans le cas 
d’une onde plane monochromatique, les fonctions 
propres de l’énergie, de la quantité de mouvement 
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et de l’hélicité sont aussi fonctions propres de la 
chiralité. En fait, 


OC Xe, Ys = — La Go Os, 
Y5O0 — à. 


On démontre alors que 


Ys u(E)(s)(p) = + su(Æ)(s)(p), 


Y5 u(Æ)(8)(— P) = + su(Æ)(s)(— P). 


(14) 
(15) 


. Soient 4 — + 1 les valeurs propres de la chira- 
lité. k et s ne sont pas indépendants, mais 


LE EE, 


où £ désigne le signe de l’énergie et € celui de la 
quantité de mouvement. Ainsi, dans les états 
(+ p, + pet (— p,— p), hélicité et chiralité sont 
soit toutes deux positives soit toutes deux néga- 
tives, mais dans les états (+ p,—p}et (— p, + p), 
elles diffèrent de signe. 

Il est à noter qu’en cas de masse propre nulle, 
il existe un certain arbitraire dans le choix de la 
définition des solutions. En effet, par exemple, 
les fonctions propres 


S H6)(p : », 1) = HO (p) e—p.rFr (16) 


représentent les mêmes ondes planes que les 
F6 (p ; r, t). Ainsi, au point de vue physique, 
on peut identifier 


Ÿ Op 0) PE (p; rt). (17) 
Mais, selon (13), 
Ÿ (6) (p 1,1) & FPE pp; r, it). (18) 
Alors 
PH (p ; 7,1) PO p; r, 1). (19) 


Cette relation nous montre que les ondes planes 
(+ p, + p,s)et (— p, — p,s) par exemple sont, 
au signe de phase près, représentées par une même 
fonction propre. 

Quant à l'identification des solutions (6), il est 
facile d’en établir les règles. Si l’on tient toujours 
compte du signe de phase dans la définition des 
ondes planes, ces règles sont les suivantes : (avec 


sS——5,k— — k) 


P: EE, p—>—p,s—>s, donc k—>k; 


T: E—E, D D ss, donc k—k ; 


CINE E, pp, ss, donc k—k; 


M(= y;)n’engendre aucun changement, comme 
nous montrent les équations (14) et (15). 

On remarque que ces règles d'identification sont 
exactement les mêmes que celles dans le cas de 
masse propre non nulle. (Cf. [1] et [2].) 
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En tenant compte de l’identité (19), on voit que 
(C) FE (p ; r, 1) & PPO(— D; », 0) 
= FEHIG(p ; r, 1). 

La transformation de jauge C' engendre ainsi un 
changement d’hélicité (s —>s) et donc un change- 
ment de chiralité aussi (x — k). L'invariance de 
jauge C n’est donc pas une invariance parfaite, si 
l’on tient compte de l’hélicité et de la chiralité. 
Mais, comme y; n’engendre aucun changement 
d'état, l’invariance de jauge M, au contraire, est 
une invariance parfaite. 

D’après les règles établies plus haut, on peut 


identifier les ondes planes engendrées par le 
groupe $ dans le schéma suivant : 


(D) FHO (D ; », 0 & FH; », à) 
2 (CPT) FO (p ; », 0 À PO p; rt); (200) 
(P) FO (p ; », 9 À PHO( p; », 0 


2 (CT) FH (p ; »,9 À FAO; rt); (206) 

(TI FÆO(p ; », 0 À FEB p; r, 0 

(CP) FHO(p; 7,0 À FPEAp;r4; (200) 
(PT) FH (p ; r, 0 À FÆG(p; r,0 

2 (C) FH (p ; r, 0 À PPE(—p;r,t. (204) 


Prenant une certaine représentation des ma- 
trices œ@, on pourra écrire explicitement les spi- 
neurs 4Æ9(— p) (voir Appendice), et l’on démon- 
trera facilement que le schéma (20) est bien vérifié. 


4. Onde sphérique. — Dans le cas de masse 
propre non nulle, il y a deux types indépendants de 
solution pour une onde sphérique monochroma- 
tique, à savoir 


DE (q : r,10) = V9; r, 0, pheFit, (21) 


où g est la fréquence de l’onde, et s = + 1 dési- 
gnent les deux types indépendants de solution 
ayant chacun une seule constante de normali- 
sation. Dans l’approximation non relativiste, le 
moment cinétique orbital est une constante de 
mouvement dont la valeur propre est { = j — s/2 
dans la solution D, Ainsi, dans les grandes 
composantes de VF, les fonctions angulaires 
sont Y##12(0, o), et dans les petites composantes, 
elles sont Y7#17(0, +). 

Quand la masse est nulle, les quatre compo- 
santes de V9 sont évidemment d’un même 
ordre de grandeur. s perd done le sens de spin dans 
ce cas. Alors, s — + 1 désigne encore deux types 
de solution, mais ceux-ci ne sont plus indépendants 
à cause de l’invariance de jauge M. En effet, quand 
la masse n’est pas nulle, l’opération M engendre 


un changement de s (ss ——s) [2]. Cela 
signifie qu’en cas de masse nulle, 


ve VE (a: 7, 6, o) & Vin” (g; r, 0, @). 


(22, 
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Comme y, est une constante de mouvement, 1l 
existe les ondes sphériques qui sont fonctions 
propres de y, aussi. Nous les désignons par 


po 0) = UP (g: 7,6, pleTrit, (23) 


où 4 — + 1 sont les valeurs propres de la chiralité. 
On a évidemment 

Un (g 3 r, 0, 9) (1 + Eva) Vin (9 r, 0, 9), 
(s — + 1 ou — 1). Finalement les fonctions spa- 
tiales satisfont aux quatre équations suivantes : 


{— Ads JS ES GISEALU : r, 0, ©) 


—{La;ii+i);m; k} Un (a; r,6,œ) (24) 


où J — — ir X d + o/2 est le moment cinétique 
total. 

Tout comme dans le cas des ondes planes, on a 
ici aussi d’autres choix de la définition des ondes 
sphériques. En fait, par exemple, les fonctions 
propres 


qe Prin Or o)et/2s) 


représentent les mêmes ondes sphériques que les 
DO (q ;r,t). Au point de vue physique, on peut 
donc identifier 


A ACL et ten (26) 
Il est facile de voir que 
UP (g rm — 0,7 + o) À Un P(g; r, 6, @. (27) 
Alors 
Sr SEE ec), (28) 
puis 
pe, 0) EM (q: à). (29) 


Ainsi, en cas de masse nulle, si l’on n’a pas à pré- 
ciser le signe de phase, les ondes sphériques d’éner- 
gie positive et celles d’énergie négative de mêmes g, 
J, m, k sont représentées par les mêmes fonctions 
propres. 

Pour l'identification des ondes sphériques engen- 
drées par le groupe $, les règles en sont déductibles 
de celles des ondes planes du paragraphe 3. Mais, 
selon (22), l’opération M (= Y;) engendre un chan- 
gement de s. Nous avons finalement, pour les ondes 
sphériques, les règles d'identification suivantes : 

P: EE, mm, s—>s, k—k: 
T': EE, mm, S—s, k—k; 
CEE MMS Ne 
ME): ESE mm ss, EX 


Aucune 
de J. 


opération n’engendre un changement 
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En ce qui concerne C, on a, à l’aide aussi de (29), 


a PM: »,0 R FEP(g; »,0 
Ze wHE) 1, . 
œ FE GRR, Ge 


La transformation de jauge C engendre donc les 


changements m — met k —k. Tout comme dans 
le cas des ondes planes, l’invariance de jauge C'n’est 
pas ici non plus une invariance parfaite. Mais l’in- 
variance de jauge M est toujours une invariance 
parfaite. \ 

Enfin, d’après les règles d'identification citées 
plus haut, on peut résumer l'identification des. 
ondes sphériques engendrées par le groupe S dans 
le schéma suivant : 


k 
(1 PNB SM UE Er 0 


= (CPT) Fe Pa: r,0 & Vie (q: r,t); (300) 

2) PP (5: 7,0 À FEW (g; r,0) 
(CT) Pie (a: 0) À PEEND (g: rt); (30h) 

(7) FM »,0 REP (q; »,0) 
(CP) F6; r,0 REP (g; rt); (806) 

Tr: », 0 À EP (0; », à 
= (0) Ph ; 7,0 & EP (G; 7, 0). (304) 


Prenant une certaine représentation des ma- 
trices æ&, on pourra écrire explicitement les fonc- 
tions spatiales WU (9 ; r, 0, o) (voir Appendice), 
et l’on démontrera alors que le schéma (30) est bien 
vérifié. 

Les solutions (21) ne nous intéressent pas dans le 
cas de masse propre nulle, car s n’a pas de signi- 
fication physique ici, et les fonctions (21) ne sont 
pas fonctions propres de y,. Mais, si l’on veut iden- 
tifier les ondes d®(g ; r, t) engendrées par P, 
T, Cet M(= y;), les règles et le schéma en sont 
exactement les mêmes que ceux dans le cas de 
masse propre non nulle. (Voir à ce propos [1]et[2].) 


. 9. Identification des états propres dans la théorie 
à deux composantes. — Pour les ondes planes, on 
peut les définir par exemple soit par 


ENS) (p ; r, 1) — vÉ)S)(p) eip.rF rt), (31) 
soit aussi par 

à HO) (p : r,t) = 06) (p) e—ip.rFpt), (32) 
où $ — + Î sont ies valeurs propres de l’hélicité 


Sp — S.p/p. Du point de vue physique, on peut 
identifier 


DHO(p ; », 1 Hp : r, #). (33) 


N° 6 
Les fonctions v satisfont aux équations 


{Sp ; 9 }0ÆO(p) = {4 p;s}vtO(p). (34) 


Il est évident que l’hélicité est positive dans les 
états d’énergie positive, et négative dans les états 
d’énergie négative. Les équations (34) nous mon- 
trent que 


v(Æ)(s)(p) R VF) (— p), (35) 
Alors 
HO (p ; », 0 À PO p; r,9 
YPO p; 7, ti). (36) 


Les ondes planes (+ p, + p,s) et (— p,—p,s) 
par exemple sont donc représentées par une même 
fonction propre, si l’on ne précise pas le signe de la 
phase dans la définition de ces ondes. 

Bien qu’il n’y ait pas d’arbitraire lié aux inva- 
riances de jauge M et C dans la théorie à deux 
composantes, le choix de la définition des ondes 
planes (4 ou d par exemple) reste arbitraire. 

Pour l'identification des solutions (8), les règles 
en sont déductibles de celles établies au para- 
graphe 3 pour les ondes planes dans la théorie à 
quatre composantes, à savoir 

CP: E—=—E, p — p, 
OS DCS 


SES 
P—=—P, ss —S$. 


Puis, à l’aide aussi de (36), on établit le schéma 
suivant : 


(1) LEO (p ; r, 1) & LES (p ; r, 1) 


2 (CPT) LE (p ; »,1) À YPO p; 7,1); (870) 
(T) YO (p ; », 0 À LE p; r,1 
2 (CP) SH (p; r,0 À YPE(p; r, 0. (376) 


Quant aux ondes sphériques, il y a aussi des 


définitions différentes mais équivalentes, par 
exemple 
jm (a 3 7,0) = Wim (g 3 7, 0, ) eFit, (38) 


ef 


dim (gi r, 0 — mm (gi r,nm— 0,7 + ei, (29) 


Alors 


AU) 


(a: v, 0) > din (gi 7,0). (40) 


Les fonctions spatiales w satisfont aux équations 
—{+q;j0+0:;m}wm(ssr 0, e) (41) 


où J — — ir X do + s/2 est le moment cinétique 
total. On démontre facilement que 


ag: r,r—0,x + +) win (g:; r, 0, +). (42) 
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Alors 
(+) A 
QU (gi r, 0 Rd (gr, 0) 


2 Qim(g; , td). (43) 


Donc, si l’on n’a pas à préciser le signe de la 
phase, les ondes sphériques (+ q, jm) et (— q, jm) 
sont représentées par une même fonction propre. 

Selon les règles d'identification établies au para- 
graphe précédent, celles pour les solutions (8) sont 
maintenant les suivantes : 


(CUS GUESS 
TE D, 


E,j—=]j,m—m; 
1), Mm— m. 
On peut alors établir le schéma suivant : 
(+) 
jm | 


(1) dé (g: r, 1) & g;r,i 


= (CPT) We (g: r,0 À Y(g;: r,t); (ka) 
(T) (a; », 4 & QE (g; 5) 
(CP) Ym(asr or bg; rt:  (4ab) 


6. Conclusion. — Aux paragraphes précédents, 
nous avons établi les règles et les schémas complets 
de l'identification des états propres engendrés parle 
groupe d’inversion $ pour les ondes planes et les 
ondes sphériques monochromatiques dans la théo- 
rie à quatre composantes ainsi que dans celle à 
deux composantes du neutrino. La masse propre 
nulle crée un certain arbitraire dans le choix de la 
définition de ces ondes. A cause de cela, la conju- 
gaison de charge C devient une invariance de jauge 
dans la théorie à quatre composantes. 

Mais cette invariance de jauge C n’est pas une 
invariance parfaite, car elle relie non pas deux 
états propres identiques mais deux ondes planes 
d’hélicité et de chiralité opposées ou deux ondes 
sphériques de moment cinétique et de chiralité 
opposés. Ainsi, sous C, 


WpENeNEe, AP, t) — (91) (p ; 7, t), 
gg rl or (g:; v,t). 


L'invariance de jauge M, au contraire, est une 
invariance parfaite (l’invariance y;). 

Si l’on ne précise pas le signe de la phase dans la 
définition des ondes planes, neutrino et antineu- 
trino peuvent être identiques dans la théorie à 
quatre composantes, si l’on ne tient pas compte de 
l’hélicité et de la chiralité. En fait, remplissant 
l'hypothèse de Majorana, les ondes planes sont les 
combinaisons 


pm, Tv t) de (0), 


p;T, 1). 


En précisant le signe de la phase dans la défi- 
nition des ondes planes, nous avons conclu, dans 
un article précédent [1], que les ondes planes en 
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théorie de Majorana ne sont pas des états propres de 
l’énergie ni de la quantité de mouvement, ce qui 
n’est pas le cas si l’on n’a pas à préciser le signe 
de la phase en cas de masse propre nulle. 

De même, les ondes sphériques en théorie de 
Majorana doivent être les combinaisons 


ph (2,0 2 FEW (G; v,0, 


qui ne sont pas des états propres de J, ni de y3. 

Dans la théorie à deux composantes, qui n’est pas 
invariante sous C, neutrino et antineutrino se dis- 
tinguent parleur hélicité opposée ; ils ne sont donc 
pas identiques. Dans la théorie à quatre compo- 
santes, un état propre de l’énergie et de la quantité 
de mouvement est en général une superposition des 
deux hélicités opposées, ce qui permet l’identifi- 
cation de l’antineutrino au neutrino selon l’hypo- 
thèse de Majorana en supposant une amplitude 
égale des deux hélicités. 

Dans la théorie à quatre composantes, les ondes 
planes existent aussi, qui sont invariantes par 
rapport à CP, si l’on ne précise pas le signe de la 
phase. Elles sont les combinaisons 


PO (D à», 1) à pe Es) —p:r,i 


quine sont pas états propres dela quantité de mou- 
vement. Elles ne peuvent pas définir des états 
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propres d’un neutrino libre. Il ne convient donc 
pas, en vue de l’identifier au neutrino, de définir 
l’antineutrino avec CP au lieu de C par rapport au 
neutrino dans la théorie à quatre composantes. 

Grâce à l’arbitraire relatif au choix de la défi- 
nition des solutions, l’invariance CPT est une inva- 
riance parfaite dans le cas de masse propre nulle, 
car l’opération CPT n’engendre aucun changement 
d’état propre, comme nous montrent les schémas 
d'identification (20), (30),(37) et (44) dans la théorie 
à quatre composantes ainsi que dans celle à deux 
composantes. Si l’on définissait l’antineutrino 
avec CPT par rapport aù neutrino, ces deux parti- 
cules seraient toujours identiques. Une telle défi- 
nition serait triviale, car elle ne nous donne pas 
d’aspect nouveau des théories du neutrino. 

Dans la présente étude sur l’identification des 
états propres sous le groupe de Lorentz et le groupe 
d’inversion G, nous n’avons considéré que la théorie 
de Dirac des particules de spin 1/2. En fait, cette 
étude et ses conclusions se généralisent facilement 
à toutes les autres théories des particules, car, 
comme la théorie de Dirac, elles constituent aussi 
des représentations du groupe de Lorentz et du 
groupe d’inversion G. Il est d’ailleurs facile de voir 
que les règles pour l’identification des états propres 
doivent être les mêmes dans toutes les théories 
relativistes des particules. 
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Prenons pour @&, o et y; les matrices 


[ e-| 


el 
OS 


=f 1 
WE ol’ 


S1, Sa, 3 étant les matrices 2-2 de Pauli. Les fonctions propres 4#®®{_L p) dansles ondes planes mono- 
chromatiques du paragraphe 3 prennent alors les expressions normées suivantes : (p,; > 0) 


ut+)(+1)(+ 1 (p) = y) (+ D(— p) 


UD (p) = D (— p) 


u—)(+ 1) —1(p) 


ut) 1(+ D (p) = gt) i)(+ D(— p) 


ut++ DD p) — 


P+Ps 
TER she P1 + ia 
2 Vp(p + ps) P+pPs : 
P +) 1Pa 
= (P1 TES iP2) 
LE A UITe eS P +Ps 
2 Vp(p + p:) Pi — Po ; 
PIE DA) 
—pEp:) 
AE! no — (P1 + ip) 
2 Vp(p + ps) P+Ps : 
P1 + 1Pa 


— (P1 — ip:) 
= Ds PEU Pr Ps 
2 Vp(p + pa) | — GES ia) |” 


RAT 7) 
Des La? 
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Puis, dans la théorie à deux composantes ($5), 


DH D(p) = DEL p) = Fr = | . Fe À 
3 1 2 

D 1){(p) = DIU (— p) — = en | 
3 3 


Dans les ondes sphériques du paragraphe 4, les fonctions spatiales sont, à une constante de norma- 
lisation près, 


Te æ à 1/2 en _ 
2 ve (0, o) & Ki (A) Jon Ÿ = PE (0, 9) 


( 15 J r US 
2r +algr) \ 2(] + 1) 2r 
1/2 5 1/2 TE, 
(2) re VÉECE vie, à 86 (2) tan VI 280, 9 
Un (ire) | s 
LIN | = 1/2 EME nee 
Qt (1) J{gr) ee ed Yi78 (0, @) + k (1) Jj+ (gr) VE Ye (0, ç) 
1/2 J+1—m +1/ 
lé l (4) 7 i(gr) VE Tr (6, p) — k 2 Jj+1(gr) Ve Fe (0, ©) 


Ici J(x) est la fonction de Bessel d’ordre impair, et Y?#12(0, ) les fonctions sphériques. 


Dans la théorie à deux composantes, les fonctions spatiales des ondes sphériques ($ 5) sont 


1/2 Te = 1/2 mm ae 
(E)" seen VE rte, à 2 (2) von VE 5e, 9 


De 1 
#ÿm (g ; r, 0, @) = q\u? NES g\1/? Î— M ym+192 
Le (1) Jj+ (gr) V7. YF (6, @) + (2) Ji(gr) = Yi (0, y) 
On sait que 
VRAUE (r — 0, 5 + @) = (— 1)" Ye (0, 9) 
donc 


j k 
1 un 0 (gr, 0, @), 


1) mm (g ; r, 0, gl. 
Manuscrit reçu le 145 décembre 1959. 


Um P(g5r,r— 0,7 + œ) = (— 


mg; r, rm — 0, x + œ) = (— 
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Résumé. — Les résultats de l’article référé comme I sont utilisés pour calculer les tenseurs de 
conductivité électrique et de pouvoir thermoélectrique dans le cas général d’un. cristal possédant 
la symétrie cubique. On montre que les effets magnétoélectriques et thermomagnétoélectriques sont 
insuffisants pour permettre la détermination de la structure de bande d’un cristal à symétrie 


TOME 24, JUIN 1960, PAGE 544. 


cubique. 


Abstract. — The results of the previous paper (I) are used to calculate the conductivity and 


thermoelectric power tensors in the general case of a crystal of cubic symmetry. 


It is shown 


that measurements of the magnetoelectric and thermomagnetoelectric effects are not sufl- 
cient to establish the band structure of a cubic crystal. 


I. Introduetion. — Dans un article précédent(1), 
nous avons calculé la contribution des porteurs 
localisés dans une « vallée » de l’espace réciproque 
à la densité de courant électrique J en présence 
d’un champ électrique E et d’un champ magné- 
tique H. Ce premier résultat nous a permis d’éva- 
luer la densité de courant électrique pour des 
modèles simples, savoir: modèle « S » et 
modèle « G ». Ensuite nous avons calculé la contri- 
bution des porteurs d’une «vallée » unique au flux 
d’énergie isotherme. En utilisant ce résultat, nous 
avons déterminé le tenseur pouvoir thermoélec- 

> 


tique Q (H) pour des modèles simples, savoir : 
modèle isotrope et modèle «S ». 
Nous nous proposons dans le présent article de 
calculer la densité de courant électrique J et le 
ES 


ea 

tenseur pouvoir thermoélectrique Q (I) dans le cas 
d’un cristal de symétrie cubique quelle que soit la 
répartition des vallées. 

Nous utiliserons les notations définies dans l’ar- 
ticle I. 

La discussion des résultats ainsi obtenus nous 
permet de déterminer les paramètres définissant la 
structure de bande, accessibles par la mesure des 
effets magnétoélectriques et thermomagnétoélec- 
triques. 

Nous montrons qu’en toute rigueur, nous ne 
pouvons pas définir totalement la structure de 
bande à partir des résultats précédents, bien que 
certaines approximations que nous préciserons per- 
mettent d’aboutir à des résultats intéressants. 

Pour effectuer nos calculs, nous supposons que le 
cristal considéré possède V vallées distinctes. Ce 
nombre V doit tenir compte du fait que les vallées 
peuvent se trouver en limite de la première zone de 

(*) Article I, J. Physique Rad., 1960, 21, 249. 


Brillouin. Par exemple, V prend les valeurs sui- 
vantes : 

dans le modèle du germanium : V = 4, 

dans le modèle du silicium : V — 6, 

dans un cristal tout à fait quelconque : V — 48 
ou 24. 

Le courant J cherché s’obtient à partir du cou- 
rant fourni par les porteurs d’une vallée en som- 
mant les contributions des différentes vallées : 


J—>;j 
L4 


la somme Y doit être étendue aux V vallées dis- 
2 


tinctes. 

Si V est inférieur au nombre total des symétries 
du cube, soit 48, chaque vallée sera dégénérée 
d'ordre 48/V et par suite, nous pouvons calculer J 
sous la forme suivante : 


J=(V/48)  Zj 


la somme Z étant étendue aux quarante-huit 
vallées que l’on déduit de l’une d’entre elles par 
application des symétries du cube. Nous consi- 
dérerons le groupe G des 48 symétries du cube 
comme produit des deux sous-groupes 


JR — (1, 01, O2 O3 Où Gi, O1 O2 O3) 


et 


R — (ee C3, C?) (1, 2) 


dans lesquels les éléments sont les suivants : 
oi : symétrie par rapport au plan x = 0, 
en rotation de 2r/3 autour de l’axe ternaire 
) 
2 : symétrie par rapport au plan Li — Lo = 0. 
Les axes de coordonnées af sont choisis suivant 
les arêtes de la maille. 
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Notons dès maintenant quelques propriétés 
simples des deux sous-groupes J et 3 : le sous- 
groupe J est commutatif, chaque élément égale 
son inverse et il n’est constitué que des quatre 
éléments (1, 61, 62, 63) et de leurs opposés (— 1, 
M1 109 Os): 

En effet, on a : 


Gi 5j = — 06% OÙ (1, j, k) — permutation de (1, 2, 3) 
et 


C10203 — — À. 


Le sous-groupe € produit des deux sous- 
sroupes JC, — (1, C;, CE) et X, — (1, Z) possède 
_6 éléments distincts. Il est isomorphe au groupe 
des permutations des trois chiffres (1.2.3). 

L'utilisation du sous-groupe * sera simple si 
l’on remarque que : 


Do 


IT. Caleul du vecteur densité de courant. — 
Nous avons montré que pour une vallée le courant 
est : 


Ch CL 


= 
2E ne 2 >> : re 5 
IPC er 


ES 
T 1 27 
RE EN FE EN = M D 

D Es ( ) El 
où n est le nombre d’électrons dans la vallée consi- 
dérée. 

> 

ES 


« : un tenseur du second ordre tel que l’énergie & 
d’un électron s’écrive 


> 
> 
26 = «k.k. 


Rapportée aux axes des ellipsoïdes représentant 
2 


ea 
les surfaces d’énergie constante, la matrice « s’écrit 
par définition : 


— 0 0 
> mi 
Ze il 
6 = jo) — 0 
Ma 
0 0 == 
M3 


où M M2, M, sont les masses effectives principales. 


> 
> 
œ —1 


Per 


— mm; Ma m3 


Es 
Cou — LC. 2 ri Fa est le carré de la pulsation 


cyclotron. , 
Pour évaluer J, nous calculons successivemen 


les termes d’ordre zéro, un et deux en H. 
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a) TERME D'ORDRE 0. — Nous appelons À la 
> 
. , . x Re 
matrice équivalente à « dans les axes du cube, et 
nous posons B — AT, 
Le courant cherché est : 


JT, = (VJa8)n < + > (2, S—' AS) E. 


En utilisant le résultat selon lequel tout tenseur 
du second ordre possédant les symétries du cube 
doit être proportionnel au tenseur unité, il nous 


suffit pour évaluer (2 SL où de calculer l’un des 
sea 


éléments diagonaux. 

Commençons par calculer la contribution du 
sous-groupe X, = (1, C3, C?). En remarquant que 
l'effet de C, sur une matrice quelconque revient à 
permuter circulairement l’ordre des indices, on 
trouve que : 

> 
> 
Z S—1A$S — Tr A — Tra 
Sec: 
La somme cherchée s’écrit alors : 


> T1 | Dr As) qi 


TEJE TS: sefc: 


La quantité Tr À étant invariante par rapport 
aux opérations 7e J + X (au nombre de 16), la 
somme cherchée s'écrit : 


D S—1 AS — 16 Tr. à. 
seG 


Finalement, le terme d’ordre zéro s’écrit : 


> 
> 
en LE 


b) TERME D'ORDRE 1. — Nous devons calculer : 
V n à : 
ARE a —1BSH/N E 
Ju Pire” S SA 
où 
DS 5 (H.$—1 BSH). 


Si nous nous limitons dans J aux termes du 
second ordre en H, l'expression ci-dessus de J, peut 
s’écrire : 


Ponts S (S—1 B$)] H A E. 


M 48 1 se& 


Le résultat précédent nous permet d’écrire 1immé- 
diatement 


_ 
e 
Fe el H \ E. 


c) TERME D'ORDRE 2. — Le terme du second 
ordre en H dans l’expression du courant J est 


donné par : 
Be DR Re 


sq 
8 P 
(H(H.E)—(H.S- BSH)S—' ASE] 


SOIL : 
Vn 
J=— <> H(HLE) 
_ MMS (HS BSH)S-1ASE 
48P sea 
Il nous faut donc évaluer la somme : 


V(H, E) = Ÿ (H.S—-1BSH) S—' ASE. 
SeG 


Le groupe G étant un groupe de transformations 
orthogonales : 
S—1 — ,$0 
où S° est la matrice transposée de S. 
Le vecteur V(H, E) peut alors s’écrire : 


V(H, E) = X (HS BSH) S—1ASE 
SeG 


= X S—1(SH.B SH) AS E. 
se& 


Si nous définissons le vecteur U(H, E) par la 
relation : 
U(H, E) = (H.BH) AE 
le vecteur V(H, E) s’exprimera par : 


V(H, E) — 2 $—1 U(SH, SE). 
Se 


Nous commencerons par calculer le vecteur 
W(JH, E) défini par : 


W{(H, E) = À 6-1 U(oH, 6E) 


oeJe 


= Y oU(oll, 6Ë) 
ee 


Le] 
puis V(H, E) sera donné par la somme : 


V(H, E) = 2 C—1 W(CH, CE). 
ce 


Nous montrons en annexe Î que : 
V{H, E) = LH(E.H) 
+ (p— v) EH? + (3—p— 6) T(E, H) 
où nous posons : 
= 16 À Aij Bij 
5 AE D 
ES 
Ed 
DENNNE ZM TEEKS e ce Re 
v= 82 AùBü 
% 


Si 
ea 
œ—1 


T(E, H) étant un vecteur qui dans les axes du 
cube a pour composantes ÆE, HŸ, E, H5, E, H2. 
At] et Bij sont respectivement les éléments des 
matrices À et B. 
En utilisant la symétrie des matrices À et B nous 
pouvons écrire que : 


ue + 2v = 16 2 Aij Bij = 162 Aij Bji = 48. 
: TR 
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Cette relation nous permet de donner une expres- 
sion de V(H, E) ne dépendant plus que des deux 
paramètres p et v: 

V(H, E) = (48 — 2v) H(E.H) 
+ (p— v) ER? + (5v — p — 48) T(E, H). 

Finalement le terme du second ordre en H 
s'écrit : 

Vn 


se Rs Er 2 
= 5 < > 2vH(EH.) + (v— o) EH 


+ (48 + p— 5) T(E, H)]. 
Si nous remarquons que dans le cristal le nombre 
total d'électrons de conduction est N — Vn, l’ex- 


pression définitive du vecteur densité de courant 
est : 


J; 


RYY 


De. E 


D Ces (2vH(H.E) + (v— +) ER 


se L48P 


+ (48 + pb —5v) T(E, H)]. 


IIT. Effets thermomagnétoélectriques. — Nous 
calculons le tenseur « pouvoir thermoélectrique » à 
partir de sa définition donnée dans l’article (1) : 


> > 
> > > _ > 
à o—'(H)ZKi—H) vu? 
QUE RTE 
ce 
où o(H) : est le tenseur de conductivité élec- 


trique en présence du champ ma- 
gnétique 1. 


Se 

+ 

K;(4H) : le tenseur du second ordre liant le 
flux d’énergie isotherme au champ 
électrique, pour les porteurs de la 
ième vallée, 


u : le potentiel électrochimique. 
#e 
Le tenseur 6 1(H)s/évalue aisément en utilisant 
la formule d’inversion valable au second ordre en H 
démontrée dans l’article I. On obtient le résultat 
sous la forme : 
E = 54H) J = p5J + p, HA J 
+ e2{L(J.H.) H + MJH + NT(J, H)] 
où T(J, H)est un vecteur dont la 1è®e composante 


est dans les axes du cube J, H,2 
et 


3 
CSS 


16, Sir Tr a 


FAR OR tiRE es" MAC) 
ME CT TU 
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3 < T2 > Tr a—1 


“4 NP< 5 >2{Tr a 
DR re 
PRÉPE En UTre 
B) <> 
P2 


16 NP < 7 >? [Tr a) 
N = 48 + p — 5. 


Nous avons calculé (article I) que pour une 
vallée : 
2e > 
WM=K(H)E — — PANES té >aE 
r? 8 1 
7 + o =. AETA r) 
T6 { + « 
nuit. pHU.E) — xE ? 
<e 
Le tenseur K (H) étant formellement identique 


> 


=<e 

au tenseur 6 (H}), la somme ZX K;(— H) étendue 
aux 48 vallées se déduit immédiatement du calcul 
précédant sous la forme : 


DK -H)r7 
UE D MC »s-UHA 
3 3P 
pe) 
ri” (2v H(V.H) + (v— +) VH? 


+ (48 + 0 — 51) T(V, H)] 


où V'est un vecteur quelconque. 
Rappelons que ce résultat n’est valable qu’au 


second ordre en H. 
Æ > 


- > 

Connaissant les tenseurs 6 —{(H) et 2 K;(— H) 

nous pouvons calculer facilement le tenseur « pou- 
> 


> 
voir thermoélectrique » Q (H) sous la forme : 
a = 
G = (Q— À) V 


+ Q HA V + L'H(H.V) + M'HV + N'T(V,H) 


où (en unités c. g. 8.) : 


ne 
De a 
Le k1 Tra—ip< 6 > | 
Dep = Tr 2 
RAR AN 10 Fy (Tr &—1)2 
a a(ri] 
,_k1 1 [v—eb SC ne 
471 = 5m | AC eptre 1] 
, AL 1 48 + p—5v 
DR pre 16 fr) 
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f(r) et g(r) étant deux fonctions assez complexes du 
paramètre + définies par : 


jt) ve io ar 

Te it 
2 Te AT OR 
SON RS 0 00 +. : 


D'autre part, notons que le potentiel électro- 


chimique y sera fonction de la densité d’états de la 
bande à laquelle appartiennent les porteurs consi- 
dérés ainsi que de la densité de porteurs et de la 
statistique utilisée. Par exemple, pour un semi- 
conducteur de type x non dégénéré (statistique de 
Maxwell-Boltzmann), nous avons : 


ŒE — kT log (NM) 
où {V est la densité d’électrons de conduction et M, 
le nombre d’états de la bande de conduction 


2m KT 
h? 


3/2 
We = 2 ( ) (naine mne)UeRTE 


IV. Interprétation des résultats précédents. — 
Dans les paragraphes précédents, nous avons 
> 


> 
montré que le tenseur de conductivité électrique © 
ÆS 


e 72 * > 0 
et le tenseur pouvoir thermoélectrique Q s’expri- 
maient à l’aide des quantités suivantes : 


Tatin Ve Autos Trio lien 
v 


Nous nous proposons maintenant de rattacher 
explicitement ces quantités aux paramètres qui défi- 
nissent la forme des vallées et leur position dans les 
axes du cube. 

Par définition « est une matrice diagonale dont 
les éléments sont les inverses des trois masses effec- 
tives principales m, m, m3. On a donc : 

1 1 1 


MG = ES — 
TA OUI UTITS 


Trait = m1 + M: + M: 


I Te TT Elle | 
+ | | | +- _— 
l l Ï Q 


PUS TT TI S UT TS TS 


TD 
I 


s[3 + 
Pour des raisons de simplicité nous poserons : 


mi 
M3 MTS 


m m m m 
Qu = + — ee 2 + —, 


Ces quantités u; ne dépendent que des rapports 
des masses effectives principales entre elles. Ce sont 
des nombres qui définissent la « forme » des vallées 
d'énergie à une similitude près. [ls ne sont pas indé- 
pendants et vérifient la relation suivante : 


2 Ua ta — (ui + ô + uà) + 1 = 0. 


o48 


Le cas isotrope pour lequel les vallées sont sphé- 
riques correspond à by = Le = Ug = 
A l’aide de ces notations : 


p = 8(3 +22 y). 


Le calcul de v fait intervenir explicitement les 
matrices À et B, précisons-les. Considérons une 
vallée arbitraire, il existe une transformation ortho- 
gonale S qui permet de passer des axes du cristal 
aux axes principaux de cette vallée. Si l’on con- 
naît S on peut déduire l’orientation de la vallée 
considérée, puis, par le jeu des opérations conser- 
vant la symétrie du cube, l’ensemble des 48 vallées. 

La matrice S représente une transformation 
orthogonale, elle doit donc vérifier la relation : 


S9 — ,$—1 
où 5 est la matrice transposée de S, c’est-à-dire 


que ses éléments s°;; sont liés à ceux de la ma- 
trice S par la relation 


0 F 
Sii — Sÿi- 
— > 


E 
Par définition, on a: À = 59 4 $,B—S0% 15, 


c’est-à-dire : 


<0 I 
RE 2. STI = » 
k M 


Aij = Ski Skj 
ki MX k 


== 2 Si Mk dk St — > mA Sti Stj- 

On vérifie les propriétés de symétrie de À et B : 
Aij = Ajü Pi = BP. 

Calculons v : £ 


mu. 2 


= 82 4iBi 82 GE 7) 


C= 


mx 


=sEh ts (ut) (Dh à). 


Mx M 


Posons : 
_ yes 
— d'kili 
"1 


où 7, k, l'représentent des nombres tous différents, 
égaux à 1, 2 ou 3. En remarquant que les relations 
d’orthogonalité permettent d’écrire 


2h 22 0 
î,k j 
il vient : 
8(5—22 vi) + 16 2 u5 y; 
"] 1 


= 24 + 16 À (u;ÿ — 1) y; 


Les paramètres y; que nous venons de définir 
sont Ceux que nous utiliserons pour caractériser la 
Cposition » des différentes vallées. En effet, si l’on 
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se donne arbitrairement trois nombres y; on peut 
essayer d’en déduire la matrice S ; le problème est 
soluble et présente en général 48 solutions dis- 
tinctes qui se déduisent l’une de l’autre par les 
opérations de symétrie du cube. 

La connaissance des y; et des x; est donc équi- 


forme des vallées. Pour obtenir la structure com- 
plète, il faut ajouter à ces paramètres sans dimen- 
sion, une des trois masses effectives principales mi 
(ou leur moyenne géométrique m = (m1 m2 Ma)/3). 

En regroupant les résultats précédents, nous 
écrivons le courant J sous la forme 


J=oE+0 HAE + L(E.H)H + MEH° - NT(E,H) 


avec (en unités c. g. s.) : 


ee ft 


Ms Me Ma 


Ne ne Los Tu + Ma + M3 
3C M1 Mo M3 
Net 1 
[ie PO IPN SRE 
OC M1 Mo Ma É Ë n (es 1) ‘il 
Nef 
M = 3e < T? ba (lui — 1) Jul 
3c Mi Mo My 
_. 
N ——— no D — — = 
M3 Mo M3 [vi sÈ(u 1} 81 
Nous poserons en outre 
= 2LSEMEENE NET A I 
; 3 : sci Ma Me Ma 


De même, pour le pouvoir thermoélectrique 


tee (os A Z 
+ L'H(H.V) + M'VH? + N'T(V, H) 


avec (en unités c. g. s.) : 


> 
= 
Q - QUHAV 


Q L'hEree 
Te 

Q a MAR Tee 
ï <Q < x >2 ] 


PT To Te 
Mi Mae + Mo Ms + MaMi | 


k 
tie [6 De ln 
| (m1 + me + m3)? 
Mai Me + Mo M3 + Mama 
— Du) f(x) 
(m; + Mo SE m3)? 
M1 Mo + Mo M3 + m3 M); é 
k 
No [Du — 5 ui — 1) y; — 3] f{x) À 


{| 
M Ma + Mo Ms + MsMi 


a] 4 
M [Ou — 


s)] A 


£ 
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où les fonctions f(+) et g(+) ont été définies précé- 
demment, k, e, c, sont respectivement la constante 
de Boltzmann, la charge de l’électron et la vitesse 
de la lumière. 


V. Discussion. — Nous allons maintenant dis- 
cuter les résultats précédents en nous attachant au 
point de vue suivant : la connaissance des effets 
magnétoélectriques et thermomagnétoélectriques 
d’un cristal cubique, permet-elle d’atteindre la 
structure de bande de ce cristal ? 

Nous supposons donc que les résultats expéri- 
mentaux nous ont permis de déterminer les coef- 
ficients 65 61 LMN Q, Q, L'M' N'etnous essayons 
d’en déduire les inconnues A, +(8), m,, M», Ma 
Y1 Y2 Ya. Nous tiendrons compte du fait que les 
effets que l’on mesure sont d’autant plus faibles et 
plus difficiles à évaluer que l’ordre en Æ est plus 
élevé et que les mesures thermiques sont moins 
précises que les mesures électriques. 

La première difficulté provient de la présence de 
la fonction +(6) qui apparaît par des valeurs 
moyennes diverses. 

Il est a priori impossible de déterminer quoi que 
ce soit si l’on ne se fixe pas un modèle pour cette 
fonction. Nous prendrons + = +, 6" ce qui corres- 
pond à des lois que l’on obtient théoriquement 
quand on étudie le « scattering » des porteurs. La 
valeur de r dépend du processus de «scattering » ; 
ces valeurs sont généralement des entiers ou semi- 
entiers compris entre — 1/2 et 3/2. Avec cette 
hypothèse, les différentes valeurs moyennes de + 
se réduisent au produit d’une fonction connue der 
par une certaine puissance de vs. 

Remarquons ensuite que les paramètres y; ne 
figurent que par la combinaison Zu — 1) vi. 
Même si l’on arrive à déterminer les u; nous 
n’obtiendrons qu’une relation linéaire entre les 
3 paramètres indépendants y; ; il est donc impos- 
sible de déterminer rigoureusement la position des 
vallées dans le zone de Brillouin. Cette relation per- 
mettra cependant de rejeter les modèles pour les- 
quels elle n’est pas vérifiée, ce qui limite le nombre 
des essais que l’on est en droit de tenter. 

La connaissance des trois paramètres L, M, N, 
nous permet d'évaluer Z(ui — 1) y: et Eu. 

En effet, on a : 


SR RC MN 

DU CT EM UN 
38L+M+N 

Dit Me —55T EMEA 


Ce sont les seules relations que l’on peut obtenir 
sans faire d’hypothèses supplémentaires. 

Pour continuer nous procéderons ainsi : nous 
commençons par fixer r, il est alors possible de 
calculer les diverses valeurs moyennes de + et nous 
montrerons que l’on peut déduire ÎV, m;, Me, Mas 
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et To, à partir des quantités 6,6, LMN et 
Qo —u /kT. Pour tester la validité de l'hypothèse 
faite sur r, il reste à reporter les valeurs trouvées 
dans les équations donnant Q, L'M'N' et de véri- 
fier si les résultats obtenus coïncident avec les 
valeurs expérimentales. 

Montrons maintenant comment on peut faire la 
détermination des principaux paramètres quand on 
Hi r connu. Formons le rapport 6,6,/02, il 
vien 


So S1 


Nec < Dre 


Te + 22W) 


Co 3 


Zu est déjà connu, le rapport des valeurs moyennes 

de + qui ne dépend que de r est connu lui aussi. 

Cette équation permet donc de calculer M. 
Formons maintenant 6/0; 


MAT als ME UT 


Ma + Mo + Ma < Dm > 


cé  ( 1 1 1 } 


On 05 ME 


N étant connu, on en déduit la valeur d’une nou- 
velle fonction homogène par rapport aux trois 
masses m. Jointe à la connaissance de Zu, on peut 
en déduire les rapports des trois masses entre elles 
(ou celle des trois ;, ce qui est équivalent). 

Considérons maintenant le terme de pouvoir 
thermoélectrique Q, — uJkT. D’après les hypo- 
thèses faites ci-dessus, il est fonction de r, de N, 
de m, m, m, et du nombre V de vallées. Si l’on fait 
un choix pour V, on peut donc en déduire le pro- 
duit des trois masses puisque À et r sont connus. 
Comme nous connaissons déjà les rapports des 
masses entre elles, nous pouvons en déduire ces 
masses elles-mêmes. 

Enfin, en reportant les résultats précédents 
dans 6, par exemple, on peut calculer +4. Nous 
avons ainsi réussi à déterminer les différents para- 
mètres en n’utilisant que six des dix résultats expé- 
rimentaux disponibles ; nous pourrons donc tester 
la validité des hypothèses faites sur r et V à l’aide 
des quatre derniers résultats expérimentaux. Par 
approximations successives, il sera donc possible 
d'atteindre r et V et par suite tous les autres coeffi- 
cients. Remarquons que la détermination de V se 
trouve facilitée puisque V doit être un entier sous- 
multiple de 48. 


Conelusions. — La discussion que nous avons 
développée précédemment fait ressortir les points 
suivants : 

19 L'analyse des phénomènes de magnéto-résis- 
tance électrique et thermique des cristaux cubiques 
ne permet pas de déterminer la position des diffé- 
rentes vallées. Elle permet tout au plus de fixer une 
relation limitant le choix des modèles possibles et 
d'atteindre lenombre des vallées distinctes V. Il n’y 
a que dans les cas où ce nombre est faible que l’on 
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peut espérer préciser le modèle de la structure de 
bande. 

29 L'analyse des phénomènes de magnéto-résis- 
tance électrique seuls, ne permet la détermination 
de la densité de porteurs N et des rapports des 
masses effectives que dans la mesure où l’on choisit 
a priori un processus de «scattering » déterminé. 

30 Si l’on considère en outre les résultats de 
magnéto-résistance thermique, on peut théori- 
quement déterminer les paramètres suivants : 


N, Tor 7; Mi, Mo, M3 el 2 (u: red 1) Yi. | 


Les résultats expérimentaux sont évidemment 
beaucoup plus délicats à interpréter quand deux 
processus de scattering viennent se superposer. 


ANNEXE 1 
Calcul du vecteur V (H, E) 


Nous commençons par calculer le vecteur : 


W{(H,E) = 2 oU(oH, 6E) — X (oH.B oil) A 6E 
oc oeJe 
= » o8 As E 
oc 
où 


oÿ = oH.B oH 
A9 = 6A6o. 


En utilisant le fait que le sous-groupe Je est 
formé par éléments (1, 61, 62, 63) et de leurs 
opposés, la somme cherchée peut s’écrire : 

W(HE)=2 


G€(1,9:,03.08) 


o8 As E. 
En partant de 
© = Bis HŸ + By: H5 + B3s H3 
+ 2B19 H1 H3 + 2B13 Hi H3 + 2B33 Ha H3 


et 
A1 A2 A3 
A1 = |A; A2 A3 
As À 32 A3 


et en remarquant que : 
— $, se déduit de w?, en remplaçant H; par 
— Hi. 


— À,,8e déduit de À, en changeant le signe des 
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éléments contenant l’indice : une fois seulement, | 
on obtient immédiatement le vecteur W(H, ËE) | 
sous la forme | 


«€ 
où G(11) est la matrice suivante : 
Gi = 8Aü (Bis HŸ + Bs2 H3 + B2s Hi) 
Gis = Gi = 16 Ai Bi; Hi H; 
Le vecteur V(H, E) est alors donné par la somme : [! 


V(H,E) = > C—:W(CH, CE) 
ce x 
où X est le sous-groupe (1, C3, C£)) (1, Z) iso- 
morphe au groupe des permutation des 3 chiffres 
(152258) ï 
Le résultat précédent pour W(H, E) permet 
d’écrire : 
V(H,E) = Z C-:G(CH)CE = M(H)E 
cex 


où M(H) est la matrice qu’il nous faut évaluer et 
qui est définie par : À 


M(H) = À C—1 G(CH) C. 
€ 


Remarquons que C1 G(CH)C se déduit de G (H) 
par la transformation suivante : 

— H est remplacé par le vecteur CH dont les 
composantes se déduisant de celle de H par appli- 
cation de la permutation C. 

— Les éléments de G(H) sont remplacés par 
ceux que l’on obtient en effectuant sur les indices 
des lignes et des colonnes la permutation CT, 

En utilisant cette méthode de formation des 
matrices CT1 G(H)C, nous obtenons immédia- 
tement : 


Mis = (p — v) H? + (38v— 0) H} 
Mi; = M = 6 H;H; 
avec 
v= 82 Au Bu 
% 


8 Tr ATrB—8TraTr «—1 
16 2 Aÿ Bi. 


F © 
I 


I 


Manuscrit reçu le 18 novembre 1959. 
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ÉTUDE THERMOMAGNÉTIQUE DU SYSTÈME CoO — Ni0 


ParMIeNWESERRES, 


Laboratoire Pierre-Weiss, Institut de Physique, Strasbourg. 


Résumé. — Les oxydes CoO et NiO sont miscibles en toutes proportions. L'étude thermo- 
magnétique du système CoO-NiO confirme ce résultat : CoO seul, ou dissous dans NiO présente 
le même comportement. A haute température la variation de y est représentée par une loi de 
Weiss : 4(7 — 0) — 3,25. Les points de Curie paramagnétiques 0 et les points de transition 7% 
varient sensiblement linéairement avec la concentration. 


Abstract. — The oxides CoO and NiO are miscible in all proportions. 


The thermomagnetic study 


of the system CoO-NiO confirms the following result : CoO alone, or dissolved in NiO exhibits the : 


same behavior. 
X(T = 6) = 3,25. 
linearly with the concentration. 


Le protoxyde de cobalt, de structure cubique à 
faces centrées, peut former avec les protoxydes de 
magnésium, de nickel et de manganèse, de même 
structure, des solutions solides en toutes propor- 
tions [1]. J. Robin [2] dans ses recherches sur la 
structure et la stabilité des oxydes mixtes a établi 
les diagrammes d’équilibre des systèmesCoO—MgO 
et CoO—NiO ; le paramètre cristallin des oxydes 
mixtes varie linéairement en fonction de la compo- 
sition cationique depuis do — 4,252 À jusqu’à 
amso = 4,204 À ou axio — 4,168 À respectivement. 

L’étude  thermomagnétique du système 
CoO—M£gO [3], [4] a montré que CoO conserve 
dans toutes les solutions solides étudiées le même 
comportement qu’il présente à l’état pur: antiferro- 
magnétisme aux basses températures, maximum 
d’aimantation pour une température Te, puis, au- 
dessus, variation de la susceptibilité magnétique 
s’exprimant par une loi de Curie-Weiss : 
xu(T — 0) — 3,26 de même constante de Curie ; 
le point de Curie paramagnétique 0 variant linéai- 
rement avec la composition des solutions. 

Dans le composé 2Co0, NiO [4], où les deux 
composants sont antiferromagnétiques (alors que 
dans le système CoO—MgO un des composants est 
diamagnétique), l’oxyde de cobalt présentait aussi 
le même comportement. J. Robin a mis très aima- 
blement à ma disposition quelques échantillons des 
oxydes mixtes cobalt-nickel ayant servi à ses 
recherches, ce dont je le remercie vivement. Alors 
que le composé 2Co0, NiO étudié précédemment 
avait été préparé au four solaire par M. Foëx, les 
composés de Robin étaient obtenus par voie 
humide, en dissolvant des quantités connues de 
nitrates isomorphes Co(NO:):, 6H,0 et Ni(NO:)» 
6H,0, en évaporant les solutions à l’étude à 160 °C 
jusqu’à siccité ; le produit pulvérisé était ensuite 
calciné à température convenable. 


At high temperature the variation of x is represented by the Weiss law : 
The paramagnetic Curie points 0 and the transition points T, vary almost 


._ Six oxydes mixtes ont été étudiés, de compo- 
sitions atomiques comprises entre 20 et 95 % de 
cobalt. 


NiO. — Pour dégager le comportement de CoO 
dans ces solutions je me suis servi de la variation 
thermique obtenue avec un oxyde de nickel [5] 
très pur entre 77 0K et 1 308 °K. La susceptibilité 
magnétique de NiO croît des basses températures 
jusqu’à un maximum très étalé autour de 
Te — 650 °K, dans un intervalle de 40° environ ; 
au-dessus x décroît indéfiniment ; entre 800 et 
1 300 0K existe une loi de Weiss 


x(T + 3 000) — 3,66. 


Singer [6] trouve dans cette région une loi de Weiss 
X(T + 2 000) = 2,65 sur un monocristal. 


CoO. — L'étude du protoxyde de cobalt par 
Mie Henry La Blanchetais [7] donne pour le maxi- 
mum d’aimantation la température 7, — 292 °K 
et au-dessus de 460 0K une loi de Weiss : 


X(T + 280) = 3,0546 


pour un échantillon calciné sous vide à 820 0C ; au 
voisinage du maximum existait une sensibilité au 
champ, dans le sens antiferromagnétique. Avec un 
échantillon fondu sous vide au four solaire, obtenu 
à partir d’un nitrate de cobalt moins pur, en 
juin 1950, le même auteur trouvait une aimantation 
plus forte, la sensibilité au champ était dans le 
sens ferromagnétique : cet oxyde contenait une 
trace de nickel réduit ; seul le point T4 était sensi- 
blement le même (291 CK). 

M. Marc Foëx m'avait envoyé un autre échan- 
tillon fondu sous atmosphère oxydante au four 
solaire de Mont-Louis en 1954. Cet oxyde n’a pré- 
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Fic. 2. — Points de Curie paramagnétiques 0 et points de transition 7% en fonction de la concentration. 


senté aucune trace de ferromagnétisme et a donné 
les constantes suivantes : 


Te = 290 °K 
et au-dessus de 450 0K : 
Y(T + 283) — 3,20. 


Le maximum de susceptibilité (4, — 71,32.107$) 
est un peu supérieur à celui trouvé par Mile La 
Blanchetais avec l’échantillon le plus pur, mais 
inférieur à celui donné par l’échantillon contenant 
du nickel. La constante de Curie est également un 
peu plus élevée. 


Oxydes mixtes. — Les oxydes mixtes étudiés 
correspondaient aux concentrations de 95, 90, 73, 
50, 33 et 20 atomes de cobalt pour 100 atomes 
cobalt + nickel. La solution solide précédemment 
étudiée [4], 2Co0, NiO s'intègre dans cette série 
et correspond à 67 % de cobalt. 

La sensibilité au champ était extrêmement faible 
et toujours dans le sens de l’antiferromagnétisme : 
seule la solution à 33 % de cobalt présentait un 
faible ferromagnétisme, dû à une très petite quan- 
tité de nickel (décrochement de la courbe, 1/7, T, 
au point de Curie du nickel : 633 0K). Je ne retien- 
drai pour cette solution que la valeur du point de 


transition 7 — 429 0K, qui ne peut pas être faussée 
par cette trace de nickel. 

Les courbes de la figure 1 traduisent les résultats 
des mesures dans le plan (1/ym, 7) pour CoO et 
pour les 5 solutions étudiées, j’ai ajouté la courbe 
obtenue précédemment pour la solution à 67 % de 
cobalt : ces 7 courbes montrent pour CoO le même 
comportement quelle que soit la dilution : 1) anti- 
ferromagnétisme d’abord, à peu près constant, ; 
2) puis maximum de 7 au point de transition Te 
qui s'élève quand la proportion de cobalt diminue ; 
3) aux températures plus élevées, apparition, après 
un intervalle d’une centaine de degrés, d’une loi de 
Curie-Weiss de même constante de Curie pour les 
concentrations supérieures à 50 % de cobalt et 
à peine plus faibles pour les solutions les plus 
diluées. 

La température 7% du point de transition varie 
linéairement (fig. 2) avec la concentration. Bacon, 
Street et Tredgold [8] ont trouvé une relation 
analogue pour des solutions. Co0O—Mn0O. Par 
extrapolation, pour une concentration en cobalt 
nulle, on obtient une valeur de 7% à peine supé- 
rieure à 500 °K qui est loin du point de transition 
de NiO (650 °K). Street et Lewis [9] ont trouvé 
la température d’anomalie du module d’Young 
à 507 9K, inférieure à la température de 523 0K de 
l’anomalie de dilatation trouvé par M. Foëx [10], 
qui est aussi la valeur limite donnée par Richardson 

&1 
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294% 


et Milligan [11] pour la température de Néel nor- 
male de l’oxyde de nickel. 


TABLEAU I 
CONCENTRATION 
ATOM. Co /100 AT. TOI — 0 — O[Te C 
Ni + Co 
100 (CoO) 290 283 0,98 3,20 
95 305 265 0,87 3,20 
90 312 258 0,83 3,26 
73 347 238 0,69 9:21 
67 358 231 0,65 3,26 
50 395 189 0,48 2,78 
39 429 _ — -- 
20 472 90 (DIS 2,70 
0 (NiO) 650 — — — 


Mme G. Assayag et H. Bizette [12] ont montré 
que la chaleur spécifique de différentes solutions 
CoO—NiO présente un maximum à une tempé- 
rature 7% qui varie linéairement en fonction de la 
concentration moléculaire, s’extrapolant pour N10O 
vers la valeur 507 0K. Ces auteurs ont trouvé le 


N° 6 


maximum de chaleur spécifique aux températures 
suivantes : 


289,7 0K pour Co0 

3110K pour A0Ni0, 90C00 
3450K pour 25Ni0, 75Co0 
376,7 0K pour 40,1Ni0, 59,9Co0. 


Ces points se placent sur la droite donnée par 
les maxima d’aimantation (fig. 2). 

Les droites de Weiss obtenues aux hautes tempé- 
ratures donnent des points de Curie qui diminuent 
également régulièrement avec la concentration en 
cobalt. Le rapport — 0/7% sensiblement égal à 
lJunité pour CoO décroit à peu près comme la 
concentration en cobalt (tableau 1, colonne 3). 

En admettant l’additivité des propriétés ma- 
gnétiques dans les solutions solides du système 
CoO—NiO, nous avons pensé dégager le com- 
portement magnétique de CoO des mesures de 
susceptibilité faites sur ces solutions. Pour des 
concentrations supérieures à 50 % de CoO, le 
moment de l’ion Co?* dans CoO se conserve et est 
égal à 5,1 us (c — 3,25) comme il était dans les 
solutions du systène CoO—MgO. Pour des dilu- 
tions plus fortes le moment de Co?* baisse et 
n’atteint plus que 4,6 us pour 20 % de Co0. 


Manuscrit reçu le 30 mars 1960. 
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RÉSONANCES ET DÉDOUBLEMENTS ROTATIONNELS DU TYPE 4 
DANS LES MOLÉCULES A SYMÉTRIE AXIALE 


Par M. L. GRENIER-BESSON, 


Laboratoire de chimie physique de la Faculté des Sciences de Paris. 


Résumé. = Certains éléments matriciels, diagonaux en w, non diagonaux en 4 et Æ, de l’hamil- 
tonien transformé du second ordre h5 sont responsables des phénomènes de résonances et de dédou- 
blements rotationnels du type /. Nous avons calculé leurs expressions en fonction des paramètres 
moléculaires : moments d'inertie à l’équilibre, constantes de force du potentiel harmonique, cons- 
tantes d'interaction entre la vibration et la rotation, coefficients cubiques du développement de 
lénergie potentielle. Le calcul à été fait pour une molécule à symétrie axiale appartenant à un 
groupe de symétrie quelconque et appliqué à plusieurs modèles particuliers. 


Abstract. — Matrix elements of k£ (second-order transformed Hamiltonian), which are diagonal 
in v,and off-diagonalin /;and Æ are responsible for rotational /-type doubling and resonance. These 
matrix elements have been computed as functions of molecular constants : equilibrium moments 
of inertia, harmonic force constants, interaction coefficients between rotation and vibration, coef- 
ficients of the anharmonic cubic potential. Calculation has been performed for an axially sym- 
metric molecule belonging to a non-specified symmetry group, with application to various parti- 
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cular models. 


Introduction. — Il est bien connu que l’énergie 


de vibration-rotation d’une molécule à symétrie 


axiale peut être écrite sous la forme d’un dévelop- 
pement 


LES DPSE ES ARE ER ONE DE 


Ds. (e. +) + BE (I +1) — K9 + FE, 


(1) 
NC LUE. (16) 
he t 
Æ: gs 55’ Re 
LE — 2 Tss (os z) (o. =) . ay Ur 

BE it 


Le (v. | #) Her en (v.+ £) x 


S 
1) K2— Dr K* + Eÿ. 
(1e) 
Ce développement est lié au développement de 


l’hamiltonien par rapport aux coordonnées nor- 
males : 


H = H, + Hi + Hi +... (2) 


E, et E, sont données par les éléments dia- 
gonaux des opérateurs /7, et H,. L'énergie d'ordre 
zéro Æ, est la somme des énergies d’un rotateur 
rigide et d’un ensemble d’oscillateurs harmoniques. 
La correction du premier ordre Æ;, correspond à 
l'interaction de Coriolis entre le moment angulaire 
total et le moment angulaire interne dû à une 
vibration dégénérée. Pour calculer la correction du 


second ordre Æ,, il est commode d'utiliser la 
méthode de Shaffer, Nielsen et Thomas [1] et de 
considérer un hamiltonien transformé (1) : 


HER RN Nh Eee (3) 


FE, est alors obtenue à partir des éléments dia- 
vonaux en v, de l’opérateur A. Ces éléments peu- 
vent être, soit diagonaux, soit non diagonaux par 
rapport aux nombres quantiques de vibration /, et 
par rapport au nombre quantique de rotation K. 

La contribution des éléments diagonaux en /, et 
K (explicitée dans l’équation (1c)) correspond aux 
trois effets suivants : termes anharmoniques de 
l’énergie de vibration ; variation des moments 
d'inertie avec les nombres quantiques vibration- 
nels ; distorsion centrifuge. 

La contribution des éléments non diagonaux en 
1, et K, que nous avons désignée par £> se traduit 
par divers effets que l’on peut grouper en trois 
catégories (?) : résonances et dédoublements rota- 
tionnels du type { ; résonances et dédoublements 
vibrationnels du type / ; résonances et dédouble- 
ments du type Æ. 

Nous allons rappeler dans la première partie de 
cet article, la forme de l’hamiltonien h, et définir 
de façon précise à quels termes de Fhamiltonien 


(1) En l’absence de résonances liées à des dégénérescentes 
accidentelles (résonance de Fermi, résonance de Coriolis...) 
la correction du troisième ordre est obtenue à partir des 
éléments diagonaux en v, de l’opérateur A5 [2], [3]. En 
revanche, si l’on veut calculer la correction du quatrième 
ordre, il est nécessaire d'effectuer une deuxième transfor- 
mation de contact selon la méthode décrite par Amat et 
Nielsen [4], [5]. 

(2) La distinction qu’il convient de faire entre dédou- 
blement et résonance sera précisée dans la troisième partie, 


correspondent les différentes résonances que nous 
venons d’énumérer. 

Dans la deuxième partie, nous nous consacrerons 
à l’étude systématique des éléments matriciels 
contribuant aux résonances et dédoublements rota- 
tionnels du type /. 

La troisième partie rappellera les résultats rela- 
tifs au cas particulier des molécules linéaires. 

Enfin, dans la quatrième partie, au cours d’un 
bref aperçu historique, nous montrerons comment 
la théorie des résonances et dédoublements rota- 
tionnels du type / peut être utilisée pour l’interpré- 
tation des résultats expérimentaux. 


PREMIÈRE PARTIE 


Hamiltonien d’ordre deux. Définitions. 

Amat, Goldsmith et Nielsen [2] ont donné pour 
l’hamiltonien A; une expression particulièrement 
adaptée au calcul des corrections d’ordre deux et 
d’ordre supérieur : 


le = ZE, OU Y Pe Po Py Ps 
Lrstrnnne 
@ 3 + 
dE à 2 Core Pa Pb ne Yov Qa qu) 1% Pe (4) 
D ab 

É a<Ù 
+2 2 ® va Pa Pv Pe Pa 

L a<b<c 


1 
+ 2 > @ Ya 5 (ga 4o Pe + Pe Ga Go) Pa 
%X abc = 

a<b 


cd 
su > (2) ab = 
abcd 2 
a<b;c<d 


ns 


abcd 
a<b<c<d 


(Ga Go Pe Pa + Pe Pa Ga 4) 


(2) V'abca Ja 4v dc qd. 


La valeur des coefficients #$Y% intervenant 
dans cette expression a été calculée par M. L. Gre- 
nier-Besson, G. Amat et S. Maës [6]. 

Les symboles intervenant dans l’équation (4) ont 
la signification suivante : 

P, avec &« = %,y,z: composante du moment 
angulaire total suivant la direction d’un axe prin- 
cipal d’inertie. 

Gas 4v.… Coordonnées normales sans dimension. 

Pas Po . Moments conjugués des ga, 4 … 

Les indices a, b … sont utilisés ici pour simplifier 
l'écriture à la place des indices so, so’, utilisés par 
Nielsen. Rappelons que 6 = 1 ou 6 = 1, 2 suivant 
que la vibration est non dégénérée ou doublement 
dégénérée. Dans la suite de cet article, toutes les 
fois qu'il sera nécessaire, nous conviendrons de 
caractériser par l’indice n les symboles correspon- 
dant à une vibration non dégénérée (HD vanter 
par lindice £ les symboles correspondant à une 
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vibration doublement dégénérée (g,, Gus Pas Pt Vtb 
Bu il 
Nous ne nous intéressons qu'aux éléments matri-| 
ciels de À; qui sont diagonaux en v,(1),(?) et, soit dia- 
gonaux, soit non diagonaux en /, et en X. Les 
seuls termes susceptibles de fournir de tels éléments } 
sont les suivants : 


hi = >, AY Pa Ps Py Ps (5) 

apYy | 
hs — 2 > (S YÉbpi pee de Yuv Ga qu) Pa PR 

LAC DE à 

; 1 

hos = 2 (2) Ya o (Ga 4v Pe Pa + Pe Pa Ga 4) 

re ‘ 

5 » (2) Y'abcd Ga Go 4e Ga. 
UAITA à 


a<b<c<d 


Nous indiquons dans le tableau I la nature des | 
effets auxquels correspondent les éléments dia- | 
gonaux et non diagonaux des opérateurs ci-dessus : 


TABLEAU I 
HamIL- CONTRIBUTIONS CONTRIBUTIONS 
TONIEN DIAGONALES NON DIAGONALES 
EN L ET À 
ho Distorsion centrifuge. Résonance et dédou- 
blement du type 
Ke 
h39 Variation des mo- Résonance et dédou- 
ments d'inertie en blement rotation- 
fonction des nom- nels du type L. 
bres quantiques 
vibrationnels. 
h23 Termes anharmoni- Résonance et dédou- 


ques de l'énergie 
de vibration. 


blement vibration- 
nels du type L. 


(1) Tous les éléments matriciels de k2 sont en outre dia- 
ones par rapport aux nombres quantiques de rotation J | 
et I: 

(2) En l’absence de résonances accidentelles, les éléments 
non diagonaux en v, de l’opérateur h2 ne contribuent qu’à 
la correction du 4° ordre. C’est pour tenñir compte de leur 
contribution que l’on effectue la deuxième transformation 
de contact mentionnée dans la note (1) de la page 555. 


DEUXIÈME PARTIE 


Résonances et dédoublements rotationnels du 
type ? dans les molécules à symétrie axiale. — 
À. PRINCIPE DU CALCUL DES ÉLÉMENTS MATRI- 
CIELS DE 22 RESPONSABLES DES RÉSONANCES ET 
DÉDOUBLEMENTS ROTATIONNELS DU TYPE /. — Les 
éléments matriciels non nuls de l’opérateur 2, Ps 
peuvent être, soit diagonaux, soit de la forme 
(K|K + 1), (K|K + 2). Les éléments diagonaux 
en v, des opérateurs Gags Et PaPr peuvent être, soit 
totalement diagonaux, soit de la forme (4[4 + 2). 
Il en résulte que l'opérateur A2 possède des élé- 


Ë 


Ne 6 


_ ments matriciels diagonaux en v,, appartenant à 
lun des sept types suivants (1) : 
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Connaissant les éléments matriciels des opéra- 
teurs P, Pa, Ga Qr, Pa Po, On peut, sans difficulté, ca] - 
culer (voir Appendice 1) chacun des éléments (6a), 


(all + 2) (6a) (6b), (6c) … (6g) de h en fonction des nombres 
(K|K +1) (6b) quantiques de rotation J, K, des nombres quanti- 
(KLK à 9°) (6e) TUES de vibration v,, l, et de coefficients À, B … G 
: \ qui ne dépendent que des £Y, et des # Yo, 
(4, Klh +2, K +1) (64) Nous réserverons le symbole A::: aux divers 
(4, Ku+2,K+1) (6e) Coeflicients intervenant dans l'élément (64), le 
A Eds symbole B::: aux divers coefficients intervenant 
(&, Alu + 2, 2 (69) dans l’élément (6b), ete... comme l’indique le Ta- 
(Uk, Klh + 2, KT 2). (6g) bleau IT. 
TABLEAU II 


TABLEAU DES COEFFICIENTS DES ÉLÉMENTS MATRICIELS DE 92 
22 


ÉLÉMENTS DIAGONAUX PAR RAPPORT 
AUX NOMBRES QUANTIQUES VIBRATIONNELS 


L_ ÉLÉMENTS DIAGONAUX PAR RAPPORT 


AUX NOMBRES 
QUANTIQUES ROTATIONNELS 


ÉLÉMENTS EN 
ÉLÉMENTS EN 


( 1 
( 1 
ÉLÉMENTS EN (Æ|Æ +9 
Ù 


] 
ÉLÉMENTS EN (K|K — 2) 
On obtient par exemple pour l’élément (6/) : 

(4, KHlh;oll + 2, À + 2) = 
Fi heŸ [J(J + 1) — K(K + 1)] [J(J +1) 
— (4 +1) (A + 2)] (0 + & + 2) (0 + 4) }72 (7) 


— 


xx FA A tt LE +7 
(2) ht + (DT pts 


32TC 


t 
2 72 (2 yis = A Yan = #? 4 pot 


YU 2 YU lat. LU x7 
(0 Vus Æ 22 0) 8 4 2() Yes, 


à tit (LT DR vtt 

+ RON) à (Van, + FT" 
VIE Tr UV wtit (CU 

© Vtt — À @} tés) LE 95 (Co) Vaste 


| 


PUY volet CUYY o TU ylit 
+ HOT Mure MOT D} (8) 


où t, et t, désignent les deux composantes d’une 
vibration doublement dégénérée. 

Les coefficients #$ Y®? et #5 Y,, sont définis dans 
la référence [2] et explicités dans la table I de la 
référence [6]. 


(1) Seuls sont écrits dans les formules (6) les nombres 
quantiques par rapport auxquels l'élément n’est pas dia- 
gonal. 


) DURS Die 
) pipe 


ÉLÉMENTS EN ÉLÉMENTS EN 


(4|k + 2) (414 — 2) 

2 PT AA AA 
Va 
E DE 
ce. ct F° c 
+, CH “ x 
cc (6 F= 


Ils dépendent des paramètres moléculaires sui- 
vants : 

a) les constantes de force du potentiel harmoni- 
que (7) ; 

b) les moments d'inertie à l’équilibre (?) : 


12e == Lou et 1 


c) les coefficients d’interaction entre la vibra- 
tion et la rotation : coefficients d'interaction de 
Coriolis C£s 3; Coefficients de variation des 
moments et produits d’inertie en fonction des 


coordonnées normales aff, A'%,,: 


d) les coefficients £yss" des termes cubiques 
du développement du potentiel anharmonique. 

Pour calculer les coefficients À, B, C... il est 
intéressant de tenir compte des propriétés parti- 
culières de ces divers paramètres. 


4. Inventaire des coefficients Kkissos"s non nuls. — 
Quand nous avons commencé ce travail, aucune 
étude systématique des kiss n'avait été faite. 


: : ne al 2 
(1) Le potentiel harmonique s'écrit : V = 52 Xa Où: 
(2 
la coordonnée normale Q, étant liée à la coordonnée nor- 
male sans dimension 4, par la relation ga — (Aa?) "? Qu. 


. en , . x (2 x 
(2) Pour simplifier, on écrira /4 à la place de Jap, L, à Ta 
L 
place de J%. 


Sachant que les coeflicients 


Ktot'o't’o” Kionn’ Kttan 
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(RS 1 


Fes ou DEEE) 


ou t= 
COURTE ANT NL S) 


sont nécessairement nuls dans le cas des molécules 
linéaires, on admettait en général qu’il en était de 
même dans le cas des molécules à symétrie axiale. 
Cette hypothèse avait été faite notamment par 
Shaffer dans son étude de la molécule pyramidale 
XY, [7], puis par M. L. Grenier-Besson, G. Amat 
et S. Maës [6]. 

Certaines anomalies observées dans les résul- 
tats de la présente étude nous ont amenés à penser 
qu’une telle généralisation n’était pas valable et 
que le potentiel cubique ainsi défini était incom- 
plet pour les molécules à symétrie axiale. Cette 
remarque a suggéré à Henry et Amat [8] d’entre- 
pendre l’étude systématique des coefficients cubi- 
ques ko sos’ €t quartiques so sa" ss" du dévelop- 
pement dela fonction potentielle, étude qui a permis 
de montrer que seuls les coeflicients kw sont 
nuls (1) pour n'importe quelle molécule à symétrie 
axiale, les coefficients Æjsrotor et kinn pouvant 
être différents de zéro pour certains modèles (?). 

Les formules générales données dans la table I 
de la référence [6] et les considérations ci-dessus 
permettent de calculerles 8 Y,, et {$ Y‘?. Les expres- 
sions ainsi obtenues, que nous avons utilisées dans 
nos calculs, sont données dans l’Appendice IT. Elles 
sont valables pour toutes les molécules à symétrie 
axiale, de ce fait elles diffèrent des expressions que 
Fon trouve dans la Table IV de la référence [6] 
valables uniquement dans le cas des molécules 
linéaires. 


2. Propriétés des A. — Amat et Henry [10] 
ont établi les relations suivantes : 


ads 
AS RUDS "so As’ 


LAURE Se 94 
te = Le (9a) 
‘48 ‘48 a a _ + 4 aF> 
4 203 (9! “ie A + > IS FT Se (9b) 
5 als 
dont nous avons également tenu compte. 
9. PTOpDriétestdes QC he bn 


pour chaque groupe de symétrie auquel la configu- 
ration d'équilibre de Ia molécule est susceptible 


(1) Il importe de souligner que les coefficients ktonn, Sont 
nuls quelle que soit « Porientation » choisie pour les axes 
ll dans l’espace de configuration des coordonnées nor- 
males t (voir ci-après paragraphe 3). 

(?) Herman et Shaffer avaient déjà signalé que le 
coefficient Ægtn pouvait être non nul dans le cas d’ure 
molécule possédant un axe d'ordre 4 [9]. 
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d’appartenir, il existe un certain nombre de rela- 
tions : 

—_ entre les coefficients a non nuls 
Céoss” non nuls 


ksos’o's"s” non nuls. 


set }» » » 


— )] ) » 


Ces relations ont été établies par Henry et Amat 


[11],[8]. Elles dépendent de « l’orientation » rela- | 


tive du système d’axes x, y, z lié à la molécule et 
° 2 4 ‘ (à 
des différents systèmes d’axes #, {, l1 2, des es- 


paces de configuration des vibrations dégénérées. | 


Des simplifications peuvent être obtenues par un 
choix convenable de cette «orientation », certains 


des coefficients aË, CE so Ksos'o/s+” devenant alors 


nuls et certains des coefficients non nuls devenant 
égaux en valeur algébrique ou en valeur absolue. 


B) RÉSULTATS GÉNÉRAUX. — 1. Nous avons 
ainsi établi que les coefficients Af:*, Ak°%, B#%, Bt, 
C#, CE, DE sont nuls pour tous les types de molé- 
cules à symétrie axiale. 

Seuls les coefficients £!, F4, G£ pouvant être 
différents de zéro, seuls peuvent exister les éléments 
matriciels suivants : 


(ve, L,, K|h52|v4, l ÊE 7. K + 2) 


(104) 


= F£ he! [J(J +1) — A(K + 1)] [J(J +4) 
— (A +1) (A + 2)] (ve + 4 + 2) (+ k) } V8 
(vs, li, Klhiolts, k +9, X +) (10b) 
= GE he{ [J(J +1) — K(K + 1)] [J(J +1) 
— (A F1) (A 7 2)] (0 + 4 + 2) (04 F 4) Ve 
(ve, 4, Klhsolve, 4 +2, KF1)=Elhc(2K F1) (10c) 


&LJ(J +1) — AK + 1)] (0e & 4 + 2) (01 + &) }P2. 


Il est important d'ajouter que : l’élément (10a) 
existe pour toute molécule à symétrie axiale à condi- 
lion toutefois que la vibration t soit du type E, ; 
l'élément (10b) n'existe que si la molécule possède un 
axe d'ordre pair N — 4,6, 8... De plus la vibration 
dégénérée t doit être du type Æ, si l’axe est d’ordre 
6, E, si l’axe est d'ordre 8 .. E,  sil’axe est 

DE 1 
d'ordre N ; l'élément (10c) n'existe que pour les 
molécules possédant un axe d'ordre impair N — 3,5, 
7... et ne possédant pas de plan de symétrie 6. De 
plus la vibration dégénérée { doit être du type Æ, si 
l’axe est d’ordre 5, E, si l’axe est d’ordre 7... 
E y_, si l’axe est d’ordre N. 


2 
Ces résultats mettent en évidence le fait que la 
forme de la matrice { À; } dépend de l’ordre de 
l’axe de la molécule. Ils ont suggéré à G. Amat [12] 
l’énoncé d’une règle (1) simple exprimant la condi- 


0) Ce sont certaines divergences apparues entre les 
résultats de nos calculs relatifs aux coefficients ASP CT 
la règle d’Amat qui ont amené à reconsidérer le problème 
des coefficients Æsos’o’s’” de l'énergie potentielle dans le 
cas des molécules à symétrie axiale. 


PE N° 6 


l tion pour qu'un élément matriciel (v, LL K 
L + Al, & + Al, K + AK) soit différent de 
Cette condition s’écrit : 


hlv, 
ZÉTO. 


ENST 2 my A4 = pN (11) 

p désigne un nombre entier positif, négatif ou 
nul, N l’ordre de l’axe de symétrie de la molécule ; 
my — { si la vibration dégénérée est du type £,, 
mi = 2, si elle est du type E,.... 

Cette règle a été démontrée en utilisant la 
théorie des groupes et généralisée aux hamilto- 
niens d'ordre supérieur. 


Remarques. — a) La définition de ! soulève 
quelques difficultés lorsque l’axe de la molécule 
est un axe de rotation-réflexion. En fait nous dési- 
gnons par { l’ordre de la rotation d’un angle 27r/N 
par rapport à laquelle sont définis les types de 
symétrie À, B, E;,, E,, … Les conventions géné- 
ralement admises correspondent à : 

M=ropour les sroupes C Cr -Cov Ds Din: 
Da (n impair), S>, (n impair). 

— 2n pour les groupes D, (n pair), S2, (n pair). 

b) Dans le présent paragraphe et dans les 
expressions générales (12) et (13), nous écrivons À, 
B, E;, E», E3, sans préciser les indices g ou u,' ou ”, 
éventuellement nécessaires pour caractériser les 
types de symétrie de certains groupes particuliers. 
Nous écrivons de même À et B sans préciser les 
types À, ou 4,, B, ou B.. 

c) Rappelons enfin que pour les groupes pos- 
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sédant un axe de symétrie d’ordre 3 ou 4 on écrit 
E au lieu de Æ,. 

2. Nous donnons ci-dessous les expressions géné- 
rales des coefficients Æ, F4, G{ différents de zéro 
en fonction des constantes moléculaires. 

L'expression générale des coefficients #4 et 
Gi est la suivante : 


FE 
(t du type Æ:) a 
Tr À do Léa)? — (of 
CG  a9 3121091, En Os) 
( du type De | 327rc1% À 2 
SL | 
TE TE 34 + DA BE 
& Dai af] + 2 D A (cn) — (Ke 
nm M— À 
(de type Æ) 
Le vel A(s1+20 
+ 2iCtn Con] \ NE | > )ÿ/4 Cane, af, 
(de type 4) 
+ kan 08, + knnts 2 + Fat, at] 
1 : ; 
Gr do > oi L2Ætn an + chtityn DE 
n 4» 
(du type B) 
Legende Dekio an)] (2) 


d désigne le symbole de Kronecker : 34 — 0 
LE T ClouiE= L dl. 

Par ailleurs la signification des symboles X et 
X' et la valeur des facteurs Ô,, d,, e qui diffèrent 
suivant les groupes de symétrie et suivant que 
l’on a affaire à F4 ou à G' sont définies dans les 
tableaux ci-dessous. 


CogrriciEnT F4 (t Du TYPE E:) 


(1) Car af (x, B = x ou y) est nul. 


N GROUPES Ô; X X’ do € 
N = 3 Ca, Cha D;, Cr D;p, D, 56 CIDIOMMONOCEO WE 1 A E 0 + 1 
| Car D) Cause detre caca ee + 4 A Contribution nulle (1). en | +1 
Der l Or nue de RE RE (0 B Contribution nulle (1). +1 + 1 
CN, CNn, DN, CN»; D; Da, SN …. + 1 A Jos 0 + 1 
N impair 
SN A es Re eines le en eue le 0 B Va 0 + 1 
N = PTS SN he Hier ass ete lee vis 2 
2 
INOEN 
3 pair 
CoerFiciENT G£ (t DU TYPE ES : N PAIR) 
Fr Et 
} i i Le 1 — 1 
cer CAO AD Ce Dee enmeeste 0 F nee ue : - it 
<a NS EN Se se aie € suce + 1 ontributio 3 (P 
CN CNr  DN> CNo, DNRE Te, 0 B 12 0 — 1 
\o St rene + 1 A ne 0 ent 
NE") 4 
8, 10... }—— 
7 Pair 
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L'expression du coefficient El est la suivante : 


—?? { 3 pe ag 


1 3onclale Me 


Es 
( du type E— 


9 


4 


QT VA . XX VA XX XZ 
AURAS er 1 di At + ds )| 
34 + Aro 7 va 2 (90 7 
+ 4 2 Céér, Ces + Gr, Gars + Cire Gr 
V'=E£t À "+ À 
n3/2 1 + 20y ; 
a 71 t Î £ ar 
ps ‘ /4 [2 [4 
AAC PANNES PER MEET Cennrr, aëx 
y& : ya ; œ2 19 
+ tv, a, + kr, as + kit, a)] (13) 


Cette expression est très générale : la seule condi- 
tion restrictive est que t soit du type Æ,, (ce qui 


2 
implique AN impair). On remarque en particulier 
que la formule (13) ne contient aucun paramètre 
ajustable analogue aux facteurs Ô,, d,, e de la for- 
mule (12) et qu'aucune restriction n’est imposée 
au type de symétrie de £’ et {”. 

Toutefois lorsqu'on calcule Æf pour les difté- 
rents groupes de symétrie, on s'aperçoit : 1) que 
le premier crochet est nul pour un certain nombre 
de groupes ; 2) que le deuxième et le troisième 
crochet ne sont différents de zéro que si t’ et t” 
appartiennent à des types de symétrie bien déter- 
minés (voir Tableau IIT). 


TABLEAU III 


t 
COEFFICIENT E+. 


TyPrE TYPE DE 
GROUPES DE SYMÉTRIE TERME EN SYMÉTRIE 
ë aB 4% p 
DE { a JE a ao DE 
na 14 
C3; Cas D; E F 0 E LE 
Gs 70) Eg F 0 Es E, 
2 6» l’3d 
Du = 0 Ey E9 
BANC, 1DF Th ==") EE HT 
10: D5a [ Eso At) E 29 Fo 
l En EU Eu E9 
C3, Cv; D; E3 = 0 E3 E; 
S D À E3g F= 0 E 39 E;9 
914 L7d E nr 
l zu = Eu E;9 
etc. 


C) APPLICATION A QUELQUES MOLÉCULES PARTI- 
CULIÈRES. — À titre d’exemple nous donnerons ici 
; ; : 
l’expression des coefficients F4, Gi, Ef pour les 
modèles moléculaires suivants: X,Y, X,YZ,.., X 1e 
pyramidal, XY, plan, X,Y,. 
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1) Molécules X,Y (ammoniac) X3YZ (halogé- 
nures de méthyle), X,YZT (cyanure de méthyle); 
X,YZTU (méthyl acétylène).… 

Ces molécules appartiennent au groupe C:,. Elles 
possèdent des vibrations non dégénérées n du 
type À, et des vibrations dégénérées t du type Æ£. 
On a 
t ns Ë (ci)? 


= 2T 


JO TMRR CESSE 


— " 
32rcl> N°? A) M — À 


EH x 

a 31? 3k dt > ent Er 

— S D itits [3/4 . 4 MATE - 3/4 
M VE M 


Al 
Gi =0 (15) 
TL À 
t — }? dt, 
Ex es 172 |: 7 (16) 
SonCLe Le O1 æ& 
34 + y 
L 4 > _ = Gt, de, | 
DÉS RE | » 
PE a es F7 
RE AR en 1 Neil 
kT, | ati PUE 
Le coefficient F4 a été calculé par Niel- 


sen [13] (1). Son résultat et le nôtre diffèrent sur 
deux points, à savoir : 

la formule (14) contient un terme supplémen- 
taire en ak ; 

les coefficients donnés par les 2 formules sont de 
signes contraires. 

Nous remarquerons à ce sujet que le coefficient 
F' n’est autre que la constante qg° de dédoublement 
du type / (voir IIe partie, paragraphe D). Nous. 
avons déjà signalé dans une publication concernant 
les molécules linéaires [14] que le signe de la cons- 
tante g° devait être changé dans les formules de la 
référence [13]. 


2) Molécule XY, plane (trifluorure de bore). — 
Cette molécule appartient au groupe D;,. Elle pos- 
sède une vibration non dégénérée nr du type 45, 
une vibration non dégénérée n' du type Ai et 
deux vibrations dégénérées t et t’ du type E' on a : 


2 34 + À 
Ft. e | t ñ es | 
NM6rC TN ARE (Gén) 
n3/2 &e a 
tien Joke ee le 17 
GE LP + lui ‘e 
CEE 200 (18) 


l’élément orienteur des bases dégénérées étant un 
plan 6, ou un axe C,, intersection d’un plan 6, et 
du plan 6; [11]. 


(*) Notons que la formule-obtenue par Nielsen n’est pas 
valable pour une molécule quelconque appartenant au 
groupe C;. Tout comme notre formule (14) elle est applicable 
aux molécules X,Y, X,Y7Z..., Le calcul général développé 
dans le paragraphe précédent montre que ces formules 
doivent être modifiées si l’on considère les molécules Cy 
susceptibles de posséder des vibrations du type A4. 
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3) Molécule X Y, pyramidale. — Cette molécule 
appartient au groupe C4. Elle possède : deux 
vibrations non dégénérées du type À, ; trois vibra- 
tions non dégénérées des types B, ou B, (1) ; deux 
vibrations dégénérées du type £. 

Si l’on choisit un plan o, pour «orienter » les 
bases dégénérées on obtient : 

Ft. — h? pe à 2 3u + 2 (a) | (19) 


39rc12 ie 2 I % (AD X— Ân 


TX Hu 
x: 23/2 : 5 k dm S : Fu 
n htin 3/4 Il n titan 3/4 


2 
167% L An (Ba) n 


t LE »> 34 + Ân fo Ve 
Le 9 (& JE 20 
16rcl? ju LE À — Àn Ds ) 
2 34 + M (er | 
FR DIE pr 
(B) M — Ân à 
73/2 TES a? 
— Pa (À be La À un À 
1 61% (B;) ñn (B:) Ân 
Et = 0. (21) 


C’est J. de Heer [15] qui, le premier, a signalé 
l’existence de deux possibilités de dédoublement 
du type / dans les molécules XY,. Il a calculé les 
coefficients F4 et Gf. mais il est difficile de retenir 
ses résultats pour plusieurs raisons : la fonction 
potentielle utilisée par de Heer est incomplète ; 
des erreurs de signe se sont glissées dans le calcul ; 
l’orientation choisie pour les coordonnées dégé- 
nérées est peu commode. 


4) Molécule X,Y, (allène) — Molécule du 
groupe D.4, elle possède des vibrations non-dégé- 
nérées n des types À,, B,, B, et des vibrations 
dégénérées du type Æ£. On a, si on choisit pour élé- 
ment orienteur (?) un plan 64: 
OMR 

16Tcl% À (B,) M — }n 


» 34 + M (cn)? | 


(B,) À4 nr: Nr 
3/2 ne Us 
T Ë >) ET 7 2 Kutn = (22) 
161% (B3) in (B:) Ân 
t 7 3 (age)? + 3 À + ÀÂn ot }? 9 
G+ Er MT VOS 2 (és (23) 
32rcli N°12 L (4) À — Àn : 


23/2 
161 


xx ay 
9 » k An DE n An. 
2 np Mahn ga nm “ht 374 


(B:) nm (B;) \n 


E = 0. (24) 

(!) B, et B, caractérisent respectivement les vibrations 
symétriques et antisymétriques par rapport à un plan © 
passant par l’axe de la molécule : il y a deux vibrations du 
type PB, et une du type B;, si l’on choisit pour plan 6, un 
plan (XY,) contenant deux noyaux Y ; il y a deux vibra- 
tions du type BP, et une du type B,, si l’on choisit pour plan 
, l’un des plans bissecteurs des plans (XY;). AE 

(2) C'est-à-dire si les coordonnées normales dégénérées 
sont choisies de telle façon que g, et g, Soient respective- 
ment symétriques et antisymétriques vis-à-vis de la 
réflexion par rapport à un plan 6. 
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D) DÉDOUBLEMENTS ET RÉSONANGES. — D'une 
manière tout à fait générale, on donne le nom de 
«résonances essentielles » ({) aux effets liés à des 
éléments matriciels (l}A,]l’) qui, comme lies 
éléments (10), sont diagonaux par rapport aux 
nombres quantiques v,. De tels éléments relient 
deux niveaux qui possèdent à l’ordre zéro des 
énergies voisines. 

On peut montrer [16], [17] que l’effet lié à un 
élément (l'}k,|[[”) donné est d’autant plus impor- 
tant que la différence ([}A']T) —([”|H'|T") est 
plus petite. Deux cas sont à distinguer : 

1) si (T|A'|T) = (1/|Æ/]T"), on se trouve en pré- 
sence d’une dégénérescence essentielle qui est 
levée par l’élément (T|k,]T'). On dit qu'il y a 
« dédoublement » ; l’effet est alors maximum et un 
élément (l|A/[") contribue à l’énergie au second 
ordre, (fig. La). 


a ; b) 
) V— 
DEMO RES 


ere 


2) Nous réserverons le nom de «résonances » 
aux cas où il n’y a pas « dédoublement ». C’est-à- 
dire aux cas où ([|A'|1) —(1"}Æ4'|[") est petit 
mais différent de zéro. L’effet est alors beaucoup 
moins important et en général l’élément (T|h:|1"} 
contribue à l’énergie à un ordre supérieur au 
second (fig. 1b). 

Les éléments matriciels (T|k[T”) étudiés dans le 
présent article sont donc susceptibles de donner 
lieu aux effets suivants : 

1) Résonances et dédoublements rotationnels du 
type L «2,2 ». 

(C'est-à-dire : AI, = +2, AK = + 2);ils sont 
dus aux éléments (104) et existent pour une molé- 
cule à symétrie axiale quelconque. Il y a dédouble- 
ment pour 4 = + 1, À — + 1 et résonance dans 
les autres cas. 

2) Résonances et dédoublements rotationnels du 
lype L«2, — 2 ». 

Ils sont dus aux éléments (10b) et n’existent 
que pour les molécules possédant un axe d’ordre 
pair. Il y a dédoublement pour 4 = +1, K = + 4, 
et résonance dans les autres cas. 

3) Résonances rotationnelles du type L «2, — 1 ». 

Elles sont liées aux éléments (10c) et n’existent 


(1) On appelle par ailleurs «résonances accidentelles » 
les effets liés à des éléments matriciels (Tlh,|1”) reliant 
deux niveaux (vs, 0...) et (Ds + As, Ds + Ans.) possé- 
dant accidentellement des énergies voisines du fait des 
valeurs particulières de certaines constantes moléculaires. 
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que pour les molécules possédant un axe d'ordre 
impair. Ce dernier effet ne peut pas donner lieu à 
un dédoublement. 


TROISIÈME PARTIE 


Cas particulier des molécules linéaires. — Les 
molécules linéaires peuvent être considérées comme 
un cas particulier des molécules à symétrie axiale 
(N = co). Sayvetz a montré qu'il fallait alors 
écrire : DER 

t 


Il en résulte immédiatement que parmi les élé- 
ments matriciels (6) seul l’élément (6/) est suscep- 
tible d’exister ; il s’écrit 
(v, ALES Dee) 

— Ph he E [+ 1) — (4 & 1) (J(J +1) 
— (le 1) (le 2 2) (ue + de +2) (ve + 4) }V2. (25) 


On trouve : 


D 


VE 


(26) 


S9rele M | 


Le coefficient F4 est identique à la constante q° 
de dédoublement du type / (1) ; compte tenu des 
relations entre les GG, ([10] eq. 194”) d’une molé- 
cule linéaire, il peut facilement se mettre sous la 
forme habituelle [13], [14] 


B° ©} 
: ) 
nl 
20} ñ On — ©} 
Nous rappelons que B, — are 
Sr Cl, 
et NOT CAY. 


QUATRIÈME PARTIE 


Utilisation des résultats précédents pour l’inter- 
prétation des spectres de vibration-rotation et de 
rotation pure. — A) MOLÉCULES LINÉAIRES. — Le 
dédoublement du type / est un effet connu depuis 
longtemps. C’est Herzherg [18] qui le premier a 
noté, dans le cas des niveaux r(/, — + 1) des molé- 
cules linéaires, la similitude entre l'effet qu'il a 
nommé dédoublement du type let le dédoublement 
du type À des niveaux d’énergie électronique des 
molécules diatomiques ; cet effet est en réalité un 
dédoublement rotationnel « 2, 2 »tel que nous l'avons 
défini dans la deuxième partie ($ D). Nielsen et 
Shaffer [19] ont publié une étude complète de ce 


(*)} Notons que la valeur en nombre d’onde du dédouble- 
ment a été écrite suivant les cas Av = qJ(J +1) [18] 
ou Av = 2q° J(J + 4) (4 + 1) [13],[14]. 
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phénomène et calculé la constante g du dédouble- 
ment (1). 

Nielsen, Amat et Goldsmith [24] firent appel à la 
résonance rotationnelle « 2, 2 » pour expliquer cer- 
taines anomalies de la constante de distorsion cen- 
trifuge des molécules linéaires. Comme nous l’avons 
déjà signalé, cette résonance est d’un ordre de 
grandeur beaucoup plus faible que le dédoublement 
rotationnel « 2, 2 » et, par suite, ne se manifeste 
que par des effets expérimentaux beaucoup plus 
fins : (le dédoublement se traduit par une modifi- 
cation de la constante d’inertie B, la résonance par 
une modification de la.constante de distorsion 
centrifuge D). 

Enfin la distinction entre résonances et dédouble- 
ments vibrationnels et rotationnels du type / a été 
définie par Amat et Nielsen [22]. 


B) MOLÉCULES À SYMÉTRIE AXIALE. — La théo- 
rie du dédoublement du type ! «2, 2 » établie dans 
le cas des molécules linéaires fut étendue par Niel- 
sen [20] aux molécules à symétrie axiale. Toutefois 
d’un point de vue expérimental, si le dédoublement 
du type / «2, 2» des molécules linéaires est un 
effet facilement observable dans le spectre infra- 
rouge de ces molécules, il n’en est plus de même 
pour les molécules à symétrie axiale par suite de la 
complexité de leur spectre de vibration-rotation. 
Jusqu'à présent, c’est surtout dans le domaine des 
microondes (spectres de rotation pure) qu’il a été 
possible d’observer cet effet. C’est aussi dans ce 
même domaine que l’on a pu observer les effets dus 
à la résonance « 2, — 1 ». 

En ce qui concerne le dédoublement et la réso- 
nance «2, — 2 »,1l n’a pas encore été possible d’en 
vérifier expérimentalement l'existence, car la 
molécule Ia plus accessible parmi celles possédant 
un axe d'ordre quatre est l’allène qui, malheureu- 
sement, ne présente pas de spectre de rotation 
pure. 


Spectres microondes. — Le dédoublement du 
type ! (2, 2 » est un effet facilement observable 
dans les transitions AJ = +1, AK — 0, si la 
molécule à symétrie axiale se trouve dans un état 
vibrationnel excité (avec [4] = 1). 

Une transition de ce type a été observée par 
Kessler, Ring, Trambarulo et Gordy [25] dans le 
cas des molécules de cyanure et d’isocyanure de 
méthyle. Le spectre a été interprété qualitativement 
par Nielsen [26] en tenant compte dans le calcul 


() Bien d’autres travaux sont à signaler relatifs au 
dédoublement du type { des molécules linéaires. Citons en 
particulier les études faites : sur la variation de q en fonc- 
tion du nombre quantique J par Nielsen [23] ; 

sur le signe de qg par Amat, Grenier-Besson et Cum- 
mins [14]; 

Sur j’ellet simultané du dédoublement du type Let de la 
résonance de Fermi du gaz carbonique par de Heer et 
Nielsen [21]. 
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des niveaux d'énergie : du dédoublement du type / 
(2, 2 » ; de la résonance du type / « 2, 2 » ; de cer- 
tains termes diagonaux de l’hamiltonien du troi- 
sième ordre. 

La théorie de Nielsen a été ensuite utilisée systé- 
matiquement pour interpréter le spectre de rotation 
pure des molécules à symétrie ternaire dans un 
état vibrationnel excité. Nous citerons en parti- 
culer l’étude du spectre du méthyl-acétylène par 
Trambarulo et Gordy [27], du trifluoro méthy!- 
acétylène par Anderson, Trambarulo, Sheridan et 
Gordy [28]. 

Cette théorie permet de rendre compte de façon 
satisfaisante de l’aspect général du spectre mais les 
valeurs ainsi obtenues pour certains paramètres 
moléculaires sont inacceptables (on trouve en parti- 
culier des coefficients &. dont la valeur absolue 
est supérieure à 1). Une théorie plus complète vient 
d’être développée par Maës (1), qui inclut dans la 
matrice de l’énergie tous les éléments diagonaux de 
l’hamiltonien transformé du troisième ordre et 
aussi les éléments de résonance du type «2, — 1 » 
que nous avons calculés dans la seconde partie de 
cet article. 

Cette théorie permet de rendre compte quanti- 
tativement des effets observés dans les spectres du 
cyanure et de l’isocyanure de méthyle. Elle devrait 
permettre également de résoudre les difficultés 
rencontrées dans les spectres des autres molécules 
à symétrie ternaire. 


Spectres infrarouges. — Jusqu'ici le dédouble- 
ment du type / «2, 2 » n’a été observé dans l’infra- 
rouge que dans le cas où il se trouve accentué par 
l’existence simultanée d’une résonance de Coriolis. 
Ce phénomène a été nommé par Nielsen [30] « dé- 
doublement du type / géant » et mis en évidence 
expérimentalement dans le spectre de l’ammoniac 
par Garing, Nielsen et Narahari Rao [29]. 

Je tiens à remercier ici M. G. Amat qui a dirigé 
les diverses étapes de ce travail. Je remercie éga- 
lement M. L. Henry qui a bien voulu se charger de 
l’étude de la molécule pyramidale XY, afin que 
nous puissions comparer nos résultats à ceux obte- 
nus par J. de Heer. 


Note ajoutée à la correction des épreuves. — Divers 
dédoublements et résonances rotationnels du type / ont été 
mis en évidence par W. S. BexEpicr dans son étude du 
spectre infrarouge de Pammoniac : 7 

1) Le dédoublement du type / «2,2 » a été observé dans 
la bande v, + v, où il n’est pas accentué par une réso- 
nance accidentelle [31]. , 

2) La résonance « 2,2 » a été observée en même temps 
que le dédoublement dans la bande v, + v4(321. 

D’autres exemples de résonance « 2,2 » ont été signalés 
dans les bandes 2v9 et 2v? [33] 


(*) Ge sont les résultats développés dans la deuxième 
partie de cet article qui ont suggéré à Maes la possibilité 
d'existence d’une résonance du type 202, —1 » dans le 
spectre du cyanure et de l’isocyanure de méthyle. 
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3) La résonance « 2,-1 », dont la présente étude a révélé 
l'existence a également été utilisée dans l’interprétation 
des spectres [34]. 

4) Le dédoublement « 2,-4 » (effet lié à l’hamiltonien 
du 3° ordre) a été observé pour les niveaux Æ — 2 [31]. 


Appendice I. — Pour calculer ses éléments ma- 
triciels, il est commode d’exprimer h;, en fonction 
des opérateurs vibrationnels V et rotationnels R 
suivants : 


: 1e De 
(n) re A EN EC 
} 9 (oi k ») 
po il pe Fe pa 
û 7 à Th h° le hi? 
7(l 1 2 pt, 2 Pt, , Pt Pt 
va = 9 (à 72 dt F ir (x Ît ie Ps Pu) 
; P} 
FU : (A1) 


Fate LP Re EMIPAP AP Pal, 


Les opérateurs PV, V, RF, R£ sont diago- 


naux. 
V®, ne possède que des éléments (4]4 + 2) et 
Riax ne possède que des éléments (X|X + AK). 
C’est ainsi que l’on a : 
(on| VU In) = 0n + (1/2) 
Us b| ph lu Lk) = (vw +1) 
ve Li] Vive, h + 2) = [{ve + 4 + 2) (ue + h)172 
(A2) 


( 
( 
(KIA |x) = x? 
( 
( 
( 


Mt ER TI 


RuolK + 2) = | 
(J(J HA) (RE 1) (f 22) 


On peut écrire : 


be % l & T0) Un) 
. = » 2 { A“ po, I. 4 + A+ Vi Ro 
ce +— nt 


+ BP Vu Rai + BE VO Ra 


+ CP Vin Rio + CP VS Rats 
+ DD V2 Ru (A3) 
+ E0 VD RE 


+ FŸ VE Res 


+ CT) PE Ro Î — opérateurs diagonaux. 
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Qu peut facilement Caen je coefficients À:: ne 8 yo ñ 2 CR 
B:::,ete.. en fonction des #5 Y, #5 Y* en na ne #4 )  Lolhonése Xa — in 
les équations (A1) et les deux relations suivantes : 0a #0} 


Dre 


mM + nN = (M + N) 


ee mit, (A4) 


2 


mM +- 


pos de #4 en fonction des $ Y® e 
GY a à été donnée dans le texte, nous Se 
ici l'expression des deux autres coefficients non 
nuls Ef et Gi. 


t 24 yrtat a2 
Ex = — ratée + rt Er (A5) 
tit 28 vtt 
D GT — ES ar, — 2 GT" 
1 À 24 yat UPA 

& US ur + OT — BY, 

V8 yrtlot A Ve xrtit 
Dre es LA re) 

td LE vrtit x 
GE ren + 22 (0 PE — 0) Fit, 


39Tc 


72 & yhts Li Yet RO we lat 


VU ytot: au 
+ D Ven + OT (Er, 


+ re) LUE Yes, + a y (A6) 


YU YY ylit) — 0: 
— D Val, — 2 GT) + QUE Yes, 


Fe rlat € y 
+ RE GNT — Ven, — 2 y)! 


Appendice II.— Coefficients °:, Y::: 
tonien As 


de l’hamil- 
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Ces formules sont valables aussi bien pour une 
vibration $, — nr non dégénérée que pour une vibra- 
tion $, — 8, — t dégénérée. On fera toutefois les 
2 remarques suivantes : les coefficients (A8) et 
(AÏO) ne peuvent exister que si la vibration 
Sa — Sy est dégénérée on posera alors s, = #, 
S, — le ; les coefficients 44, sont nuls dans (A7) et 
(A9) si la vibration s, — n est non dégénérée. 


Manuscrit reçu le 5 février 1960, 
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LETTRES À LA RÉDACTION 


MISE EN ÉVIDENCE D'UN PROCESSUS 
DE CONTAMINATION ISOTOPIQUE 
DANS LES 
SÉPARATEURS ÉLECTROMAGNÉTIQUES 
D'ISOTOPES 


Par R. BERNAS, 
Laboratoire de Physique Nucléaire. 
Faculté des Sciences, Orsay. 
et C. CASSIGNOL, 
Commissariat à l'Énergie Atomique, Saclay. 
L'expérience décrite dans cette note à permis de 


mettre en évidence un phénomène responsable de la 
contamination d’un isotope donné par l’isotope de 


46 


masse immédiatement supérieure, lors de la séparation 
électromagnétique des isotopes. 

Le dispositif expérimental utilisé était le séparateur 
à double déviation magnétique réalisé par l’un de nous 
au Laboratoire de Physique Nucléaire de la Faculté des 
Sciences à Orsay [1] (fig. 1). 

Cet appareil est constitué par la mise en cascade 
de deux étages de déviation magnétique. Le premier 
étage est un secteur symétrique à 600 en champ uni- 
forme, identique au Séparateur de Saclay dans sa 
version primitive [2], suivi d’un analyseur inhomo- 
gène du type Siegbahn et Svartholm à indice 0,5 et à 
déviation de 1800 ; les points images et objet de cet 
analyseur étant à l’extérieur du champ. L'image de la 
fente de la source d’ions, donnée par le premier étage, 
pour une masse déterminée, constitue l’objet pour le 
second étage. Une plaque métallique portant une 


( 
FD 


_ 
\\\ EE 


ER = Schéma des principaux éléments du séparateur à double déviation m 
a 7e . # " 11: Û à . x . 
1 Source d'ions. — 4 Électrode accélératrice. — 5-8 Enceinte à vide et canalis 
Fa pare à Se : AR : 
PE HA nAqNe à 60 (champ homogène). — 12 Diaphragme intermédiaire sit 
du 197 analyseur el permettant la sélection d’un faisceau isotopique. — 18 
à champ inhomogène d’indice 0,5. — 20 Collecteur final. 


agnélique utilisé dans cette expérience. 
alion de pompage. — 9 Ana- 
ué dans le plan de focalisation 
Analyseur magnétique semi-circulaire 
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fente convenablement taillée sépare les deux étages. 

Elle autorise le passage d’un isotope déterminé dans 

le deuxième étage, à exclusion de tous les autres. 
L'expérience a porté sur le Cadmium. Le premier 


- étage étant réglé sur un isotope donné, le second étage, 


par balayage du champ magnétique et enregistrement 
du courant fournissait un spectre de masse du faisceau 
admis. On avait ainsi, directement, mais limité aux 
composantes chargées, le spectre de la contamination 
isotopique du premier étage fonctionnant seul. 
Dans les conditions normales de fonctionnement 
(p = 1.105 mmH£g) on enregistre, si le premier étage 
est aligné sur l’un des isotopes 110, 111, 112 ou 113 un 
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spectre analogue à celui de la figure 2 (relevé pour le 
113Cd). Le pic satellite n'existe pas si l’on aligne le 
premier étage sur les isotopes 114 ou 116. La hauteur 
du pic principal et celle du satellite sont proportion- 
nelles aux abondances isotopiques de l’isotope admis 
et de l’isotope de masse immédiatement supérieure. 
De plus si l’on aligne le premier étage par interpolation 
sur l’un des isotopes absents 109 ou 115, on peut cons- 
tater la présence à l’extrémité du second étage de pics 
dus respectivement aux ions du Cadmium 110 ou 116. 

Si l’on fait croître la pression dans le premier étage 
(en injectant de l’argon) jusqu’à environ 5.10—5 mm 
de mercure l'intensité du pic principal et celle du pic 
satellite diminuent et on voit apparaître et croître deux 
autres pics (fig. 2). L'identification effectuée comme 
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précédemment montre qu'il s’agit de faisceaux d'ions 
de Cadmium 112 et 114 que l’on doit attribuer à un 
phénomène de diffusion sans changement de charge 
dans le premier étage. 

On constate, de plus, que la place de ces pics corres- 
pond très exactement à la dispersion du deuxième 
étage, alors que le premier pic satellite de Cadmium 114 
est à une distance supérieure, 

Cette position apparemment aberrante est attri- 
buée à une incidence du faisceau de 11Cd nettement 
différente de celle des autres faisceaux transmis dans 
le deuxième étage, associée à la présence d’un défaut 
d'ouverture dans cet étage au moment de ces expé- 
riences. La traînée de Cadmium 113 vers les pics 114 
est due à ce défaut d'ouverture. 

Si l’on intercepte les faisceaux d'ions au voisinage 
du diaphragme intermédiaire par des plaques métal- 


HTC 


liques À ou B disposées comme sur la figure 4, on 
aboutit au résultat paradoxal au premier abord, qu’en 
approchant A progressivement, jusqu'à la limite 
(visible) du faisceau de H3Cd on fait disparaître du 
spectre de la figure 3 le Cd le plus éloigné sans 
qu'aucune atténuation des autres pics soit détectable. 
Si l’on pénètre plus avant dans le faisceau de H#Cd les 
autres pics diminuent à leur tour et disparaissent. 

Les phénomènes sont identiques si la plaque B est 
approchée du faisceau de l’autre côté du diaphragme. 

L'intérêt pratique que cette observation présente 
pour l’amélioration de la pureté des isotopes séparés 
par voie électromagnétique est certain. 

Dans le cas des appareils à un seul étage, on placera 
le diaphragme d'ouverture en avant du collecteur et du 
côté des isotopes plus légers. 

L'application aux appareils composites est aisée, Si 
elle n’est pas indispensable pour les appareils à double 
déviation magnétique dont le deuxième étage redis- 
perse les masses elle peut toutefois apporter une amé- 
lioration. Par contre, dans le cas du séparateur de 
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Saclay qui comporte un deuxième étage électrosta- 
tique [3] inapte à redisperser les masses, le contrôle de 
l'ouverture du faisceau admis dans le deuxième étage 
doit permettre des progrès sensibles. 

Des résultats obtenus récemment, lors de la sépara- 
tion des isotopes du mercure, confirment cette prévi- 
sion. 

Une explication de ce phénomène a été proposée par 
l’un de nous [4]. 


Lettre reçue le 30 mars 1960. 
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ÉNERGIE DE VIBRATION-ROTATION 
DES MOLÉCULES POLYATOMIQUES. 
NOTE SUR LES TABLES DES COFFFICIENTS 
DE L'HAMILTONIEN TRANSFORMÉ 
DU SECOND ORDRE 


Par M. L. GRENIER-BESSON, G. AMAT 
et S. MAES, 


Laboratoire de Chimie Physique 
de la Faculté des Sciences de Paris. 


Nous avons donné dans une publication anté- 
rieure [{]les valeurs des coefficients #f- y9- af Ted. 
Ut figurant dans les hamiltoniens du second 
Ordre: lo, he, ds 

Dans les Tables I, IT et III de la référence [1], ces 
coefficients sont exprimés sous une forme générale, 
valable pour n'importe quelle molécule. Mais la 
Table IV a été établie en faisant certaines hypothèses 
sur la nature des coefficients £ape non nuls du potentiel 
cubique. Ces hypothèses, explicitées dans le Tableau A 
de [1], sont certainement valables pour les molécules 
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linéaires et les molécules asymétriques. Elles étaient 
également en accord avec la forme généralement uti-e 
lisée pour le potentiel cubique des molécules à symétrie 
ternaire (voir par exemple [2] pour le modèle XY;). 
Nous avions admis que le Tableau A et, par conséquent, 
la Table IV étaient valables pour toutes les molé- 
cules à l'exception des molécules à symétrie sphérique 
(modèles XY, tétraédrique et XY, octaédrique (°)). 
Cependant, certaines difficultés rencontrées récem- 
ment par M. L. Grenier-Besson dans l’étude des réso- 
nances rotationnelles du type / [4] nous ont amenés à 
penser que la question du potentiel cubique des molé- 
cules à symétrie axiale devait être reprise et faire 
l’objet d’une étude générale pour les différents groupes 


de symétrie. Cette étude a été effectuée par L. Henry | 


en collaboration avec G. Amat [5] : elle a montré que 
des termes avaient été omis dans les potentiels cubiques 
utilisés précédemment par divers auteurs pour les 
molécules à symétrie ternaire. 

En ce qui concerne la référence [1], la Table VI est 
valable pour les molécules linéaires et pour les molé- 
cules asymétriques (2). Le texte doit donc être modifié 
de la façon suivante : 

— Page 787, première colonne, 5 lignes avant la 
fin, supprimer : « molécule à symétrie axiale ou une ». 

— Page 787, deuxième colonne, 10 lignes avant la 
fin, lire : « molécules à symétrie axiale ou sphérique ». 

(1) Pour ces molécules voir [3]. 

(2) Une Table analogue à la Table IV peut être établie 
pour les molécules à symétrie axiale ou sphérique en par- 
tant des Tables I, II, III et en utilisant pour les kan les 
résultats donnés dans les références [5] ou [3]. 


Lettre reçue le 28 mai 1960. 
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